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Abstract 

An introductory survey of the theoretical ideas and calculations and the ex- 
perimental results which depart from Landau Fermi-hquids is presented. Common 
themes and possible routes to the singularities leading to the breakdown of Landau 
Fermi liquids are categorized following an elementary discussion of the theory. Sol- 
uble examples of Singular Fermi liquids include models of impurities in metals with 
special symmetries and one-dimensional interacting fermions. A review of these is 
followed by a discussion of Singular Fermi liquids in a wide variety of experimental 
situations and theoretical models. These include the effects of low-energy collective 
fluctuations, gauge fields due either to symmetries in the hamiltonian or possible 
dynamically generated symmetries, fluctuations around quantum critical points, the 
normal state of high temperature superconductors and the two-dimensional metallic 
state. For the last three systems, the principal experimental results are summarized 
and the outstanding theoretical issues highlighted. 



Contents 



1 Introduction 



1.1 Aim and scope of this paper 



1.2 Outline of the paper 



4 
4 
7 



2 Landau's Fermi-liquid 



2.1 Essentials of Landau Fermi-liquids 



2.2 Landau Fermi-liquid and the wave function renormalization Z 



2.3 Understanding microscopically why Fermi-liquid Theory works 



11 
17 



^ Present and permanent address 



Preprint submitted to Elsevier Preprint 



1 February 2008 



2.4 Principles of the Microscopic Derivation of Landau Theory 



2.5 Modern derivations 



2.6 Routes to Breakdown of Landau Theory 



3 Local Fermi-Liquids & Local Singular Fermi-Liquids 



3.1 The Kondo Problem 



3.2 Fermi-liquid Phenemenology for the Kondo problem 



3.3 Ferromagnetic Kondo problem and the anisotropic Kondo problem 



3.4 Orthogonality Catastrophe 



3.5 X-ray Edge singularities 



3.6 A Spinless Model with Finite Range Interactions 



3.7 A model for Mixed- Valence Impurity 



3.8 Multi-channel Kondo problem 



3.9 The Two-Kondo-Impurities Problem 



4 SFL behavior for interacting fermions in one dimension 



4.1 The One Dimensional Electron Ga^ 



4.2 The Tomonaga-Luttinger Model 



4.3 One-particle Spectral functions 



4.4 Thermodynamics 



4.5 Correlation Function^ 



4.6 The Luther-Emery Model 



4.7 Spin charge separation 



4.8 Spin-charge Separation in more than one-dimension? 



4.9 Recoil and the Orthogonality Catastrophe in Id and highei 



4.10 Coupled One-dimensional Chains 



4.11 Experimental observations of one-dimensional Luttinger liquid| 



behavior 



22 
27 
29 
34 
34 
38 
39 
40 
41 
43 
45 
47 
51 
55 
57 
59 
62 
63 
64 
65 
67 
69 
72 
76 

76 



2 



5 Singular Fermi-liquid behavior due to gauge fields 



5.1 SFL behavior due to coupling to the electromagnetic field 



5.2 Generalized gauge theories 



6 Quantum Critical Points in fermionic systems 



6.1 Quantum critical points in ferromagnets, antiferromagnets, and 



charge density waves 



6.2 Quantum critical scaling 



6.3 Experimental Examples of SFL due to Quantum-Criticality: Open 



Theoretical Problems 



6.4 Special complications in heavy fermion physics| 



6.5 Effects of impurities on Quantum Critical Point; 



7 The High-Tc Problem in the Copper-Oxide Based Compounds] 



7.1 Some basic features of the high-Tc materials 

7.2 Marginal Fermi Liquid behavior of the normal state) 



7.3 General Requirements on a Microscopic Theory 



7.4 Microscopic Theory 



8 The Metallic State in Two-Dimensions 



8.1 The two-dimensional Electron Gas 



79 
79 
82 
84 

85 



96 
101 
103 
104 
105 
108 
115 
117 
123 
124 



.2 Noninteracting Disordered Electrons: Scaling Theory of Localization 126 

130 



8.3 Interactions in Disordered Electrons 



8.4 Finkelstein Theory 



8.5 Compressibility, Screening length and a Mechanism for Metal- 



Insulator Transition 



8.6 Experiments 



8.7 Discussion of the Experiments in light of the theory of Interacting 



135 



137 



138 



144 



Disordered Electrons 



.8 Phase Diagram and Concluding Remarks 



150 



3 



9 Acknowledgements 



References 



152 
152 



1 Introduction 



1.1 Aim and scope of this paper 



In the last two decades a variety of metals have been discovered which dis- 
play thermodynamic and transport properties at low temperatures which are 
fundamentally different from those of the usual metallic systems which are 
well described by the Landau Fermi-liquid theory. They have often been re- 
ferred to as Non-Fermi-liquids. A fundamental characteristic of such systems 
is that the low-energy properties in a wide range of their phase diagram are 
dominated by singularities as a function of energy and temperature. Since 
these problems still relate to a liquid state of fermions and since it is not a 
good practice to name things after what they are not, we prefer to call them 
Singular Fermi-liquids (SFL). 

The basic notions of Fermi-liquid theory have actually been with us at an 
intuitive level since the time of Sommerfeld: He showed that the linear low 
temperature specific heat behavior of metals as well as their asymptotic low 
temperature resisitivity and optical conductivity could be understood by as- 
suming that the electrons in a metal could be thought of as a gas of non- 
interacting fermions, i.e., in terms of quantum mechanical particles which do 
not have any direct interaction but which do obey Fermi statistics. Mean- 
while Pauli calculated that the paramagnetic susceptibility of non- interacting 
electrons is independent of temperature, also in accord with experiments in 
metals. At the same time it was understood, at least since the work of Bloch 
and Wigner, that the interaction energies of the electrons in the metallic range 
of densities are not small compared to the kinetic energy. The rationalization 
for the qualitative success of the non-interacting model was provided in a mas- 
terly pair of papers by Landau [144,145] who initially was concerned with the 
properties of liquid ^He. This work introduced a new way of thinking about the 
properties of interacting systems which is a cornerstone of our understanding 
of condensed matter physics. The notion of quasiparticles and elementary ex- 
citations and the methodology of asking useful questions about the low-energy 
excitations of the system based on concepts of symmetry, without worrying 
about the myriad unnecessary details, is epitomized in Landau's phenomen- 
logical theory of Fermi-liquids. The microscopic derivation of the theory was 
also soon developed. 
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Our perspective on Fermi- liquids has changed significantly in the last two 
decades or so. This is due both to changes in our theoretical perspective, and 
due to the experimental developments: on the experimental side, new materi- 
als have been found which exhibit Fermi-liquid behavior in the temperature 
dependence of their low temperature properties with the coefficients often a 
factor of order 10'^ different from the non-interacting electron values. These 
observations dramatically illustrate the power and range of validity of the 
Fermi-liquid ideas. On the other hand, new materials have been discovered 
whose properties are qualitatively different from the predictions of Fermi- 
liquid theory (FLT). The most prominently discussed of these materials are 
the normal phase of high-temperature superconducting materials for a range of 
compositions near their highest Tc. Almost every idea discussed in this review 
has been used to understand the high-T^ problem, but there is no consensus 
yet on the solution. 

It has of course been known for a long time that FLT breaks down in the 
fiuctuation regime of classical phase transitions. This breakdown happens in 
a more substantial region of the phase diagram around the quantum critical 
point (QCP) where the transition temperature tends to zero as a function of 
some parameter, see Fig. 1. This phenomenon has been extensively investi- 
gated for a wide variety of magnetic transitions in metals where the transition 
temperature can be tuned through application of pressure or by varying the 
electronic density through alloying. Heavy fermion with their close competi- 
tion between states of magnetic order with localized moments and itinerant 
states due to Kondo-effects appear particularly prone to such QCP's. Equally 
interesting are questions having to do with the change in properties due to 
impurities in systems which are near a QCP in the pure limit. 

The density-density correlations of itinerant disordered electrons at long wave- 
lengths and low energies must have a diffusive form. In two-dimensions this 
leads to logarithmic singularities in the effective interactions when the inter- 
actions are treated perturbatively. The problem of finding the ground state 
and low-lying excitations in this situation is unsolved. On the experimental 
side the discovery of the metal-insulator transition in two dimensions and the 
unusual properties observed in the metallic state makes this an important 
problem to resolve. 

The one-dimensional electron gas reveals logarithmic singularities in the effec- 
tive interactions even in a second-order perturbation calculation. A variety of 
mathematical techniques have been used to solve a whole class of interacting 
one-dimensional problems and one now knows the essentials of the correlation 
functions even in the most general case. An important issue is whether and 
how this knowledge can be used in higher dimensions. 

The solution of the Kondo problem and the realization that its low temper- 
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parameter 



Fig. 1. Schematic phase diagram near a Quantum Critical Point. The parameter 

along the x-axis can be quite general, like the pressure or a ratio of coupling con- 
stants. Whenever the critical temperature vanishes, a QCP, indicated with a dot 
in the figure, is encountered. In the vicinity of such a point quantum mechanical, 
zero-point fluctuations become very important. However when Tc is finite, critical 
slowing down implies that the relevant frequency scale goes as cj ~ \Tc — T\'^^ , 
dwarfing quantum effects; the standard classical critical methodology then applies. 
An example of a phase diagram of this type for MnSi is shown in Fig. 34 below. 

aturc properties may be discussed in the language of FLT has led in turn 
to the formulation and solution of impurity models with singular low energy 
properties. Such models have a QCP for a particular relation between the 
coupling constants; in some examples they exhibit a quantum critical line. 
The thermodynamic and transport properties around such critical points or 
lines are those of local singular Fermi-liquids. Although the direct experimen- 
tal relevance of such models (as of one-dimensional models) to experiments 
is often questionable, these models, being soluble, can be quite instructive in 
helping to understand the conditions necessary for the breakdown of FLT and 
associated quasiparticle concepts. The knowledge from zero-dimensional and 
one-dimensional problems must nevertheless be used with care. 

A problem which we do not discuss but which belongs in the study of SFL's 
is the Quantum Hall Effect problem. The massive degeneracies of two-dimen- 
sional electrons in a magnetic-field leads to spectacular new properties and 
involves new fractional quantum numbers. The essentials of this problem were 
solved following Laughlin's inspired variational calculation. The principal rea- 
son for the omission is firstly that excellent papers reviewing the developments 
are available [199,71,107] and secondly that the methodology used in this prob- 
lem is in general distinct from those for discussing the other SFL's which have 
a certain unity. We will however have occasions to refer to aspects of the quan- 
tum Hall effect problem often. Especially interesting from our point of view is 
the weakly singular Fermi- liquid behavior predicted in the i> = 1/2 Quantum 
Hall Effect [113]. 
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One of the aspects that we want to bring to the foreground in this review is the 
fact that SFL's all have in common some fundamental fcatTircs which can be 
stated usefully in several different ways, (i) they have degenerate ground states 
to within an energy of order ksT. This degeneracy is not due to static external 
potentials or constraints as in, for example the spin-glass problem, but degen- 
eracies which are dynamically generated, (^iij Such degeneracies inevitably lead 
to a breakdown of pcrturbative calculations because they generate infra-red 
singularities in the correlation functions, (in) If a bare particle or hole is added 
to the system, it is attended by a divergent number of low energy particle- 
hole pairs, so that the one-to-one correspondence between the one-particle 
excitation of the interacting problem and those of the non-interacting prob- 
lem, which is the basis for FLT, breaks down. 

On the theoretical side, one may now view Fermi-liquid theory as a forerunner 
of the Renormalisation Group ideas. The renormalisation group has led to a 
sophisticated understanding of singularities in the collective behavior of many- 
particle systems. Therefore these methods have an important role to play in 
understanding the breakdown of FLT. 

The aim of this paper is to provide a pedagogical introduction to SFL's, fo- 
cused on the essential conceptual ideas and on issues which are settled and 
which can be expected to survive future developments. Therefore, we will not 
attempt to give an exhaustive review of the literature on this problem or of all 
the experimental systems which show hints of SFL behavior. The experimental 
examples we discuss have been selected to illustrate both what is essentially 
understood and what is not understood even in principle. On the theoretical 
side, we will shy away from presenting in depth the sophisticated methods 
necessary for a detailed evaluation of correlation functions near QCP — for 
this we refer to the book by Sachdev [210] — or for an exact solution of lo- 
cal impurity models (see, e.g., [118,212,244]). Likewise, for a discussion of the 
application of quantum critical scaling ideas to Josephson arrays or quantum 
Hall effects, we refer to the nice introduction by Sondhi et al. [232]. 



1.2 Outline of the paper 

The outline of this paper is as follows. We start by summarizing in section 2 
some of the key features of Landau's FLT — in doing so, we will not attempt 
to retrace all of the ingredients which can be found in many of the classical 
textbooks [194,36]; instead our discussion will be focused on those elements 
of the theory and the relation with its microscopic derivation that allow us 
to understand the possible routes in which the FLT can break down. This is 
followed in section 3 by the Fermi-liquid formulation of the Kondo problem 
and of the SFL variants of the Kondo problem and of two-interacting Kondo 
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impurities. The intention here is to reinforce the concepts of FLT in a different 
context as well as to provide examples of SFL behavior which offer important 
insights because they are both simple and solvable. We then discuss the prob- 
lem of one spatial dimension {d — 1), presenting the principal features of the 
solutions obtained. We discuss why d — 1 is special, and the problems en- 
countered in extending the methods and the physics to d > 1. We move then 
from the comforts of solvable models to the reality of the discussion of possible 
mechanisms for SFL behavior in higher dimensions. First we analyze in section 
5 the paradigmatic case of long range interactions. Coulomb interactions will 
not do in this regard, since they are always screened in a metal, but transverse 
electromagnetic fields do give rise to long-range interactions. The fact that as a 
result no metal is a Fermi-liquid for sufficiently low temperatures was already 
realized long ago [120] — from a practical point of view this mechanism is not 
very relevant, since the temperatures where these effects become important 
are of order 10~^^ Kelvin; nevertheless, conceptually this is important since it 
is a simple example of a gauge theory giving rise to SFL behavior. Gauge the- 
ories on lattices have been introduced to discuss problems of fermions moving 
with the constraint of only zero or single occupation per site. We then dis- 
cuss in section 6 the properties near a quantum critical point, taking first an 
example in which the ferromagnetic transition temperature goes to zero as a 
function of some externally chosen suitable parameter. We refer in this section 
to several experiments in heavy fermion compounds which are only partially 
understood or not understood even in principle. We then turn to a discussion 
of the marginal Fermi-liquid phenemenology for the SFL state of copper-oxide 
High - Tc materials and discuss the requirements on a microscopic theory that 
the phenemenology imposes. A sketch of a microscopic derivation of the phen- 
emenology is also given. We close the paper in section 8 with a discussion of 
the metallic state in = 2 and the state of the theory treating the diffusive 
singularities in d = 2 and its relation to the metal-insulator transition. 



2 Landau's Fermi-liquid 

2. 1 Essentials of Landau Fermi-liquids 

The basic idea underlying Landau's Fermi-liquid theory [144,145,194,36] is 
that of analyticity, i.e. that states with the same symmetry can be adiabati- 
cally connected. Simply put this means that whether or not we can actually 
do the calculation we know that the eigenstates of the full Hamiltonian of 
the same symmetry can be obtained perturbatively from those of a simpler 
Hamiltonian. At the same time states of different symmetry can not be ob- 
tained by "continuation" from the same state. This suggests that given a hard 
problem which is impossible to solve, we may guess a right simple problem. 
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The low energy and long wavelength excitations, as well as the correlation 
and the response functions of the impossible problem bear a one-to-one corre- 
spondence with the simpler problem in their analytic properties. This leaves 
fixing only numerical values. These are to be determined by parameters, the 
minimum number of which is fixed by the symmetries. Experiments often pro- 
vide intuition as to what the right simple problem may be: for the interacting 
electrons, in the metallic range of densities, it is the problem of kinetic en- 
ergy of particles with Fermi statistics (If one had started with the opposite 
limit, just the potential energy alone, the starting ground state is the Wigner 
crystal — a bad place to start thinking about a metal!). If we start with 
non-interacting fermions, and then turn on the interactions, the qualitative 
behavior of the system does not change as long as the system does not go 
through (or is close to) a phase transition. Because of the analyticity, we can 
even consider strongly interacting systems — the low energy excitations in 
these have strongly renormalized values of their parameters compared to the 
non-interacting problem, but their qualitative behavior is the same of that of 
the simpler problem. 

The heavy fermion problem provides an extreme example of the domain of 
validity of the Landau approach. This is illustrated in Fig. 2, which shows the 
specific heat of the heavy fermion compound CeAl^. As in the Sommerfeld 
model, the specific heat is linear in the temperature at low T, but if we write 
Cy ~ 7T at low temperatures, the value of 7 is about a thousand times as 
large as one would estimate from the density of states of a typical metal, using 
the free electron mass. For a Fermi gas the density of states A^(0) at the Fermi 
energy is proportional to an effective mass m*: 



with kp the Fermi wavenumber. Then the fact that the density of states at 
the chemical potential is a thousand times larger than for normal metals can 
be expressed by the statement that the effective mass m* of the quasiparticles 
is a thousand times larger than the free electron mass m. Likewise, as Fig. 3 
shows, the resistivity of CeAl^ at low temperatures increases as T^. This also 
is a characteristic sign of a Fermi-liquid, in which the quasiparticle lifetime r 
at the Fermi surface, determined by electron-electron interactions, behaves as 
r ~ However, just as the prefactor 7 of the specific heat is a factor 

thousand times larger than usual, the prefactor of the term in the resistivity 
is a factor 10^ larger — while 7 scales linearly with the effective mass ratio 
m*/m, the prefactor of the term in the resistivity increases for this class 

^ In heavy fermions, at least in the observed range of temperatures, the transport 
lifetime determining the temperature dependence of resistivity is proportional to 
the single-particle lifetime. 
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Fig. 2. Specific heat of CeAl^ at low temperatures from Andres et al. [28]. The 
slope of the linear specific heat is about 3000 times that of the linear specific heat 
of say Cu. But the high-temperature cut-off of this linear term is smaller than that 
of Cii by a similar amount. The rise of the specific heat in a magnetic field at low 
temperatures is the nuclear contribution, irrelevant to our discussion. 
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Fig. 3. Electrical resistivity of CeAl^ below 100 mK, plotted against T^. Prom 
Andres et al. [28]. 

of Fermi-liquids as (m*/m)^. 



It should be remarked that the right simple problem is not always easy to 
guess. The right simple problem for hquid ^ife is not the non-interacting 
Bose gas but the weakly interacting Bose gas (i.e. the Bogoliubov problem 
[44,146]). The right simple problem for the Kondo problem (a low-temperature 
local Fermi liquid) was guessed [184] only after the numerical renormalization 
group solution was obtained by Wilson [273]. The right simple problem for 
two-dimensional interacting disordered electrons in the '"metallic" range of 
densities (section 8 in this paper) is at present unknown. 

For SFL's the problem is different: usually one is in a regime of parameters 
where no simple problem is a starting point — in some cases the fluctua- 



10 




Fig. 4. Bare-particle distribution at T = for a given spin-direction in a transla- 
tionally invariant Fermi-system with interactions (full line) and without interactions 
(dashed line). Note that the position of the discontinuity, i.e. the Fermi wavenumber 
kp, IS not renormalized by interactions. 

tions between solutions to different simple problems determines the physical 
properties, while in others even this dubious anchor is lacking. 



2.2 Landau Fermi-liquid and the wave function renormalization Z 



Landau theory is the forerunner of our modern way of thinking about low- 
energy effective Hamiltonians in complicated problems and of the renormali- 
sation group. The formal statements of Landau theory in their original form 
are often somewhat cryptic and mysterious — this both reflects Landau's style 
and his ingenuity. We shall take a more pedestrian approach. 

Let us consider the essential difference between non-interacting fermions and 

an interacting Fermi-liquid from a simple microscopic perspective. For free 
fermions, the momentum states |k) are also eigenstates with energy eigenvalue 



of the Hamiltonian. Moreover, the distribution of particles is given by the 
Fermi-Dirac function for the thermal occupation n^^, where a denotes the spin 
label. At T = the distribution jumps from 1 (all states occupied within the 
Fermi sphere) to zero (no states occupied within the Fermi sphere) at|k| = kp 
and energy equal to the chemical potential /x. This is illustrated in Fig. 4. 

A good way to probe a system is to investigate the spectral function; the spec- 
tral function A{k,u!) gives the distribution of energies uj in the system when 
a particle with momentum k is added or removed from it (Remember that re- 
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Fig. 5. (a) The non-interacting spectral function A{k,uj) at fixed as a function of 
u. (b) The spectral function of single-electron excitations in a Fermi-liquid at fixed 
/c as a function oi lo. If ^A{k,uj) is normalized to 1, signifying one bare particle, the 
weight under the Lorentzian, i.e. the quasi-particle part, is Z. As explained in the 
text, at the same time Z is the discontinuity in Fig. 4. 

moving a particle excitation below the Fermi energy means that we add a hole 
excitation). As sketched in Fig. 5(a), for the non- interacting system, Ao(k, a;) 
is simply a delta- function peak at the energy ek, because all momentum states 
are also energy eigenstates: 



^o(k, uj) — S{uj — (ek — fj)) , for a; > , 



= Im — 

TT CU 



1 



1 



ImG'°(k,a;) . 



(3) 
(4) 



Here S is small and positive; it reflects that particles or holes are introduced 
adiabatically, and is taken to zero at the end of the calculation for the pure 
non-interacting problem. The first step of the second line is just a simple 
mathematical rewriting of the delta function; in the second line the Green's 
function G° for non-interacting electrons is introduced. More generally the 
single-particle Green's function G(k, a;) is defined in terms of the correlation 
function of particle creation and annihilation operators in standard textbooks 
[182,4,208,159]. For our present purposes, it is sufficient to note that it is re- 
lated to the spectral function ^(k, a;), which has a clear physical meaning and 
which can be deduced through Angle Resolved Photoemission Experiments : 



dx 



A{k, x) 



oj — iJ, — X + i6 sgn(a; — /x) 



(5) 



A(k,uj) thus is the spectral representation of the complex function G{k,uj). 
Here we have defined the so-called retarded Green's function which is especially 
useful since its real and imaginary parts obey the Kramers- Kronig relations. In 
the problem with interactions G(k, u) will differ from G°(k, a;). This difference 
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Fig. 6. Schematic illustration of the perturbative expansion (8) of the change of 
wave function as a result of the addition of an electron to the Fermi sea due to 
interactions with the particles in the Fermi sea. 

can be quite generally defined through the single-particle self-energy function 
E(k,a;): 

(G(k,a;))-^ = (G°(k,a;))"' - E(k,a;) . (6) 

Eq. (5) ensures the relation between G(k, u) and A(k, a;) 

A{k,uj)^ --lmG{k,uj) . (7) 

TT 

With these preliminaries out of the way, let us consider the form of A(k,uj) 
when we add a particle to an interacting system of fermions. 

Due to the interaction (assumed repulsive) between the added particle and 
those already in the Fermi-sea, the added particle will kick particles from below 
the Fermi-surface to above. The possible terms in a perturbative description 
of this process are constrained by the conservation laws of charge, particle 
number, momentum and spin. Those which are allowed by the aforementioned 
conservation laws are indicated pictorially in Fig. 6, and lead to an expression 
of the type 

ki,k2,k3 o"i, 0-2,0-3 

X 5k,ki-k2+k3^(t7;t7i, (72,0-3)1^^) H — ■ (8) 

Here the cj^'s and Ck's are the bare particle creation and annihilation operators, 
and the dots indicate higher order terms, for which two or more particle- 
hole pairs are created and 6{a; ai, (T2, as) expresses conservation of spin under 
vector addition. The multiple-particle-hole pairs for a fixed total momentum 
can be created with a continuum of momentums of the individual bare particles 
and holes. Therefore an added particle with fixed total momentum has a wide 
distribution of energies. However if Zk defined by Eqn. (8) is finite, there will 
be a well-defined feature in this distribution at some energy which is in general 
different from the non-interacting value fi^k"^ / {2m). The spectral function in 
such a case will then be as illustrated in Fig. 5. It is useful to separate the 
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well-defined feature from the broad continuum by writing the spectral function 
as the sum of two terms, A(k,uj) = A,,fjh{k, u) + AincohO^,^)- The single- 
particle Green's function can be expressed as a sum of two corresponding 
terms, G{k,u) = GcohiK'^) + Gj„cofe(k, a;). Then 

G^niKco)^ Z'^., , (9) 



which for large lifetimes Tk gives a Lorentzian peak in the spectral density 
at the quasiparticle energy ek = Ck — A*. The incoherent Green's function is 
smooth and hence for large Tk corresponds to the smooth background in the 
spectral density. 

The condition for the occurrence of the well-defined feature can be expressed 
as the condition on the self-energy E(k, cu) that it has an analytic expansion 
about a; = and k = k^;' and that its real part be much larger than its 
imaginary part. One can easily see that were it not so the expression (9) 
for Gcoh could not be obtained. These conditions are necessary for a Landau 
Fermi-liquid. Upon expanding E(k, a;) in (12) for small u and small deviations 
of k from k^r and writing it in the form (9), we make the identifications 

Ck = Ck^k^k , - = ^k ImE(k^^, cu^O) , (10) 

Tk 



where 




ui=0,k=kF 



From Eq. (8) we have a more physical definition of Z]^: is the projection 
amplitude of iV'k^^^) onto the state with one bare particle added to the ground 
state, since all other terms in the expansion vanish in the thermodynamic limit 
in the perturbative expression embodied by (8), 

^k/' = (v^r^i4iv^^). (12) 



In other words, Zk is the overlap of the ground state wavefunction of a system 
of interacting iV ± 1 fermions of total momentum k with the wave function of 
N interacting particles and a bare-particle of momentum k. Zk is called the 
quasiparticle amplitude. 

The Landau theory tacitly assumes that Zk is finite. Furthermore it asserts 
that for small a; and k close to ki?-, the physical properties can be calculated 
from quasiparticles which carry the same quantum numbers as the particles. 
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i.e. charge, spin and momentum and which may be defined simply by the 

\<^')=ll\r). (13) 



creation operator 7^^: 



Close to hp, and for T small compared to the Fermi-energy, the distribution 
of the quasiparticles is assumed to be the Fermi-Dirac distribution in terms of 
the renormalized quasiparticle energies. The bare particle distribution is quite 
different. As is illustrated in Fig. 4, it is depleted below kp and augmented 
above hp, with a discontinuity at T = whose value is shown in microscopic 
theory to be Zk. A central result of Fermi-liquid theory is that close to the 
Fermi energy at zero temperature, the width 1/rk of the coherent quasiparticle 
peak is proportional to (ek — /i)^ so that near the Fermi energy the lifetime 
is long and quasiparticles are well-defined. Likewise, at the Fermi energy 1/rk 
varies with temperature as T^. From the microscopic derivation of this result, 
it follows that the weight in this peak, Zk, becomes equal to the jump Z in nko- 
when we approach the Fermi surface: Z]^^ Z for h ^ hp. For heavy fermions, 
as we already mentioned, Z can be of the order of 10^^. But as long as Z is 
nonzero, one has Fermi-liquid properties for temperatures lower than about 
ZEp. Degeneracy is effectively lost for temperatures much higher than ZEp 
and classical statistical mechanics prevails. 

An additional result from microscopic theory is the so-called Luttinger the- 
orem, which states that the volume enclosed by the Fermi-surface does not 
change due to interactions [182,4]. The mathematics behind this theorem is 
that with the assumptions of FLT, the number of poles in the interacting 
Green's function below the chemical potential is the same as that for the 
non-interacting Green's function. Recall that the latter is just the number of 
particles in the system. 

Landau actually started his discussion of the Fermi-Liquid by writing the 
equation for the deviation of the (Gibbs) free-energy from its ground state 
value as a functional of the deviation of the quasiparticle distribution function 
?T,(k, cr) from the equilibrium distribution function nQ{k,a) 

6n(k, a) = n{k, a) — noCk, a) , (14) 



as follows: 



1 1 

G = Go + — XK^k - /i)5'^ka + ^TTJ hk',aa'SnkaSnk'a' ^ (15) 

k,(T kk',(TCT' 



^ It is an unfortunate common mistake to think of the properties in this regime as 
SFL behavior. 
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Note that (ek — /x) is itself a function of 5n\ so the first term contains at least 
a contribution of order [Snf' which makes the second term quite necessary. In 
principle, the unknown function fick^aa' depends on spin and momenta. How- 
ever, spin rotation invariance allows one to write the spin part in terms of two 
quantities, the symmetric and antisymmetric parts and Moreover, for 
low-energy and long- wavelength phenomena only momenta with k k^r play 
a role; if we consider the simple case of ^He where the Fermi surface is spher- 
ical, rotation invariance implies that for momenta near the Fermi momentum 
/ can only depend on the relative angle between k and k'; this allows one to 
expand in Legendre polynomials Pi{x) by writing 

^^(0)/^!.',..' ""-^"^ E FrPiii^ ■ k') . (16) 

1=0 



From the expression (15) one can then relate the lowest order so-called Landau 
coefficients Fq and and the effective mass m* to thermodynamic quantities 
like the specific heat Cy, the compressibility k., and the susceptibility %: 

^ = -=(1 + ^..,!!!:. + (17) 

CyQ m kq m Xo m 



Here subscripts refer to the quantities of the non-interacting reference sys- 
tem, and m is the mass of the fermions. For a Galilean invariant system (like 
^He), there is a a simple relation between the mass enhancement and the 
Landau parameter Ff, and there is no renormalization of the particle current 
j; however, there is a renormalization of the velocity: one has 

j = k/m , V = k/m* , !!!l=fi + ^) . (18) 

m \ 6 J 



The transport properties are calculated by defining a distribution function 
5n(ko"; r, t) which is slowly varying in space and time and writing a Boltzmann 
equation for it [194,36]. 

It is a delightful conceit of the Landau theory that the expressions of the 
low-energy properties in terms of the quasiparticles in no place involve the 
quasiparticle amplitude Z\^. In fact in a translationally invariant problem as 
liquid ^ife, Z-^ cannot be measured by any thermodynamic or transport mea- 
surements . A masterly use of conservation laws ensures that Z's cancel out 
in all physical properties (One can extract Z from measurement of the mo- 
mentum distribution. By neutron scattering measurements, it is found that 
Z pa 1/4 [108] for He' near the melting line). This is no longer true on a 
lattice, in the electron-phonon interaction problem [198] or in heavy fermions 
[250] or even more generally in any situation where the interacting problem 
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contains more than one type of particle with different characteristic frequency 
scales. 

2.3 Understanding microscopically why Fermi-liquid Theory works 

Let us try to understand from a more microscopic approach why the Landau 
theory works so well. We present a qualitative discussion in this subsection and 
outline the principal features of the formal derivation in the next subsection. 

As we already remarked, a crucial element in the approach is to choose the 
proper non-interacting reference system. That this is possible at all is due 
to the fact that the number of states to which an added particle can scatter 
due to interactions is severely limited due to the Pauli principle. As a result 
non-interacting fermions are a good stable system to perturb about: they 
have a finite compressibility and susceptibility in the ground state, and so 
collective modes and thermodynamic quantities change smoothly when the 
interactions are turned on. This is not true for non-interacting bosons which 
do not support collective modes like sound waves. So one cannot perturb about 
the non-interacting bosons as a reference system. 

Landau also laid the foundations for the formal justification of Fermi Liquid 
theory in two and three dimensions. The fiurry of activity in this field fol- 
lowing the discovery of high-Tc phenomena has led to new ways of justifying 
Fermi-liquid theory (and understanding why the one-dimensional problem is 
different). But the principal physical reason, which we now discuss, remains 
the phase space restrictions due to kinematical constraints. 

We learned in section 2.2 that to be able to define quasiparticles, it was nec- 
essary to have a finite Z^^, and that this in turn needed a self-energy function 
'E{kF,uj) which is smooth near the chemical potential, i.e. at a; = 0. Let us 
first see why a Fermi gas has such properties when interactions are calculated 
perturbatively. 

In Fig. 7 we show the three possible processes that arise in second order 
perturbation theory for the scattering of two particles with fixed initial energy 
cu and momentum q. Note that in two of the diagrams. Fig. 7(a) and 7(b) the 
intermediate state has a particle and a hole while the intermediate state in 
diagram 7(c) has a pair of particles. 

We will find that, for our present purpose, the contribution of diagram 7(a) 
is more important than the other two. It gives a contribution 
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(a) (b) (c) 






+ 1 J"'^ P2 + q P2 pT"^ P2 + q P2 ^!^k P'^ + l 

+k + q 

Fig. 7. The three second-order processes m a perturbative calculation of the correc- 
tion to the bare interaction in a Fermi-liquid 

Here g is a. measure of the strength of the scattering potential (the vertex 
in the diagram) in the hmit of small q. The denominator ensures that the 
largest contribution to the scattering comes from small scattering momenta 
q: for these the energy difference is linear in q, -Ek+q — -Ek ~ Q ■ Vk, where 
Vk is a vector of length vp 'm the direction of k. Moreover, the term in the 
numerator is nonzero only in the area contained between two circles (for d = 2) 
or spheres (for d = 3) with their centers displaced by q — here is the phase 
space restriction due to the Pauli principle! This area is also proportional to 
q ■ Vk, and so in the small q approximation we get from diagram 7(a) a term 
proportional to 

/ "-"v ., f . (20) 

cij — q ■ Vk + «o aek 



Now we see why we diagram 7(a) is special. There is a singularity at a; = q • Vk 
and its value for small u and q depends on which of the two is smaller. 
This singularity is responsible for the low energy-long wavelength collective 
modes of the Fermi liquid in Landau theory. At low temperatures, df /de^ = 
— ^(ck — f^), so the summation is restricted to the Fermi surface. The real part 
of (19) therefore vanishes in the limit qvp/uj —>■ 0, while it approaches a finite 
limit for uj —>■ 0. The imaginary part in this limit is proportional tcQ u, 

Imx{(l,uj) = g'^TTN{0) for u < qvp , (21) 



while /mx(q, tu) = for uj > vpQ- This behavior is sketched in Fig. 9i(h). 
Explicit evaluation yields for the real part 



i?ex(q,c^)=/iV(0) 



1 + In 

qvp 



UJ 



LU + qvp 



(22) 



which gives a constant (leading to a finite compressibility and spin suscep- 
tibility) at UJ small compared to qvp. For diagram 7(b), we get a term uj — 

^ This behavior implies that this scattering contribution is a marginal term in the 
renormalization group sense, which means that it affects the numerical factors, but 
not the qualitative behavior. 



18 




Fig. 8. (a) Restriction on allowed particle-hole excitations in a Fermi-sea due to 
kinematics. The plasmon mode has been drawn for the case d = 3. (b) The absorp- 
tive part of the particle-hole susceptibility (in the charge, current and spin channels) 
for UJ < qvp in the Fermi gas. 

(i?p^_P2+k+q — E-k) in the denominator. This term is always finite for general 
momenta pi and p2, and hence the contribution from this diagram can always 
be neglected relative to the one from 7(a). Along similar lines, one finds that 
diagram 7(c), which describes scattering in the particle-particle channel, is 
irrelevant except when pi = — P2, when it diverges as In a;. 

Of course, this scattering process is the one which gives superconductivity. 
Landau noticed this singularity but ignored its implication^ Indeed, as long as 
the effective interactions do not favor superconductivity or as long as we are at 
temperatures much higher than the superconducting transition temperature, 
it is not important for Fermi-liquid theory. 

Let us now look further at the absorptive spectrum of particle-hole excitations 
in two and three dimensions, i.e., we examine the imaginary part of Eq. (19). 
When the total energy uj of the pair is small, both the particle and the hole have 
to live close to the Fermi surface. In this limit, we can make any excitation with 
momentum q < 2kF- For fixed but small values of q, the maximum excitation 
energy is ~ qvp', this occurs when q is in the same direction as the main 
momentum k of each quasiparticle. For q near 2kF the maximum possible 
energy is u = vplq — 2k f\. Combining these results, we get the sketch in Fig. 
8(a), in which the shaded area in the u-q space is the region of allowed particle- 
hole excitation^. From this spectrum one can calculate the polarizability, or 
the magnetic susceptibility. 

The behavior sketched above is valid generally in two and three dimensions 
(but as we will see in section 4, not in one dimension). The important point 

^ Attractive interactions in any angular momentum channel (leading to supercon- 
ductivity) are therefore relevant operators. 

^ In the presence of long-range Coulomb interactions one gets in addition to the 
particle-hole excitation spectrum a collective mode with a finite plasma frequency 
as (; — > in d = 3 and a w ~ y/q behavior in d = 2. 
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Fig. 9. The single-particle self-energy diagram in second order. 

to remember is that the density of particle-hole excitations decreases linearly 
with uj for uj small compared to qvp- We shall see later that one way to undo 
Fermi-liquid theory is to have uj vary as k"^ in two dimensions or ~ fc^ in 
three dimensions. 

We can now use Im x(q, uj) to calculate the single-particle self-energy to second 
order in the interactions. This is shown in Fig. 9 where the wiggly line denotes 
x(q, v) which in the present approximation is just given by the diagram of Fig. 
7(a). 

For the perturbative evaluation of this process, the intermediate particle with 
energy-momentum (cj + z/), {k + q) is a free particle. Second order perturbation 
theory then yields an imaginary part or decay rate 



Im E(k, uj) = 




in three dimensions for k ^ k^. In two dimensions, the same process yields 
Im E(k^, uj) ~ uj"^ \n{Ep/uj). 

The uj"^ decay rate is intimately related to the analytic result (22) for Im x(g, uj) 
exhibited in Fig. (8). As may be found in text books, the same calculation for 
electron-phonon interactions or for interaction with spin waves in an antifer- 
romagnetic metal gives ImS(kir,ti;) ~ (uj/uJc)^, with ujc is the phonon Debye 
frequency in the former and the characteristic zone-boundary spin-wave fre- 
quency in the latter. 

The real part of the self-energy may be obtained directly or by Kramers- 
Kronig transformation of (23). It is proportional to uj. Therefore, if the quasi- 
particle amplitude Z^p is evaluated perturbativelyjj] 

Zk.^l- 2g'N{0)/EF . (24) 

Thus in a perturbative calculation of the effect of interactions the basic ana- 
lytic structure of the Green's function is left the same as for non-interacting 
fermions. The general proof of the validity of Landau theory consists in show- 
ing that what we have seen in second order in g remains true to all orders in g. 

This quantity has been precisely evaluated by Galitski [100] for the model of a 
dilute Fermi gas characterized by a scattering length. 
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The original proofs [4] are self-consistency arguments — we will consider them 
briefly in section 2.4. They assume a finite Z in exact single-particle Green's 
functions and show effectively that to any order of perturbation, the polar- 
izability functions retain the analytic structure of the non-interacting theory, 
which in turn ensures a finite Z . 

In one dimension, phase space restrictions on the possible excitations are cru- 
cially different^, since the Fermi surface consists of just two points in the 
one-dimensional space of momenta — see Fig. 10. As a result, whereas in 
d = 2 and = 3 a continuum of low-energy excitations with finite q is possi- 
ble, at low energy only excitations with small k or k ^ 2kF are possible. The 
subsequent equivalence of Fig. 8 for the one- dimensional case is the one shown 
in Fig. 11. Upon integrating over the momentum k with a cut-off of 0{kp) 
the contribution from this particle-hole scattering channel to Rex(g, tu) is 

/ dk -4 ^ ln[(cj + qvF)/EF] . (25) 



(Note that (25) is true for both q << kp and \q — 2kF\ « kp.) This in turn 
leads to a single particle self-energy calculated by the process in Fig. 9 to 
be Re E(ki?, tu) ~ t^lnti; giving a hint of trouble. Also the Cooper (particle- 
particle) channel has the same phase space restrictions, and gives a contribu- 
tion also proportional to Intu. The fact that these sngular contributions are of 
the same order, leads to a competition between charge/spin fluctuations and 
Cooper pairing fluctuations, and to power law singularities. Also, the fact that 
instead of a continuum of low energy excitations as in higher dimensions, the 
width of the band of allowed particle hole excitations vanishes as — 0, is 
the reason that the properties of one-dimensional interacting metals can be 
understood in terms of bosonic modes. We will present a brief summary of 
the results for the single-particle Green's function and correlation functions in 
section 4.9. 

In special cases of nesting in two or three dimensions, one can have situations 
that resemble the above one- dimensional case. When the non-interacting Fermi 
surface in a tight binding model has the square shape sketched in Fig. 12(a) 
— this happens for a tight-binding model with nearest neighbor hopping on 
a square lattice at half-filling — a continuous range of momenta on opposite 
sides of the Fermi surface can be transformed into each other by one and the 
same wavenumber. This so-called nesting leads to log and log^ singularities 
for a continuous range of k in the perturbation theory for the self energy 
S(k, cu). Likewise, the partially nested Fermi surface of Fig. 12(h) leads to 
charge density wave and antiferromagnetic instabilities. We will come back to 

^ It might appear surprising that they are not different in any essential way between 
higher dimensions. 
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Fig. 10. Single-particle energy in one dimension, in the approximation that 
the dispersion relation is linearized about kp. Note the Fermi surface con- 
sists of just two points. The spectrum of particle-hole excitations is given by 
uj{q) = e{k + q) — e{k) = kpg/m. Low-energy particle-hole excitations are only 
possible for q small or for q near 2kF- 



Fig. 11. Phase space for particle-hole excitation spectrum in one dimension. Com- 
pare with the same in higher dimensions, Fig. 8. For linearized single-particle kinetic 
energy = ±VF{k — kp), particle-hole excitations are only possible on lines going 
through k = and k = 2kF- 

these in sections 2.6 and 6. 

2.4 Principles of the Microscopic Derivation of Landau Theory 

In this section we will sketch how the conclusions in the previous section 
based on second-order perturbation calculation are generalized to all orders in 
perturbation theory. This section is slightly more technical than the rest; the 
reader may choose to skip to section 2.6. 

We follow the microscopic approach whose foundations were laid by Landau 
himself and which is discussed in detail in excellent textbooks [184,194,36,4]. 
For more recent methods with the same conclusions see [222,121]. Our empha- 
sis will be on highlighting the assumptions in the theory so that in the next 
section we can summarize the routes by which the Fermi-liquid theory may 
break down. These assumptions are usually not stated explicitly. 




k 
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Fig. 12. (a) The nested Fermi surface on obtained in a tight binding model on a 
square lattice with nearest neighbor hopping; (b) A partially nested Fermi surface 
which leads to charge density wave or antiferromagnetic instabilities 




P5 P2-k 

Fig. 13. Diagrammatic representation of Eq. (27). 

The basic idea is that due to kinematic constraints, any perturbative process 
with n particle-hole pairs in the intermediate state gives contributions to the 
polarizability proportional to (ou/Ef)"'. Therefore the low energy properties 
can be calculated with processes with the same "skeletal" structure as those 
in Fig. 7, which have only one particle-hole pair in the intermediate state. So 
one may concentrate on the modification of the four-legged vertices and the 
single-particle propagators due to interactions to all orders. Accordingly the 
theory is formulated in terms of the single particle Green's function G{p) and 
the two-body scattering vertex, 

r(pi, P2, Pi + k, p2 -k) =r(pi, p2, k) (26) 

Here and below we use, for sake of brevity, p, etc. to denote the energy- 
momentum four- vector (p, uj) and we suppress the spin labels. The equation 
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for r is expanded in one of the two particle-hole channels as[^ 



-if T'^^){p,,q,k) G{q)G{q + k)T{q,p2,k)-0^^ , (27) 

where T^^^ is the irreducible part in the particle hole channel in which Eq. 
(27) is expressed. In other words, T^^^ can not be split up into two parts by 
cutting two Green's function lines with total momentum k. So F*^^) includes 
the complete vertex in the other (often called cross-) particle-hole channel. 
The diagrammatic representation of Eq. (27) is shown in Fig. 13. In the sim- 
plest approximation F^^-* is just the bare two-body interaction. Landau theory 
assumes that F*^^) has no singularities!^ An assumption is now further made 
that G{p) does have a coherent quasiparticle part at |p| — Pf and ~ 0, 

G{p) = —^^ + G,nc , (28) 

u-ep + id sgn(ep) 



where ip is to be identified as the excitation energy of the quasiparticle, Z its 
weight, and Ginc the non-singular part of G (The latter provides the smooth 
background part of the spectral function in Fig. 5(b) and the former the sharp 
peak, which is proportional to the 5 function for Cp = Cp — fi. It follows [182,4] 
from (28) that 

G{q)G{q + k) = 3 -6{u)6{\q\ - p^) + 0(g) , (29) 

Vf — Vr, ■ k 



for small k and cu, and where v and {v + uj) are frequencies of the two Green's 
functions. Note the crucial role of kinematics in the form of the first term which 
comes from the product of the quasiparticle parts of G; 0(g) comes from the 
scattering of the incoherent part with itself and with the coherent part and is 
assumed smooth and featureless (as it is indeed, given that G^nc is smooth and 
featureless and the scattering does not produce an infrared singularity at least 
perturbatively in the interaction). The vertex F in regions close to k k^ and 
a; ~ is therefore dominated by the first term. The derivation of Fermi-liquid 
theory consists in proving that the equations (27) for the vertex and (28) for 
the Green's function are mutually consistent. 

^ To second order in the interactions the correction to the vertex in the two possible 
particle-hole channels has been exhibited in the first two parts of Fig. 7. 
^'^ The theory has been generalized for Coulomb interactions [194, 184, 4]. The general 
results remain unchanged because a screened short-range interaction takes the place 
of r^^^ .This is unlikely to be true in the critical region of a metal-insulator transition, 
because on the insulating side the Coulomb interaction is unscreened. 
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The proof proceeds by defining a quantity r^(pi,p2, k) through 



r'^(pi,P2,^) = r«(pi,P2,fc) 



(30) 



r'^ contains repeated scattering of the incoherent part of the particle- hole pairs 
among itself and with the coherent part, but no scattering of the coherent part 
with itself. Then, provided the irreducible part of T^^'' is smooth and not too 
large, is smooth in k because 4>{q) is by construction quite smooth. 

Using the fact that the first part of (29) vanishes for VF\k\/u; — > 0, and com- 
paring (27) and (30) one can write the forward scattering amplitude 



lim 



lim r{pi,p2,k) 



= r-(pi,p2) . 



(31) 



This is now used to write the equation for the complete vertex F in terms of 



r{pi,p2,k)^r^{pi,P2) 

+ 7^r^ / r-(pi,g)r(g,p2,fc) '''''^j n,, (32) 

(27r)"^'t;F J a; — Vq • k 

where in the above \q\ = Pf and one integrates only over the solid angle Qg. 

Given a non-singular V^, a non-singular T is produced (unless the denominator 
in Eq. (32) produces singularities after the indicated integration-the Landau- 
Pomeranchuk singularities discussed below). The one-particle Green's function 
G can be expressed exactly in terms of F — see Fig. 14. This leads to Eq. (28) 
proving the self-consistency of the ansatz with a finite quasiparticle weight Z. 
The quantity Z'^T^ is then a smooth function and goes into the determination 
of the Landau parameters. 

The Landau parameters can be written in terms of the forward scattering am- 
plitude. In effect they parametrize the momentum and frequency independent 
scattering of the incoherent parts among themselves and with the coherent 
parts so that the end result of the theory is that the physical properties can 
be expressed purely in terms of the quasiparticle part of the single-particle 
Green's function and the Landau parameters. No reference to the incoherent 
parts needs to be made for low energy properties. For single-component trans- 
lational invariant fermions (like liquid '^He) even the quasiparticle amplitude 
Z disappears from all physical properties. This last is not true for renormal- 
ization due to electron-phonon interactions and in multi-component systems 
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Fig. 14. Diagram for the exact single-particle self-energy in terms of tlie exact vertex 
r and the exact single-particle Green's function. 



such as heavy fermions. Special simplifications of the Landau theory occur in 
such problems and in other problems where the single-particle self-energy is 
nearly momentum independent [168,264,101,250,175]. 

As we also mentioned, the single-particle self-energy S can be written exactly 
in terms of the vertex F: the relation between the two is represented diagram- 
matically in Fig. 14. The relations between E and F are due to conservation 
laws which Landau theory, of course, obeys. But the conservation laws are 
more general than Landau theory. It is often more convenient to express these 
conservation laws as relations between the self-energy and the three-point ver- 
tices, Aa(pe; qu) which couple external perturbations to either the density (the 
fourth component, a — 4) or the current density in the a — (1, 2, 3) direction 
The diagrammatic representation of the equation for A is shown in Fig. 15. 
The following relations (Ward identities) have been proven for translationally 
invariant problems: 



lim { - ^ 0, ? ^ o) A„(pe; qu)^^- ^E(p, e) , (a = 1, 2, 3) (33) 
\q J m opa 

]ira{-^0, g ^ 0^ A4(pe; qu;) = 1 + , (34) 

0, e ^ O) A,(pe; qo;) = ^ - /-E(p, e) , (a = 1, 2, 3) (35) 
/ m apa 

lim ^ 0, e ^ o) A4(pe; qcu) = 1 + ^^^^ ■ (36) 




A relation analogous to (34) is derived for fields couphng to spin for the case 
that interactions conserve spin. The total derivative in (35) and (36) [rather 
than the partial derivative in (33) and (34)] represents that dH/dfi is the 
variation in S when e is changed to e + de together with /x to /i + d/i, and 
dTi/dpa represents the variation when the momentum p as well as the Fermi- 
surface is translated by dp^- 
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A = A + XG G I 

Fig. 15. Vertex for coupling to external perturbations at energy-momentum {uj,q); 
A is the bare vertex. 

Eq. (36) is an expression of energy conservation, and Eq. (34) of particle num- 
ber conservation. Eqs. (33) and (34) together signify the continuity equation. 
Eq. (35) represents current conservationf"^. 

In Landau theory, the right hand sides in Eqs. (33)-(36) are expressible in 
terms of the Landau parameters. These relations are necessary to derive the 
renormalization in the various thermodynamic quantities that we quoted in 
Eqs. (17) and (18) as well as the Landau transport equation. Needless to say, 
any theory of SFL must also be consistent with the Ward identities. 



2.5 Modem derivations 

The modern derivations of Fermi-liquid theory start as well by assuming the 
existence of a Fermi surface. Kinematics then inevitably leads to similar con- 
siderations as above. Instead of the division into coherent and incoherent part 
made in Eq. (28), the renormalization group procedures are used to systemat- 
ically generate a successively lower energy and small momentum Hamiltonian 
with excitations of particle ever closer to the Fermi surface. The calculations 
are done either in terms of fermions [222] or newly developed bosonization 
methods in arbitrary dimensions [121]. The end result is equivalent to Eqs. 
(28), (30) and (32). These methods may well turn out to be very important 
in finding the structure of SFL's and in systematizing them. 

These derivations do the calculation in arbitrary dimension d and conclude 



The Ward identity Eq. (35) does not hold for impure system where the Fermi- 
surface cannot be defined in momentum space. Since energy is conserved, a Fermi 
surface can still be defined in energy space, and hence the other Ward identities 
continue to hold. This point is further discussed in section 7. 
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that the forward scattering amphtude is 



r"(pi,P2,A;) 



( 



k 



) 



fiPi,P2,k) , 



(37) 



where / is a smooth function of all of its arguments. In one-dimension the 
forward scattering amplitude has a logarithmic singularity, as we noted earlier. 

We can rephrase the conceptual framework of Landau Fermi-liquid theory in 
the modern language of Renormalization Group theory [222]. As we discussed, 
in Fermi-Liquid theory one treats a complicated strongly interacting fermion 
problem by writing the Hamiltonian Ti as 



In our discussion Tisimpie was the non-interacting Hamiltonian. The non-inter- 
acting Hamiltonian is actually a member of a 'line' of fixed-point Hamiltonians 
a.* all of which have the same symmetries but differ in their Landau parame- 
ters Fi'"" etc. The F^'s, obtained from the forward scattering in Landau theory 
are associated with marginal operators and distinguish the properties of the 
various systems associated with the line of fixed points. Landau fermi-liquid 
theory is first of all the statement of the domain of attraction of this line of 
fixed points. The theory also establishes the universal low temperature prop- 
erties due to the "irrelevant" operators generated by Tirest due to scattering 
in channels other than the forward channel. Landau theory does not estab- 
lish (at least completely) the domain of attraction of the "critical surface" 
bounding the domain of attraction of the Fermi-liquid fixed line from those of 
other fixed points or lines. If Tirest were to generate a "relevant" operator — 
i.e., effective interactions which diverge at low energies (temperatures) — the 
scheme breaks down. For example, attractive interactions between fermions 
generate relevant operators-they presage a transition to superconductivity, a 
state of different symmetry. But if we stay sufficiently above Tc, we can usually 
continue using Landau theoryp^ 



We note that in a Renormalization Group terminology, all Landau parameters 
/kk',(T(T' originating from forward scattering (i.e. zero momentum transfer), are 
"marginal operators" [143,222]. All other operators that determine finite tempera- 
ture observable properties are "irrelevant". Thus, in a "universal" sense, condensed 
matter physics may be deemed to be an "irrelevant" field. So much for technical 
terminology! 
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2. 6 Routes to Breakdown of Landau Theory 



From Landau's phenomenological theory, one can only say that the theory 
breaks down when the physical properties — specific heat divided by temper- 
ature, [^compressibility, or the magnetic susceptibility — diverge or when the 
collective modes representing oscillations of the Fermi-surface in any harmonic 
and singlet or triplet spin combinations become unstable. The latter, called 
the Landau-Pomeranchuk singularities, are indeed one route to the breakdown 
of Landau theory and occur when the Landau parameters Fi''^ reach the crit- 
ical value — (2/ + 1). A phase transition to a state of lower symmetry in then 
indicated. The new phase can again be described in Landau theory by defining 
distribution functions consistent with the symmetry of the new ground state. 

The discussion following Eq. (8) in section 2.2 allows us to make a more general 
statement. Landau theory breaks down when the quasiparticle amplitude 
becomes zero; i.e. when the state c\\iIj^) and iV'k^'*'^) ^^e orthogonal. This can 
happen if the series expansion in Eq. (8) in terms of the number of particle- 
hole pairs is divergent. In other words, addition of a particle or a hole to the 
system creates a divergent number of particle-hole pairs in the system so that 
the leading term does not have a finite weight in the thermodynamic limit. 
From Eq. (11) connecting the Z's to S's, this requires that the single-particle 
self-energy be singular as a function of a; at A; ~ /ci?. This in turn means that 
the Green's functions of SFL's contain branch cuts rather than the poles unlike 
Landau Fermi-liquids. The weakest singularity of this kind is encountered in 
the borderline "marginal Fermi-liquids" wherep^ 



If a divergent number of low-energy particle-hole pairs is created upon addition 
of a bare particle, it means that the low-energy response functions (which 
all involve creating particle-hole pairs) of SFL's are also divergent. Actually 
the single-particle self-energy can be written in terms of integrals over the 
complete particle-hole interaction vertex as in Fig. 14. The implication is that 



' The specific heat of a system of fermions can be written in terms of integrals 
over the phase angle of the exact single-particle Green's function ([4]). Given any 
singularity in the self-energy, C^/T is never more singular than InT. This accounts 
for the numerous experimental examples of such behavior we will come across. 

To see why this is the borderline case, note that a requisite for the definition 
of a quasiparticle is that the quasiparticle peak width = 2S should vanish 
faster than linear in w, the quasiparticle energy. Thus S ~ u; is the first power for 
which this is not true. The lv\'b.{uJc./<^) term in Eqn. (39) is then dictated by the 
Kramers-Kroning relation. 




(39) 
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the interaction vertices are actually more divergent than the single-particle 
self-energy. 

Yet another route to SFL's is the case in which the interactions generate 
new quantum numbers which are not descriptive of the non-interacting prob- 
lem. This happens most famously in the Quantum Hall problems and in one- 
dimensional problems (section 4) as well as problems of impurity scattering 
with special symmetries (section 3). In such cases the new quantum numbers 
characterize new low-energy topological excitations. New quantum numbers 
of course imply Z = 0, but does Z = imply new quantum numbers. One 
might wish to conjecture that this is so. But there is no proof of thi^^. 

In the final analysis all breakdowns of Landau theory are due to degeneracies 
leading to singular low-energy fluctuations. If the characteristic energy of the 
fluctuations is lower than the temperature, a quasi-classical statistical me- 
chanical problem results. On the basis of our qualitative discussion in section 
2.3 and the sketch of the microscopic derivation in section 2.4, we may divide 
up the various routes to breakdown of Landau theory into the following (not 
necessarily orthogonal) classes: 

(i) Landau-Pomeranchuk Singularities: Landau theory points to the possibility 
of its breakdown through the instability of the collective modes of the Fermi- 
surface which arise from the solution of the homogeneous part of Eqn.(32). 
These collective modes can be characterized by the angular momentum d. of 
oscillation of the Fermi-surface and whether the oscillation is symmetric "s" 
or anti-symmetric "a" in spin. The condition for the instability derived from 
the condition of zero frequency of the collective modes are [194,36] 

F/ < -{21 + 1) , Fl < -(2£ + 1) . (40) 



The £ = conditions refer to the divergence in the compressibility and the 
(uniform) spin-susceptibility. The former would in general occur via a first- 
order transition, so is uninteresting to us. The latter describes the ferromag- 
netic instability. No other Landau-Pomeranchuk instabilities have been exper- 
imentally identified. But such new and exotic possibilities should be kept in 
mind. Thus, for example, an F{*-instability corresponds to the Fermi- velocity 
— > 0, a instability to a "d-wawe-like" instability of the particle- hole excita- 
tions on the Fermi-surface etc. Presumably these instabilities are resolved by 
reconstruction of the Fermi-surface with (patches) of energy gaps. Coupling of 
the damped transverse-excitations of charged-fermions to zero-point fluctua- 
tions of the electromagnetic-fields produces an SFL which we study in section 
5.1. The microscopic interactions necessary for the Landau-Pomeranchuk in- 

It would indeed be a significant step forward if such a conjecture could be proven 
to be true or if the conditions in which it is true were known. 
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stabilities and the critical properties near such instabilities have not been well 
investigated, especially for fermions with a lattice potential. 

It is also worth noting that some of the instabihties are disallowed in the 
limit of translational invariance. Thus, for example, time-reversal breaking 
states, such as the "anyon-state" [148,59] cannot be realized because in a 
translationally invariant problem the current operator cannot be renormalized 
by the interactions, as we have learnt from Eqs. (18), (33). 

(ii) Critical regions of Large Q-Singularities: Landau theory concerns itself 
only with long wavelength response and correlations. A Fermi-liquid may 
have instabilities at a nonzero wave- vector, for example a charge-density wave 
(CDW) or spin-density wave (SDW) instability. Only a microscopic calculation 
can provide the conditions for such instabilities and therefore such conditions 
can only be approximately derived. An important point to note is that they 
arise perturbatively from repeated scattering between the quasiparticle parts 
of G while the scattering vertices (irreducible interactions) are regular. The su- 
perconductive instability for any angular momentum is also an instability of 
this kind. In general such instabilities are easily seen in RPA and/or matrix 
calculations. 

Singular Fermi-liquid behavior is generally expected to occur only in the crit- 
ical regime of such instabilities [112,156]. If the transition temperature Tc is 
finite then there is usually a stable low temperature phase in which unstable 
modes are condensed to an order parameter, translational symmetry is bro- 
ken, and gaps arise in part or all of the Fermi-surface. For excitations on the 
surviving part of the Fermi-surface, Fermi-liquid theory is usually again valid. 
The fluctuations in the critical regime are classical, i.e. with characteristics 
frequency a;/^ <^ ksTc- 

If the transition is tuned by some external parameter so that it occurs at 
zero temperature, one obtains, as illustrated already in Fig. 1, a Quantum 
Critical Point (QCP). If the transition is approached at T = as a function 
of the external parameter, the fluctuations are quantum-mechanical, while if 
it is approached as a function of temperature for the external parameter at its 
critical value, the fluctuations have a characteristic energy proportional to the 
temperature. A large region of the phase diagram near QCP's often carries 
SFL properties. We shall discuss such phenomena in detail in section 6. 

(Hi) Special Symmetries: The Cooper instability at g = 0, Fig. 1(c), is due to 
the "nesting" of the Fermi-surface in the particle-particle channel. Usually in- 
dications of flnite g-CDW or SDW singularities are evident pertubatively from 
Fig. T (a) or Fig. 1(h) for special Fermi-surfaces, nested in some q-direction 
in particle-hole channels. One-dimensional fermions are perfectly nested in 
both particle- hole channels and particle-particle channels (Figs. 7 (a) -(c)) and 
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hence they are both logarithmically singular. Pure one- dimensional fermions 
also have the extra conservation law that right going and left going momenta 
are separately conserved. These introduce special features to the SFL of one- 
dimensional fermions such as the introduction of extra quantum numbers. 
These issues are discussed in section 4.9. Several soluble impurity problems 
with special symmetries have SFL properties. Their study can be illuminating 
and we discuss them in section 3. 

(iv) Long- Range Interactions: Breakdown of Landau Fermi-liquid may come 
about through long-range interactions, either in the bare Hamiltonian through 
the irreducible interaction or through a generated effective interaction. The 
latter, of course, happens in the critical regime of phase transitions such as 
discussed above. Coulomb interactions will not do for the former because of 
screening of charge fluctuations. The fancy way of saying this is that the 
longitudinal electromagnetic mode acquires mass in a metal. The latter is not 
true for current fluctuations or transverse electromagnetic modes which due 
to gauge invariance must remain massless. This is discussed in section 5.1, 
where it is shown that no metal at low enough temperature is a Fermi-liquid. 
However, the cross-over temperature is too low to be of experimental interest. 

An off-shoot of an SFL through current fluctuations is the search for ex- 
tra (induced) conservation laws for some quantities to keep their fluctuations 
massless. This line of investigation may be referred to generically as gauge the- 
ories. Extra conservation laws imply extra quantum numbers and associated 
orthogonality. We discuss these in section 5.2. The one-dimensional interact- 
ing electron problem and the Quantum Hall effect problems may be usefully 
thought of in these terms. 

(v) Singularities in the Irreducible Interactions: In all the possibilities dis- 
cussed in (i)-(iii) above the irreducible interactions F^^^ defined after Eq. 
(27) are regular and not too large. As noted after Eq. (30) this is necessary to 
get a regular F*^. When these conditions are satisfied the conceivable singu- 
larities arise only from the repeated scattering of low-energy particle-hole (or 
particle-particle) pairs) as in Eq. (32) or its equivalent for large momentum 
transfers. 

A singularity in the irreducible interaction of course invalidates the basis of 
Landau theory. Such singularities imply that the parts of the problem con- 
sidered harmless perturbatively because they involve the incoherent and high 
energy parts of the single-particle spectral weight as in Eq. (30) are, in fact, 
not so. This is also true if F^-*^^ is large enough that the solution of Eq. (30) is 
singular. Very few investigations of such processes exist. 

How can an irreducible interaction be singular when the bare interaction is 
perfectly regular? We know of two examples: 



32 



E 




1- 



Fig. 16. A model band structure for a solid with more than one atom per unit 
cell. Actually in CuBq where excitonic singularities have been invoked to produce 
a ferromagnetic state [280] there are three equivalent points in the Brillouin zone 
where the conduction band minima and the valence band maxima occur. 

In disordered metals the density correlations are diffusive with characteristic 
frequency uj scaling with q^. The irreducible interactions made from the dif- 
fusive fluctuations and interactions are singular in c? = 2. This gives rise to a 
new class of SFL's which are discussed in section 8. One finds that in this case 
the singularity in the cross-particle-hole channel (the channel different from 
the one through which the irreducibihty of F^^^ is defined) feeds back into a 
singularity in F*^^''. This is very special because the cross-channel is integrated 
over and the singularity in it must be very strong for this to be possible. 

The second case concerns the particle singularities in the irreducible interac- 
tions because of excitonic singularities. Usually the excitonic singularities due 
to particle-hole between different bands occur at a finite energy and do not 
introduce low energy singularities. But if the interactions are strong enough 
these singularities occur near zero frequency. In effect eliminating high energy 
degrees of freedom generates low energy irreducible singular vertices. 

Consider, for example, the band-structure of a solid with more than one atom 
per unit cell with (degenerate) valence band maxima and minima at the same 
points in the Brillouin zone, as in Fig. (16) . 

Let the conduction band be partially filled, and the energy difference between 
fi and the valence band marked Eq in Fig. 16 be much smaller than the attrac- 
tive particle- hole interactions V between states in the valence (v) band and 
the conduction band (c). For any finite V excitonic resonances form from scat- 
tering between v and c states, as in the X-ray edge problem to be discussed 
in section 3.5. For large enough V such resonances occur at asymptotically 
low energy so that the Fermi-liquid description of states near the chemical 
potential in terms of irreducible interaction among the c-states is invalid. The 
effective irreducible interactions among the low-energy states integrate over 
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the excitonic resonance and will in general be singular if the resonance is near 
zero-energy. 

Such singularities require interactions above a critical magnitude and are 
physically and mathematically of an unfamiliar nature. In a 2-band one- 
dimensional model, exact numerical calculations have established the impor- 
tance of such singularities [237,234]. 

Recently it has been found that CaB^ or SrE^ with low density of trivalent 
Eu or quadrivalent Ce ions substituting for (Co, Sr) is a ferromagnet [280]. 
The most plausible explanation [288,33,35] is that this is a realization of the 
excitonic ferromagnetism predicted by Volkov et al. [263]. The instability to 
such a state occurs because the energy to create a hole in the valence band and 
a particle in the conduction band above the fermi-energy goes to zero if the 
attractive particle-hole (interband) interactions are large enough. This problem 
has been investigated only in the mean-field approximation. Fluctuations in 
the critical regime of such a transition are well worth studying. 

Excitonically induced singularities in the irreducible interactions are also re- 
sponsible for the Marginal Fermi-liquid state of Cu— O metals in a theory to 
be discussed in section 7. 



3 Local Fermi-Liquids &; Local Singular Fermi-Liquids 

In this section we discuss a particular simple form of Fermi-liquid formed by 
electrons interacting with a dilute concentration of magnetic impurity. Many 
of the concepts of Fermi-liquid theory are revisited in this problem. Variants of 
the problem provide an interesting array of soluble problems of SFL behavior 
and illustrate some of the principal themes of this article. 



3.1 The Kondo Problem 

The Kondo problem is at the same time one of the simplest and one of the 
most subtle examples of the effects of strong correlation effects in electronic 
systems. The experiments concern metals with a dilute concentration of mag- 
netic impurities. In the Kondo model one considers only a single impurity; the 
Hamiltonian then is 

n = tY^ cl^Cja + JS • cjaco , (41) 

{ij) 
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where {cia, c]^) denote the annihilation and creation operators of a conduction 
electron at site i with projection a in the ^-direction of spin a. The second 
term is the exchange interaction between a single magnetic impurity at the 
origin (with spin S = 1/2) and a conduction electron spin. 

When the exchange constant J > the system is a Fermi-liquid. Although 
not often discussed, the ferromagnetic (J < 0) variant of this problem is one 
of the simplest examples of a singular Fermi-liquid. 

There are two seemingly simple starting points for the problem: (i) J = 0: 
This turns out to describe the unstable high temperature fixed pointf^. The 
term proportional to J is a marginal operator about the high temperature 
fixed point because as discovered by Kondo[138] in a third order perturbation 
calculation, the effective interaction acquires a singularity ~ J^/t"^ ln{t/uj). (ii) 
t = 0: The perturbative expansion about this point is well behaved. This turns 
out to describe the low temperature Fermi-liquid fixed point. One might be 
surprised by this, considering that typically the bare t/ J is of order 10"^^. But 
such is the power of singular renormalizationsP^. 

For the reader unfamiliar with reading a renormalization group diagram like 
that of Fig. 17(b) or 18, the following explanation might be helpful. The flow in a 
renormalization group diagram signifies the following. The original problem, with 
bare parameters, corresponds to the starting point in the parameter space in which 
we plot the flow. Then we imagine "integrating out" the high energy scales (e.g. 
virtual excitations to high energy states); effectively, this means that we consider 
the system at lower energy (and temperature) scales by generating effective Hamil- 
tonians with new parameters so that the low energy properties remain invariant. 
The "length" along the flow direction is essentially a measure of how many energy 
scales have been integrated out — typically, as in the Kondo problem, this decrease 
is logarithmic along the trajectory. Thus, the regions towards which the flow points 
signify the effective parameters of the model at lower and lower temperatures. Fixed 
points towards which all trajectories flow in a neighborhood describe the universal 
low temperature asymptotic behavior of the class of models to which the model 
under consideration belongs. When a fixed point of the flow is unstable, it means 
that a model whose bare parameters initially lie close to it flows away from this 
point towards a stable fixed point; hence it has a low-temperature behavior which 
does not correspond to the model described by the unstable fixed point. A fixed line 
usually corresponds with a class of models which have some asymptotic behavior, 
e.g. an exponent, which varies continuously. 

A particularly lucid discussion of the renormalization procedure may be found 
in [143]. Briefly, the procedure consists in generating a sequence of Hamiltonians 
with successively lower energy cut-offs that reproduce the low-energy spectrum. All 
terms allowed by symmetry besides those in the bare Hamiltonian are kept. The 
coeffcients of these terms scale with the cut-off. Those that decrease proportionately 
to the cut-off or change only logarithmically, are coefficients of marginal operators, 
those that grow/decrease (algebraically) of relevant /irrelevant operators. Marginal 
operators are marginally relevant or marginally irrelevant. Upon renormalization the 
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The interaction between conduction electrons and the locahzed electronic level 
is not a direct spin interaction. It originates from quantum-mechanical charge 
fluctuations that (through the Pauli principle) depend on the relative spin ori- 
entation. To see this explicitly it is more instructive to consider the Anderson 
model [19] in which 

(ij) " 

+ E(^^i.Co,. + /i.c.) (42) 

fe,(T 

The last term in this Hamiltonian is the hybridization between the localized 
impurity state and the conduction electrons, in which spin is conserved. In the 
particle- hole symmetric case, = —U /2 is the one-hole state on the impurity 
site in the Hartree-Fock approximation and the one-particle state has the 
energy U/2. 

Following a perturbative treatment in the limit p^, ^ ^ 1 the Anderson model 
reduces to the Kondo Hamiltonian with an effective exchange constant Jg// ~ 

{vytf/u. 

The Anderson model has two simple hmits, which are illustrated in Fig. 17: 

(i) V = 0: This describes a local moment with Curie susceptibility % ~ /Xg/T. 
This limit is the correct point of departure for an investigation for the high 
temperature regime. As noted one soon encounters the Kondo divergences. 

(ii) U = 0: In this limit the impurity forms a resonance of width F ~ V^/t 
at the chemical potential which in the particle-hole symmetric case is half- 
occupied. The ground state is a spin singlet. This limit is the correct starting 
point for an examination of the low temperature properties (T <^ Tk)- A 
temperature independent contribution to the susceptibility and a linear con- 
tribution to the specific heat (~ A^(0)T/F) are contributed by the resonant 
state. Hence the name local Fermi-Liquid. 

The conceptually hard part of the problem was to realize that (ii) is the cor- 
rect stable low temperature fixed point and the technically hard problem is to 
derive the passage from the high-temperature regime to the low-temperature 
regime. This was first done correctly by Wilson [273] through the invention of 
the Numerical Renormalization Group (and almost correctly by Anderson and 



flow is to the strong coupling J = oo flxed point, sec Fig. 17. The terms generated 
from t ^ serve as irrelevant operators at this fixed point; this means that they do 
not affect the ground state but determine the measurable low-energy properties. 
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Fig. 17. (^aj Hartree-Fock excitation spectrum of the Anderson model in the two 
limits of zero hybridization, V = and zero interaction, U = 0. (b) Renormalization 
group flow of the Kondo problem. 

Yuval [22,23] by analytic methods). The analysis showed that under Renormal- 
ization Group scaling transformations the ratio {J/t) increases monotonically 
as illustrated in Fig. 17(b) — continuous RG flows are observed from the high 
temperature extreme (i) to the low temperature extreme (ii) and a smooth 
crossover between between the two regimes occurs at the Kondo temperature 



Because all flow is towards the strong-coupling fixed point, universal forms 
for the thermodynamic functions are found. For example, the specific heat Cy 
and the susceptibility x scale as 



where the fs are universal scaling functions. 

An important theoretical result is that compared to a non-interacting reso- 
nant level at the chemical potential, the ratio of the magnetic susceptibility 
enhancement to the specific heat enhancement. 




(43) 





(44) 
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for spin 1/2 impurities at T -C Tk is precisely 2 [273,184]. In a noninteracting 
model, this ratio, nowadays called the Wilson ratio, is equal to 1, since both x 
and C„ are proportional to the density of states A^(0). Thus the Wilson ratio 
is a measure of the importance of correlation effects. It is in fact the analogue 
of the Landau parameter Fq of Eq. (18). 

3.2 Fermi-liquid Phenernenology for the Kondo problem 

Following Wilson's solution [273], Nozieres [184] showed that the low-tempe- 
raturc properties of the Kondo problem can be understood simply through 
a (local) Fermi-liquid framework. This is a beautiful example of the applica- 
tion of the concept of analyticity and of symmetry principles about a fixed 
point. We present the key arguments below. For the application of this line of 
approach to the calculation of a variety of properties we refer the reader to 
papers by Nozieres and Blandin [184,185]. 

The properties of a local impurity can be characterized by the energy-dependent 
s-wave phase shift 5cr(e), which in general also depends on the spin of the con- 
duction electron being scattered. In the spirit of Fermi-liquid theory the phase 
shift may be written in terms of the deviation of the distribution function 5n(e) 
of conduction electrons from the equilibrium distribution, 

5a{.e) = 5o(e) + Y.^aa\e,e')5n,\e') + ... . (46) 

e'a' 

About a stable fixed point the energy dependence is analytic near the chemical 
potential (e = 0), so that we may expand 

6o{t) = 6o + ae + ... , Saa'{e,e') = 4>aa' + ■■■ ■ (47) 

Just as the Landau parameters are expressed in terms of symmetric and an- 
tisymmetric parts, we can write 

0TT = <Pii^<P' + r , , <f>u = <t>U ■ (48) 

Taken together, this leaves three parameters a, (j)'^ and (j)"" to determine the 
low-energy properties. Nozieres [184] showed that in fact there is only one 
independent parameter (say a which is of 0{1/Tk), with a prefactor which 
can be obtained by comparing with Wilson's detailed numerical solution). To 
show this note that by the Pauli principle same spin states do not interact, 
therefore [184] 

= 0- + 0« = . (49) 
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Secondly a shift of the chemical potential hy S/j, and a simultaneous increase in 
Sn by N{(})Sfi should have no effect on the phase shift, since the Kondo-effect 
is tied to the chemical potential. Therefore according to 46) and (47) 

[a + N(0)(l)']5ii = , =^ 0' = -a/N(0) . (50) 



Thirdly, one may borrow from Wilson's solution that the fixed point has So — 
7r/2. This expresses that the tightly bound spin singlet state formed of the 

impurity spin and conduction electron spin completely blocks the impurity 
site to other conduction electrons; this in turn implies maximal scattering 
and phase shift of 7r/2 for the effective scattering potential [184]. In other 
words, it is a strong-coupling fixed-point where one conduction electron state 
is pushed below the chemical potential in the vicinity of the impurity to form 
a singlet resonance with the impurity spin. One may now calculate all physical 
properties in term of a. In particular, one finds 6C^/C^ = 2a/{7rVN{0) and a 
similar expression for the enhancement of x,such that the Wilson ratio of 2. 



3.3 Ferromagnetic Kondo problem and the anisotropic Kondo problem 



The ferromagnetic Kondo problem provides us with the simplest example of 
SFL behavior. We will discuss this below after relating the problem to a general 
X-ray edge problem in which the connection to the so-called orthogonality 
catastrophe is clearer. As discussed in Sec. 2, orthogonality plays an important 
role in SFLs generally. 

we start with the anisotropic generalization of the Kondo Hamiltonian, which 
is the proper starting model for a perturbative scaling analysis [21,98], 



H^tY^ cl^Cja + [J± {S'^cliCkn + S cl^Ck'i) + 

{ij),<7 k,k' 

+J^S' (cl^Ck',^ - cliCk',i)] . (51) 

Long before the solution of the Kondo problem, perturbative Renormalization 
Group for the effective vertex coupling constants J± and as a function of 
temperature were obtained [21,98]. The scaling relation between them is found 
to be exact to all orders in the J's: 

{Jl - Jl) = const, (J+ = J_) . (52) 



In the flow diagram 18 we show the scaling trajectories for the anisotropic 
problem. In the antiferromagnetic regime the flows continuously veer towards 
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Fig. 18. Renormalization group flows for the Kondo problem, displaying the line of 
critical points for the "ferromagnetic" problem and the flow towards the fixed point 
J± = J* ^ for the "antiferromagnetic" problem. 

larger and larger {J±, J^) values; at the resulting fixed point ( J*±, J*z) ~ 
(oo, oo) singlets form between the local moment and the conduction electrons. 

The 'ferromagnetic' regime spans the region satisfying the inequalities < 
and \Jz\ > J±- Upon reducing the bandwidth the coupling parameters flow 
towards negative Jz values. Observe the line of fixed points on the negative Jz 
axis. Such a continuous line is also seen in the Kosterlitz-Thouless transition 
[139] of the two-dimensional XY model. Moreover, in both problems continu- 
ously varying exponents in physical properties are obtained along these lines 
(in fact, the renormalization group flow equations of the Kondo model for 
small couphng are mathematically identical to those for the XY model) . This 
is an instance of a zero temperature Quantum Critical line. The physics of the 
Quantum Critical line has to do with an "Orthogonality Catastrophe" which 
we describe next. Such orthogonalities are an important part of the physics of 
SFL's generally. 

3.4 Orthogonality Catastrophe 

As we saw in section 2, a Fermi-liquid description is appropriate so long as 
the spectrum retains a coherent single particle piece of finite weight Z > 0. So 
if by some miracle the evaluation of Z reduces to an overlap integral between 
two orthogonal wave functions then the system is a SFL. 

In the thermodynamic {N — > 00) limit, such a miracle is more generic than 
might appear at first sight. In fact, such an orthogonality catastrophe arises 
if the injection of an infinitely massive particle in more than one dimension 
produces an effective finite range scattering potential for the remaining N 
electrons [20] (see Sec. (4.9). Such an orthogonality is exact only in the ther- 
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modynamic limit: The single particle wave functions are not orthogonal. It 
is only the overlap between the ground state formed by their Slater determi- 
nantsp^which vanishes as tends to infinity. 



More quantitatively, if the injection of the additional particle produces an 
s-wave phase shift 60 for the single particle wave functions (all of them), 



then an explicit computation of the Slater determinants reveals that their 
overlap diminishes as 



Here \iPn) is the determinant Fermi sea wave function for particles and 
l?/'^) is the wavefunction of the system after undergoing a phase shift by the 
local perturbation produced by the injected electronP^ 

Quite generally such an orthogonality {Z = 0) arises also if two particle 
states of a system possess different quantum numbers and almost the same 
energy. These new quantum numbers might be associated with novel topolog- 
ical excitations. This is indeed the case in the Quantum Hall Liquid where 
new quantum numbers are associated with fractional charge excitations. The 
SFL properties of the interacting one-dimensional fermions (Section 4) may 
also be looked on as due to orthogonality. Often orthogonality has the effect of 
making a quantum many-body problem approach the behavior of a classical 
problem. This will be one of the leitmotifs in this review. We turn first to a 
problem where this orthogonality is well understood to lead to experimental 
consequences, although not at low energies. 

3.5 X-ray Edge singularities 

The term X-ray edge singularity is used for the line shape for absorption 
in metals by creating a hole in an atomic core-level and a particle in the 
conduction band above the chemical potential. In the non-interacting particle 
description of this process, the absorption starts at the threshold frequency 

The results hold true also for interacting fermions, at least when Fermi-liquid 
description is valid for both of the states 

19 Through the Friedel sum rule Sq/tt has a physical meaning; it is the charge that 
needs to be transported from infinity to the vicinity of the impurity in order to 
screen the local potential [78]. 




(53) 



(54) 
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Fig. 19. Absorption lineshape for transitions between a lower dispersionless level 
and a conduction band, (a) for zero interaction between the conduction electrons 
and local level, and (h) for small interaction. 

as sketched in Fig. 19. In this Fermi edge reflecting the density 

of unoccupied states in the conduction band is expected to be visible the 
spectrum. 

However, when a hole is generated in the lower level, the potential that the 
conduction electrons see is different. The relevant Hamiltonian is now 

U^e, (Sd - ^) + E ^Ack + ^ E V{k, k') (4c,, - {dU - , (55) 

where spin indices have been suppressed. The operators {d, d'^) annihilate or 
create holes in the core level, which is taken to be dispersionless. The flrst two 
terms in the Hamiltonian represent the unperturbed energies of the core hole 
and the free electrons. The last term depicts the screened Coulomb interaction 
between the conduction electrons and the hole in the core level. 

As a consequence of the interactions, the line shape is quite different. This is 
actually an exactly solvable problem [183]. There are two kinds of effects, (a) 
excitonic — the particle and the hole attract, leading to a shift of the edge 
and a sharpening of the edge singularity — and (h ) an orthogonality effect of 
the type just discussed above, which smoothens the edge irrespective of the 
sign of the interaction. This changes the absorption spectrum to that of Fig. 
V^(h) in the presence of interactions. The form of the singularity is [158,183] 

A{uj) ^{uj- uoy^^^^. (56) 

The exponent Sq /tt^ is a consequence of the orthogonality catastrophe overlap 
integral; the exponent (— 25o/7r) is due to the excitonic particle hole interac- 
tions. If the hole has finite mass we have a problem with recoil which is not 
exactly solvable, but we know the essential features of the solution. As we will 
discuss later in section 4.9, recoil removes the singularity in two and three di- 
mensions and the absorption edge acquires a characteristic width of the order 
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of the dispersion of the hole band. If the hole moves only in one dimension, 
the singularity is not removed. 

3.6 A Spinless Model with Finite Range Interactions 

A model, which is a generalization of the Ferromagnetic Kondo problem and in 
which the low-energy physics is dominated by the orthogonality catastrophe, 
is given by the following Hamiltonian : 

k,l y^k 

+ ig<7l7.,,-i)(^'^-i) . (57) 

The operators (7,7^) are the annihilation and creation operators of spinless 
conduction electrons with kinetic energy e^, while as before the d operators cre- 
ate or annihilate electrons in the localized level. The local chemical potential 
has been set to zero (e^ = 0) and the Hamiltonian is particle-hole symmetric. 
The new index / is an orbital angular momentum index (or a channel index). 
Hybridization conserves point-group symmetry, so the localized orbital hy- 
bridizes with only one channel {I — 0). By contrast, the impurity couples to 
all channels via the interaction Vi. As we are summing over all moments {k,k') 
this interaction is local. 

This problem may be mapped onto the anisotropic Kondo model [103]. Indeed 
the transformation 

2Vi ^ V2J,,i - 2nvF{V2 - l)Si,o , (58) 

n^Yl ^kcla,lCk,a,l + ^J±,o{^''sf + h.C.) + ^ Jz,lSzSz,l- (59) 
k,(T,l I 

Here a is short distance cutoff. In the resulting (anisotropic multi-channel) 

Kondo Hamiltonian the spin operators S and S/ portray charge excitations of 
the local orbital and conduction band; the spin index in the resulting Kondo 
Hamiltonian should now be regarded as a charge label. Physically, this map- 
ping is quite natural. The impurity may or may not have an electron, this is 



\/27ra ' 

produces 



43 



akin to having spin up or spin down. Similarly, the kinetic hybridization term 
transforms into a spin flip interaction term of the form (S'^rf^ + h.c). As Vi 
couples to the occupancy of the impurity site, we might anticipate to scale 
with Vi- The additional correction {—2ttvf{V^ — l)'5/,o) originates from the 
subtle transformation taking the original fermionic system into an effective 
spin model. 

This problem has been solved by renormalization group methods. See Fig. 20. 
But simple arguments based on the X-ray edge singularity, orthogonality and 
recoil give the correct qualitative physics. When t = 0, the problem is that 
of the X-ray edge Hamiltonian (with = 0). When t is finite, the charge 
at the impurity orbital fluctuates (the impurity site alternately empties and 
fills). This generates, in turn, a fluctuating potential. The X-ray absorption 
spectrum is the Fourier transform of the particle-hole pair correlator 

AH - (7t(i)d(t)cit(0)7(0)),. (60) 

This quantity, which is none other than the hybridization correlation function, 
should display the X-ray edge characteristics for large frequencies (a; > ^eff) 
where the effect of recoil is unimportant: 

A(u;)^A,{u;/Wy; 7 = -— + E 4' (61) 

TT ^ TT 

The threshold frequency Ag// is determined by the recoil energy. W is the 
bandwidth. The bare hybridization width Aq ~ t'^/W. The exponent in the 
singularity contains an exitonic shift {—2So/Tr) as well as an orthogonality 
contribution (Z^^^f/TT^). The recoil is cut off by Ae//. For co < Aeff the 
electron gas becomes insensitive to the change in the potential. As the X-ray 
edge singularity is cut off at a; = 0{Aeff), self-consistency implies that 

Ae// = A(a; = Ae//) . (62) 

This leads to the identification 

Ae// = W{Ao/Wf , P = 1/(1 - 7) , (63) 

so that for 

7 < 1 , Ae// ^0 asW ^oo . (64) 

For 7 < 1, a singular Fermi-liquid emerges in which the hybridization of the 
localized d-orbital with the electron gas scales to zero at zero frequency. The 
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actual value of 7 determines the singular properties at low energy or temper- 
ature. In the single channel problem such a scenario occurs if the potential Vq 
is sufficiently attractive. On mapping to the spin problem we find that this re- 
gion corresponds to the singular Fermi-liquid Ferromagnetic Kondo problem. 
The scaling of the hybridization to zero corresponds, in the spin-model, to 
J-t — > 0. In a renormalization group language, the flows will impinge on the 
line of fixed points (J^ = 0, < 0). In this regime, we recover, once again, a 
continuous set of exponents. 

If the number of channels is large enough, the orthogonality catastrophe as- 
sociated with the change in the number of particles on the impurity site is 
sufficiently strong to drive the hybridization to zero even for the case of re- 
pulsive interactions V or antiferromagnetic Jz^i- 

In the singular regime various correlation functions may be evaluated [103]. 
For instance, the Green's function of the localized impurity 



The orthogonality induced by the fiuctuation in the occupation number of the 
impurity site leads to the decay of the correlation function with a nonuniversal 
exponent. Because of this orthogonahty catastrophe the system behaves as a 
singular Fermi-liquid. In the vicinity of the t — fixed point, the self energy 
due to hybridization 

E ~ u^-^. (66) 



A line of critical points for 7 < 1 is found. This bears a resemblance to 
the Kosterhtz-Thouless phenomenon [139]. The analogue to the emergence 
of vortices in the Kosterlitz-Thouless transition are instantons — topological 
excitations which are built of a succession of spin flips in time on the impurity 
site. 



3.7 A model for Mixed- Valence Impurity 



We next consider a slightly more realistic model [190,225,229] 



k,a,l 



+ Yl ^kk'l Ud - ^) Y^^kalCk'al 
k,k',l ^ V<y 



(67) 



45 



Fig. 20. The Renormalization flow diagram for the model with finite range interac- 
tions according to [103]. The initial values are Jxfl = 0, Jl,o = 0.5, and Jx,i varies 
from to 1 in increments of 0.05. When J^^i becomes large enough (i.e. when V is 
large enough) the flow veers from the usual Kondo fixed point to a zero hybridization 
(•^±,0 = 0) singular Fermi liquid. 

In this model both spin and charge may be altered on the impurity site. This 
enhanced number of degrees of freedom implies that the states need to be 
specified by more quantum numbers. This also allows, a priori, for higher 
degeneracy. In the following screening channels, all of equal strength V, 
will be assumed. In the U —>■ oo limit the spectrum of the impurity site may 
be diagonalized. The two lowest states are 



t n\ in i\ 77 
Tj' U) = |U, 1), With energy = - 



4 2 ' 



Cl|0) = |cr, 0), with energy 




1/2 

(68) 



where in the bras, the first entry is the charge and spin of the impurity spin 
and the second one is the compensating charge in the screening channels. 
Other states are elevated by energies OiV / \/N). The states satisfy the Friedel 
screening sum rule by having a small phase shift ^ in each of the channels. 
The parameter V can be tuned to produce a degeneracy between the two states 
{E^ = E(^) — the mixed valence condition. The enhanced degeneracy produces 
Singular Fermi-liquid like behaviorf^. A perturbative calculation for small u 



'^^ There is a singularity only in the local charge response at the impurity, not in the 
magnetic response. In this respect, the results of ([229]) are not completely correct 
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yields a self-energy 



H [culncu + iusgn{u)] + 0[{u!hiu!y] . 



(69) 



As a speculative note we remark that this physics might be of relevance to 
quantum dot problems. Quantum dots are usually described in terms of the 
Anderson model. However, there are certainly other angular momentum chan- 
nels whereby the local charge on the dot and the external environment can 
interact. As the external potential in the leads is varied one is forced to pass 
through a potential in which this mixed valence condition must be satisfied. 
At this potential the aforementioned singular behavior should be observed. 

3.8 Multi-channel Kondo problem 

Blandin and Nozieres [185] invented the multichannel Kondo problem and 
gave convincing arguments for a local singular Fcrmi-liquid behavior of it. 
Since then it has been solved by a multitude of sophisticated methods. For an 
overview of these and of applications of the multichannel Kondo problem, we 
refer to [212]. 

The multichannel Kondo problem is the generalization of the Kondo problem 
to the case in which the impurity spin has arbitrary spin S and is coupled to 
n "channels" of conduction electrons. The Hamiltonian is 



Here (. is the channel index. Degeneracy, the key to SFL behavior, is enforced 
through equal antiferromagnetic couphng J > for all the channels. When 
the couplings to the various channels are not all the same, at low enough 
temperatures a crossover to local Fermi-liquid behavior in the channel with 
the largest always occurs [7]. This crossover temperature is in general quite 
large compared to Tk because channel asymmetry is a relevant perturbation 
about the symmetric fixed point. Therefore in comparing this theory with SFL 
behavior in experiments, one should ensure that one is above the crossover 
temperature. 

The simple Kondo problem is the case 2S — 1 — n. In this case at low 
temperatures a singlet state is formed of the impurity state and the conduction 
electron electrons in the appropriate channel. In the general multi-channel case 
in which 25* = n, the physics is essentially the same, since there are exactly 
the right number of conduction electron channels to compensate the impurity 
spin at low temperatures. Thus, at low temperatures an effective spin state 



n 

rimcK -tZ^l^C^a ^ja + JO ■ Cq GCq . 



(70) 
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Fig. 21. In a multichannel problem, a (Hund's rule coupled) spin S is compensated 
by total spin 5i (= n/2) in n different conduction electron channels leaving a net 
uncompensated spin as shown, or exact compensation or over-compensation when 
special properties may arise as discussed in the text. 

is formed again, and the properties of the compensated Kondo problem 2S = n 
at finite temperature is that of crossover from a weakly interacting problem 
above the Kondo temperature to a strongly interacting problem below Tx- 

The physics of the underscreened Kondo problem 2S > n is differentPn. In 
this case, there are not enough conduction electron spins to compensate the 
impurity spin. As a result, when the temperature is lowered and the effective 
coupling increases, only a partial compensation of the impurity spin occurs by 
conduction electrons with spin opposite to it. As depicted in Fig. 21, a net spin 
in the same direction as the impurity spin then remains at the impurity site. 
Since the conduction electrons with their spin in the same directions as the 
impurity spin can then still make virtual excitations by hopping on that site 
while the site is completely blocked for conduction electrons with opposite 
spin, a net ferromagnetic coupling remains between the remaining effective 
spin and the conduction electrons. As a result, the low temperature physics of 
the underscreened Kondo problem is that of the ordinary ferromagnetic Kondo 
problem. To be more precise, the approach to the fixed point is analogous to 
that in the ferromagnetic Kondo problem along the boundary = —J±, 
because the impurity must decouple (become pseudo-classical) at the fixed 
point. 

In the overscreened Kondo problem 2S < n, there are more channels than 
necessary to compensate the impurity spin. At low temperatures, all n chan- 
nels tend to compensate the impurity spin due to the Kondo effect. Channel 
democracy now causes an interesting problem. As the effective interaction 
J scales to stronger values, a local effective spin with direction opposite to 
the impurity spin results. This effective spin must have an effective antifer- 
romagnetic interaction with the conduction electrons, since now the virtual 
excitations of conduction electrons with spins opposite to the effective local 
spin lower their energy. This then gives a new Kondo problem with a new 

Since it is hard to imagine that the angular momentum states of the impurity are 
larger than of the conduction electron states about the impurity, such models may 
be regarded of purely theoretical interest. See however Sec (6.5) 
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Fig. 22. Flow diagram for the degenerate over-screened Kondo problem exhibiting 
a critical point. 

effective interaction, and so on. Of course, in reality one does not get a suc- 
cession of antiferromagnetic Kondo problems — the net effect is that a new 
stable finite-J fixed point appears. As sketched in Fig. 22, the renormalized 
effective interaction flows to this fixed point both from the strong-coupling as 
well as from the weak-coupling side. One can understand this intuitively from 
the above picture: if one starts with a large initial value of J, then in the next 
order of perturbation theory about it, the interaction is smaller, since in per- 
turbation theory the effective interaction due to virtual excitations decreases 
with increasing J. This means that J scales to smaller values. Likewise, if we 
start from small J, then initially J increases due to the Kondo scaling, but 
once J becomes sufficiently large, the first effect which tends to decrease J 
become more and more important. The finite-J fixed point leads to nontrivial 
exponents for the low-temperature behavior of quantities like the specific heat. 



For n = 2, the power law behavior on the right hand side is replaced by a In T 
term. 

Another way of thinking about the problem is in the spirit of Wilson's renor- 
malization groTip: Consider the problem of two channels interacting with a 
5* = 1/2 impurity. The conduction states can be expressed as linear combina- 
tions of concentric orbitals of conduction electrons around the impurity. These 
successive orbitals peak further and further away from the impurity. Consider 
first the exchange coupling of the orbitals in each of the two channels peaking 
at the impurity site. Each of them has S = 1/2. Only one linear combination 
of the two channels, call it red, can couple, while the other (blue) does not. So 
after the singlet with the impurity is formed, we are left with & S = 1/2, color 
blue problem. We must now consider the interaction of this effective impurity 
with the next orbital and so on. It is obvious that to any order we will be 
left with a spin 1/2 problem in a color. Conformal field theory methods first 
showed that the ground state is left with 1/2 In 2 impurity. A nice application 
of the bosonization method [83] identifies the red and blue above as linear 
combinations of the fermions in the two channels so that one is purely real, 
the other purely imaginary. The emergence of new types of particles — the 
Majorana fermions in this case — often happens at Singular Fermi-liquids. 

As a detailed calculation confirms [7], the Ji = J2 fixed point is unstable, 
and the flow is like sketched in Fig. 24. This means that the Ji = J2 fixed 




n-2 



(71) 
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Fig. 23. Effective sfiells in Energy space in Wilson's metliod. Tlie first sliell integrates 
over a fraction A of the top of the band, the next shell A of the rest, and so on. In 
the two-channel 5 = 1/2 problem, a S = 1/2 effective impurity is left at every stage 
of interpretation. 

point is a Quantum Critical Point: in the T — J1/J2 phase diagram, there is a 
critical point at T = 0, Ji/ J2 = 1. Moreover, it confirms that asymmetry in the 
couplings is a relevant perturbation, so that the SFL behavior is unstable to 
any introduction of differences between the couplings to the different channels. 
The crossover temperature below which two-channel behavior is replaced 
by the approach to the Kondo fixed point is [7] 



where J = ( + J2)/2. 

The overscreened Kondo problem illustrates that the SFL behavior is associ- 
ated with the occurrence of degeneracy: the critical point requires degeneracy 
of the two orthogonal channels. 

An interesting application of the two-channel Kondo problem is obtained by 
considering the spin label to play the role of the channel index, while the 
Kondo coupling is in the orbital angular momentum or crystal field states 
for impurities at symmetry sites in crystals [135,69,212]. Another possibility 
that has been considered is that of scattering of conduction electrons of two- 
level tunneling centers with different angular momentum states [287]. In this 
case, the tunneling model translates into a model with x and z coupling only, 
but this model flows towards a Kondo-type model with equal x and y spin 
coupling. For both types of proposed applications, one has to worry about the 
breaking of the symmetry, and about the question of how dilute the system 
has to be for a Kondo-type model to be realistic. Interesting results in the 
tunneling conductance of two metals through a narrow constriction, shown in 
Fig. 25 appear to bear resemblance to the properties expected of a degenerate 




(72) 
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Fig. 24. Flow in temperature-anisotropy plane for the two-channel Kondo problem 
with coupling constant Ji and J2. 

two-channel Kondo effect [205,70], but this interpretation is not undisputed 
[274,206]. Applications to impurities in heavy fermions will be brieffy discussed 
in section 6.4. 



3.9 The Two-Kondo-Impurities Problem 

In a metal with a finite concentration of magnetic impurities, an important 
question is what the (weak) interaction between the magnetic impurities does 
to the Kondo physics — that the effect might be substantial is already clear 
from the fact that the Kondo effect is seen in logarithmic corrections about 
the high-temperature local moment fixed point while the RKKY interaction 
between the moments mediated by the conduction electrons occur as a power 
law correction!^. Stated simply ~ j^jyj exp(— ^^^^j ) while the RKKY in- 
teraction I ~ J'^N{0). The existence of mixed- valence and heavy fermion 
metals makes this much more than an academic question. The question of the 
competition between these two effects, and in particular whether long range 
magnetic order can arise, was first posed by Varma [248] and by Doniach [76], 
who gave the obvious answer that RKKY interactions will be ineffective only 
when the Kondo temperature below which the local spin at each impurity is 
zero is much larger than the RKKY interaction I. Considering that JN{0) is 

Stated technically, the Kondo Hamiltonian is a marginal operator while the 
RKKY operator is a relevant operator about the local moment fixed point: In a 
perturbation calculation the interaction produces corrections of 0{1/T) compared 
to a ln(r) correction of the Kondo effect. 
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Fig. 25. (a) Differential conductance G{V, T) as a function of voltage V in measure- 
ments on metal point contacts by Ralph et al. [205], for various temperatures ranging 
from 0.4 K (bottom curve) to 5.6 K (upper curve). Note the VV type behavior de- 
veloping as the temperature decreases, (b) The zero bias {V = 0) conductance as a 
function of temperature for three different samples shows a G(0, T) — G{0, 0) ~ s/T 
behavior. The scaling behavior as a function of voltage and temperature is consistent 
with two-channel Kondo behavior [205]. 

usually » 1, this is unlikely for 5" = 1/2 problems, but for large S, as encoun- 
tered typically in rare earths and actinides, it is possible in some cases. But 
the vast majority of rare earths and actinide compounds show magnetic order 
and no Heavy fermion magnetic behavior. 

The two-Kondo-impurity problem serves as a first step to understanding some 
of the physics of heavy fermions. We will however only consider the S — 1/2 
problem and work with unrealistic JN{0) so that the competition between 
the Kondo effect and RKKY is possible. 

In this subsection we will summarize the results [128,129,8,9,229] for the two 

Kondo impurity problem; like the models we discussed above, this one also has 
a quantum critical point at which SFL behavior is found. However as for other 
impurity problems, an unrealistic symmetry must be assumed for a QCP and 
attendant SFL behavior. 
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The two-Kondo-impurity Hamiltonian is defined as 

^ = ^ E icfca + J [Si ■ a(rO + S2 • a{r2)] . (73) 



k.a 



In this form, the problem has a symmetry with respect to the midpoint be- 
tween the two impurity sites ri and r2, and hence one can define even (e) 
and odd (0) parity states relative to this point. In the approximation that 
the A;-dependence of the couphngs is neglected [128] the two Kondo impurity 
Hamiltonian can then be transformed to 



'H — 'Hnon— interacting "I" (Si + S2) • JeC][./g(TCfee + Jo(^\'o^^ko 

kk' 

+ Y.Jm (Si - S2) • cl^acko + c^^acke ■ (74) 



kk' 



The couphng between the spin and orbital channels generates an effective 
RKKY interaction 

Ti-RKKY — Io{Je, Jo, Jm) Si • S2 (75) 



between the two impurity spins, with /q = 21n2(J| + — 2J^) ior t — 1. 

The main results of a numerical Wilson-type renormalization group treatment 
of this model are the following: 

(i) For ferromagnetic coupling /q > or for a small antiferromagnetic coupling 
Iq > —2.24:Tk, where Tk is the Kondo temperature of the single-impurity 
problem, one finds that there is a Kondo effect with 

(Si • S2) 7^ , (76) 



unless Iq is very small, \Io/Tk\ ^ 1. Since for uncorrelated impurity spins 
(Si • S2) = 0, (76) expresses that although in the RG language the RKKY 
interaction is an irrelevant perturbation, it is quite important in calculating 
physical properties due to large "corrections to scahng" . Another feature of the 
solution in this regime is the fact that the phase shift is 7r/2 in both channels. 
This means that at the fixed point the even-parity channel and the odd-parity 
channel have independent Kondo effects, each one having one electron pushed 
below the chemical potential in the Kondo resonance. As discussed below, this 
is due to particle- hole symmetry assumed in the model — without it, only the 
sum of the phase shifts in the two-channels is fixed. 

(ii) There is no Kondo effect for Jq < —2.24 T^. In this case, the coupling 
between the impurities is so strong that the impurities form a singlet among 
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themselves and decouple from the conduction electrons. There is no phase shift 
at the fixed point. Also, in this case the total spin Stot = 0, but the impurity 
spins become only singlet like, (Si ■ S2) ~ —3/4 for very strong coupling, 
/o <S —2.24 Tk- So again there are important "corrections to scahng". 

(Hi) The point Jq = —2.24 is a true critical point, at which the staggered 
susceptibility ((Si — S2)^)/T diverges. Moreover, at this point the specific heat 
has a logarithmic correction to the linear T dependence, ~ TlnT, while 
the impurity spin correlation function (Si • S2) becomes equal to —1/4 at this 
value. 



Although the approximate Hamiltonian (74) has a true quantum critical point 
with associated SFL behavior, we stress that the analysis shows that this crit- 
ical point is destroyed by any cj^gCj^o coupling. A coupling of this type is not 
particle-hole symmetric; since the approximate Hamiltonian (74) is particle- 
hole symmetric, these terms are not generated under the renormalization 
group flow for (74). The physical two-Kondo-impurity problem (73), on the 
other hand, is not particle-hole symmetric. Therefore, the physical two Kondo 
impurity problem does not have a true quantum critical point — in other 
words, when the two-Kondo-impurity Hamiltonian (73) is approximated by 
(74) by ignoring k-independent interactions, relevant terms which destroy the 
quantum critical point of the latter Hamiltonian are also dropped. This has 
been verified in explicit analysis keeping the symmetry breaking terms [229] . 

An illuminating way to understand the result for the two-impurity Kondo 
problem, is to note that the Hamiltonian can be written in the following form: 



n 
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Hr 







n°s=i 



(77) 



where S = Si -|- S2 is the total impurity spin. 

In this representation, the Hamiltonian Timix couples the S* = and S = 1 
state, and the following interpretation naturally emerges: for large antiferro- 
magnetic values of Iq, Tls=o is lower in energy than Tis=i, while for large 
ferromagnetic coupling Iq, the converse is true. The two- Kondo impurity cou- 
pling can thus be viewed as one in which by changing Iq, we can tune the 
relative importance of the upper left block and the lower right block of the 
Hamiltonian. In general, the two types of states are mixed by Hmix, but at the 
fixed point "Hj^j^. 0. This impfies that there is a critical value of Iq/Tk where 
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the S = and S = 1 states are degenerate, and where SFL behavior occurs. 
At this critical value, the impurity spin is a linear combination of a singlet and 
triplet state with (Si • S2) = — 1/4 (i.e. a value in between the singlet value 
-3/4 and the triplet value 1/4) and the singular low-energy fluctuations give 
rise to the anomalous specific heat behavior. 

Within this scenario, the fact that the susceptibility x is divergent at the crit- 
ical point signals that a term H • (Si — S2) lifts the spin degeneracy. Moreover, 
the leading irrelevant operators about the fixed point are all divergent at the 
critical point — of course, this just refiects the breakdown of the Fermi-liquid 
description. 

The reason for T~lmix — ^ is as follows. It arises from the last term in (74) which 
is not particle-hole symmetric because under even-odd interchange both the 
spin term and the fermion electron terms change sign. At the Kondo fixed 
point, the leading operators must all be biquadratic in fermions. An Timix in 
that case would be of the form cl^^Ck'ao and and such a term by itself would 
break particle hole symmetry. 

In the two-Kondo impurity problem, one encounters again the essentials of de- 
generacy for quantum critical points and the need for (unphysical) constraints 
to preserve the singularity. Once again new types of quantum numbers can be 
invoked in the excitations about the QCP. 

From the point of view of understanding actual phenomena for problems with 
a moderate concentration of impurities or in reference to heavy fermion com- 
pounds, the importance of the solution to the two-Kondo impurity problem is 
the large correction to scaling found in the Wilson-type solution away from the 
special symmetries required to have a QCP. These survive quite generally and 
must be taken into account in constructing low-energy effective Hamiltonians 
in physical situations. 



4 SFL behavior for interacting fermions in one dimension 

A variety of elegant mathematical techniques, including exact solutions in 
certain nontrivial limits, have been employed to analyze the problem of inter- 
acting electrons in one dimension. We discuss the essential features of these 
models here and their SFL properties. We also discuss whether the methods 
and the results can be extended to higher dimensions. 

Some cases where the models solved are experimentally realized include the 
edge states of Quantum Hall liquids and quasi-one dimensional organic and 
inorganic compounds [81,82]. In the latter case, the asymptotic low energy 
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properties are, however, unlikely to be those of the ID models because of 
the inevitable coupling to the other dimensions which proves to be a relevant 
perturbation. Nonetheless, data on carbon nanotubes [73] discussed in sec- 
tion 4.11 show clear evidence of one- dimensional interacting electron physics. 
Several one-dimensional spin chains problems can also be transformed into 
problems of one- dimensional fermions [210]. 

Our aim in this chapter is to highlight the physics leading to the Singular 
Fermi Liquid behavior and to present the most important results obtained. 
Detailed reviews of the technical steps in the various solutions as well as 
numerical calculations may be found in [109,231]. As pointed out already in 
section 2.6, logarithmic singularities appear in the second order perturbative 
calculation of the vertices in the ID problem (both in the particle-hole and 
particle-particle channels). We present the T = phase diagram as obtained 
by summing those singularities by perturbative RG. This will be followed by a 
discussion of the results of the exact solution of the Tomonaga-Luttinger and 
the more general model for special values of the coupling constants (along the 
so called "Luther-Emery" line [155]). 

One special feature of ID physics is that the low energy excitations can be de- 
scribed by either fermions or bosons. The bosonic description of the Tomonaga- 
Luttinger model is especially attractive and will be presented below. 

A distinctive feature of one dimensional physics is that single-particle as well 
as multiple-particle correlation functions are expressible in terms of indepen- 
dent charge and spin excitations, which, in general, propagate with different 
velocitiesp^ This feature has been shown to arise due to extra conservation 
laws in one dimensions [50,167]. As we shall discuss, an extension of charge- 
spin separation to higher dimensions is unlikely because there are no such 
conservation laws. 

The one dimensional singularities may be also seen as a manifestation of the 
orthogonality catastrophe [20] that we discussed in section 3.4. We shall see 
that this feature disappears in higher dimension due to the effects of recoil. 



Even a Fermi-liquid displays distinct energy scales for charge and spin (particle- 
hole) fluctuations because of the difference in the Landau parameters in the spin- 
symmetric and spin-antisymmetric channels. The phrase spin-charge separation 
should therfore be reserved for situations, as in one-dimension, where the single- 
particle excitations separate into objects carrying charge alone and which carry 
spin alone. 
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4-1 The One Dimensional Electron Gas 

The Hamiltonian describing the low energy dynamics reads 

^ T~^0 ~l~ forward ~l~ "^backward 

+ n 

Umklapp • 

(78) 

In (78) Tio is the free electron Hamiltonian 

T^o = J2 ^kclck ■ (79) 

k 

In one dimension the particles move in two branches, either to the left or the 
right. Define the charge density operators pr and the spin density operators 
Sr for the two branches, r = ±, 

Pr = ^a=±li^l,aA,a , = ^^a,a'^laT^,a'i^r,a' , (80) 

where is a Pauli-matrix. The Fourier components of the particle density 
operators read 

PrAl) = Yl 4,a,k+qCr,a,k = ■ (§1) 

k 



These bilinears may be explicitly shown to obey Bose commutation relastions. 
If we linearize about the two Fermi-points, the energy of a particle-hole pair 
created by Pr,aiq), 

er,k+q - er,k = TVpq , (82) 



is independent of k in one dimension (provided the band cutoff is taken to in- 
finity). States created by Pr,a{Q) are linear combinations of individual electron- 
hole excitations all of which have the same energy and are therefore eigenstates 
of Hq. It follows that for q > 0, Pr=+,a{q) [Pr=-,a{q)] is a raising [lowering] op- 
erator. Consequently TIq may be expressed in terms of the density operators: 

Ho = YprA^q)PrA-'>^q) ■ (§3) 

^ r=± q>0 

Upon separating the densities on a given branch into charge and spin pieces, 

Pra{x) = ^[Prix) + <jS',{x)] , (84) 
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the free Hamiltonian may be expressed as a sum in the spin and charge degrees 
of freedom: 

^o = E^ok]+E^om. (85) 

r r 

It follows that in the noninteracting problem spin and charge have identi- 
cal dynamics and propagate in unison. Once interactions are introduced, the 
electron will "fractionaliz^' and spin and charge dynamics will, in general, 
differ. 

We now discuss the various interaction terms in the Hamiltonian. The piece 
describing the forward scattering events {kp, c; —kp, cr') — > {kp, cr; —kp, cr') 
and {kp, cr; kp, a') — > {kp, a; kp, a') reads 

'H forward = 'H2 + 'H4 , (86) 

with 

^2 = y E E 9r'p+A<l)P-A-<l) , (87) 

H4^^Y.129f[p+AQ)p+A-Q)+P-AQ)p-,^'i-Q)] ■ (88) 

9 <7cr' 

The operators pr^a involve a creation and annihilation operator on the same 
branch. We also define operators which involve operators on opposite 
branches, 

Pa = Yl 4,a,kC-r,a,k+q ■ (89) 
k 

In terms of these the Hamiltonian describing backscattering interactions (the 
scattering event {+kp, a; —kp, a') {—kp, cr; kp, cr') and its reverse) becomes 

Uhackwards = '^1 = i/l E E Ql"' Pt (q) Pa'(-(l) ■ (90) 

g aa' 

and the Umklapp term reads 

^3 = ^T.T.9r'[pU<i)pt'i-<i) + Pai<i)p-A-<i)] ■ (91) 

g aa' 

In all these expressions, the coupling constants may be spin dependent: 

gf = 9i\\^aa' + 9i±Sa-a' ■ (92) 
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Fig. 26. Pictorial representation of the low-energy interaction terms in the 
one-dimensional problem. After [169]. The and "-" points are a shorthand 
for the two Fermi points k = kp and {—kp) respectively. 

As the terms Hi\\ and H2\\ describe the same process we may set gi\\ =0 with 
no loss of generality. Umklapp processes are important only when Akp is a 
reciprocal lattice vector so that all scattering particles may be near the Fermi 
points. If the one dimensional electronic system is sufficiently incommensurate 
then these processes may be neglected. 

The condition for spin rotation invariance [7^, 5] =0 reads 

921. - gi± = 5-211 - 5111 • (93) 



A guide to the behavior of the Hamiltonian (78) is obtained from the pertur- 
bative Renormalization Group flow equations. In one dimension, there are no 
truly ordered phases, but at T = there is algebraic long range order. One 
may thus determine a "phase diagram" according to which susceptibilities di- 
verge as T — >• 0: the one associated with singlet superconductivity (SS), triplet 
superconductivity (TS), charge density wave (CDW) state at 2kp, and spin 
density wave {SDW) state at 2kp. The expressions for these susceptibilities 
are given in Eq. (Ill) below, and the resulting "phase diagram" is shown in 
Fig. 27. 



4.2 The Tomonaga-Luttinger Model 



With forward scattering alone in (78) and after linearizing the kinetic energy 
we obtain the Tomonaga-Luttinger (T-L)model. In terms of fermion field op- 
erators ipr,a{,x) and ipl^{x), and density and spin operators pr{x) and S^{x) in 
real space, the T-L model is 



T-L 



dx 



-Vp J2 r'^l,J9x1pr,a + 92,cPr{x)p- 



\X 



+ X! g4:,cPr{x)pr{x) + 2 ^ g2,sSl{x)S%{x) + gi^sSl{x)Sl{x) 



(94) 
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Fig. 27. Phase diagram for one-dimensional interacting fermions in the {Kc,gi) 
plane with = in Eqn. (92) and Kc given by Eq.l02 . SDW and {CDW) 
in the upper left quadrant indicates that both the spin density wave susceptibility 
and the charge density wave susceptibility diverge as T — > 0, but that the spin 
density wave susceptibility diverges a factor In^T faster than the charge density 
wave susceptibility — see Eq. (111). Other sectors are labeled accordingly. Prom 
[218]. 



where 



9i\\ + 9i± 



9i\\ - 9i± 



(95) 



Note that the g2,s term is the only term which break SU (2) spin symmetry. 
The T-L model is exactly solvable. After all, as previously noted, the (Dirac- 
like) kinetic energy Hamiltonian Tio is also quadratic in the density operators. 
So the Hamiltonian is readily diagonalized by a Bogoliubov transformation 
whereupon the Hamiltonian becomes a sum of two independent (harmonic) 
parts describing noninter acting charge and spin density waves: the charge and 
spin density waves are the collective eigenmodes of the system. 

The simplest way to describe the TL model and to explicitly track down these 
collective modes is via the bosonization of the electronic degrees of freedom. 
The reader should be warned that many different conventions abound in the 
literature. The bosonic representation of the fermionic fields proceeds by writ- 
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ing [109,180,210] 



— lim 



V2 



(96) 



where a is a short distance regulator. ^ra{x) satisfies 



(97) 



The "Klein factors" are chosen such that the proper fermionic anticommu- 
tation relations are reproduced. The exponential envelope exp[i^ra{x)] repre- 
sents the slow bosonic collective degrees of freedom which dress the rapidly 
oscillating part Fra explikpx] describing the energetic particle excitations near 
the Fermi points. 

$ may be written in terms of the bosonic fields 0c,s and their conjugate mo- 
menta dxOc s 



(98) 



In terms of the new variables, the charge and spin densities read 

P(^) = E Pr{^) = {l^-^c , S'{X) = 5;(X) = \[^dxCt>s . (99) 



In the (^c,s) 0c,s) representation the Tomonaga-Luttinger Hamiltonian becomes 
a sum of two decoupled oscillators describing the gapless charge and spin 
density wave eigenmodes 



T-L 



v=c,s 



i^.(9A)' + 



{dx<Pu) 
K„ 



T-L I njT-L 



(100) 



The velocities of the collective charge and spin modes are easily read of by 
analogy to a harmonic string: 



c,s c,s 

+ — ) - K—, 

TT TT 



(101) 



Likewise, the moduli determining the power-law decay of the correlations are 

(102) 



^J t^vf + gT + gT ' 
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In section 4.7 we shall show how the above expressions for the spin and charge 
density wave velocities follow simply from the conservation of left and right 
moving particles in the T-L model. 

As previously noted, the charge and spin velocities are degenerate in the non- 
interacting model. When interactions arc introduced, the charge and spin ve- 
locities {vc and Vs) as well as the energy to create spin and charge excitations 
{{vs/Kg and Vc/ Kc respectively) become different. The charge constant Kc is 
less than 1 for repulsive interactions, which elevates the energy of the charge 
excitations, while Kc is greater than 1 for attractive interactions. 

4-3 One-particle Spectral functions 

The zero temperature spectral functions of the T-L model are 

A{k, u)^{cu- v,{k - A;ir))2>-V2|^ _ y^(k - kpW^-^^^ {v, > Vs) 
A{k,iu) ^{iu- v,{k - kF)y'~^^^\u} - v,{k - A;f)|2t--^/2 ^iQ3^ 

These spectral functions are sketched in Fig. 28 for the case g2 — and in 
Fig. 29 for the general case. Notice that the quasiparticle pole in A{k, cu) in a 
Landau Fermi-liquid is smeared out to produce a branch cut extending from 
the spin mode excitation energy to the charge mode excitation energy. These 
branch cuts split into two in an applied magnetic field, see Fig. 30. These 
results are exact for small cu and small \k — kp]- Another manifestation of the 
SFL behavior is the behavior of the momentum distribution function 

^fc ~ ''T'kp ~ const X sign(A; — kF)\k — , (104) 

where 

^c,s^I{Kc,s + K-^-2) . (105) 

In contrast to a Fermi liquid, the expression for does not exhibit a step-like 
discontinuity at the Fermi points. Note that the exponent 27c is non-universal 
(as usual, an outcome of a line of critical points). 

The single particle density of states 

N{u;) , (106) 

vanishes at the Fermi surface. 



62 















w 


{vck - w) 



Vfk Vck 



Fig. 28. The zero temperature spectral function A^kju) = lm{G^^^i{k,u;)} as a 
function of 00 for the case (52 = 0, g4 ^ 0) — the "one-branch Luttinger hquid" 
in which the spin and charge velocities differ but for which the the correlation 
exponents keep their canonical Kc = 1 according to (102). In the figure Vc > Vg and 
A; > are assumed. Note the inverse square root singularities. This is a consequence 
of iCc = 1 which makes 7c = 0. After Voit [262]. 




—Vck —Vgk Vgk Vck 



Fig. 29. The generic (§2 7^ 0, 54 7^ 0) zero temperature spectral function. Note the 
broader range of non-trivial singularities near Vgk and zizVck. Here both the effect of 
spin-charge velocity difference and the emergence of non-trivial exponents is visible. 
After Voit [262] 

4.4 Thermodynamics 

The contributions of the independent charge and spin modes appears inde- 
pendently in most physical quantities. 

The specific heat coefficient is found to be 




(107) 
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Fig. 30. The energy distribution curve (the spectrum A(k, uj) at fixed A;) as a function 
of u) in the presence of a magnetic field. The dashed fine is the zero-field result of Fig. 
28. The magnitude of the Zeeman splitting is enhanced with respect to {vc — Vs)k 
for clarity. Prom Rabello and Si [204]. 

where 7 = 7o for the non-interacting system. 

The spin susceptibihty and the compressibility are readily computed 

XO = Vf/Vs , k/kq^ VpKc/Vc (108) 



where xo and kq are the susceptibility and compressibility of the noninteracting 
gas. The Wilson ratio, already encountered in our discussion of the Kondo 
problem in section 3.2, 

= = (109) 

7/70 Vc + Vs 



deviates from its Fermi Liquid value of unity by an amount dependent on the 
relative separation between the spin and charge velocities. 



4-5 Correlation Functions 



Since the Hamiltonian is separable in charge and spin and as tjj is a product 
of independent charge and spin degree of freedom all real space correlation 
functions are products of independent charge and spin factors. We will only 
focus on the most important correlation functions in the illustrative examples 
below, and refer for a summary of the various exact expressions to [188]. 

The most important feature of the large distance behavior of the charge and 
spin correlators is their algebraic decay at zero temperature: 
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,,,,,,, Kc „ COs(2kF(x — x')) ^ _o/2 I 

{p{x)p{x')) ~ '-—^ + ^i,c I \|,^^ In X - 



X 



cos{AkF{x - x')) 

/^/ V ^/ /vx 1 „ cos(2A;f(x — x')) , 1/9, /I 

^ X • 5 ^ ^, + Si,,^^^— ^ lnV2 x - + ... 

(7r(x — x'))^ \x — x'\^+^'= 

at asymptotically long distances and Kg = 1. For not too repulsive interac- 
tions, so that Kc < 1, the 2kF fluctuations are dominant. We have previously 
seen that such a CDW/SDW instability may arise due to the special 2kF 
nesting wavevector in one dimension. The AkF modulations are due to Umk- 
lapp scattering. The amplitudes {Bid {Bi,s} are non-universal while the 
exponents are determined by the stiffness of the free charge and spin fields. 

While the above expressions are for Kg — 1, in the general case Kg 1 the 
spin correlator decays asymptotically with the exponent {Kg + K^). 

At nonzero temperatures, it is found that the Fourier transforms of these 
correlation functions scale as 



Xss^T^'^-'llnTl-'/' , XTS ^ T^'^-'llnTl'/' . (Ill) 

An example of the spectral density as a function at finite temperatures is 
shown in Fig. 31. 

The results above lead in principle to clear experimental signatures. X-rays, 
which couple to the charge density waves, should peak at low temperatures 
with intensities given by 

hk, ~ hk, - T^^^-' . (112) 



The NMR probe couples to the spin degrees of freedom and the theoretically 
computed nuclear relaxation time scales as 



Ti ~ T-^^ . (113) 



4-6 The Luther-Emery Model 



The Luther Emery model extends the Tomonaga-Luttinger Hamiltonian by 
including the backward scattering interactions parametrized by Tii, which 
scatter from (+A;f, ct; —kF, c') to {—kp, o", kF, cr') and vice versa. The Umklapp 
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Fig. 31. Left panel: Momentum distribution curves at u; = (i.e., the spectrum at 
fixed a; = as a function of k) for a spin rotationally invariant Tomonaga-Luttinger 
liquid, plotted as a function of v^k/T. Right panels: Energy distribution curves at 
k = (the spectrum at fixed k = 0) as a function of lo/T. In both panels, Vc/vs = 3 
and 7c = in (a), 'jc = 0.25 in (b), and 7c = 0.5 in (c). From Orgad [188]. 

processes, important (at T = 0) only when Akp is a reciprocal lattice, continue 
to be discarded. 

The backward scattering term 

ni= dx giYl i^l,a=+l^-r,a=-lA,a=-l^-r,a=l (114) 
r=±l 

written in terms of the bosonic variables introduces a nontrivial sine-Gordon 
like interaction: 



with rescaled values of the spin and charge velocities and stiffness constants. 
When gi > (repulsive interactions), gi is renormalized to zero in the long 
wavelength hmit. Since along the RG flow trajectories Kg — l^ gi/ {'kVs), this 
means that Kg renormalizes to f . The physics corresponding to this case is in 
the Tomonaga-Luttinger model that we just discussed. 

When the backscattering interactions are attractive, gi < 0, spin excitations 



2 , (5x0.^' 



+ -^||y^cos(V8^0,) (115) 
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are no longer gapless, and a spin gap of magnitude 



As ~ — 
a 



9i 



l/{2-2Ks 



(116) 



is dynamically generated. The attractive backscattering leads to the formation 
of bound particle-hole pairs which form a CDW. The spin correlation length 
is then finite, 



6 = ^. 



;ii7) 



Gaps in the charge spectrum develop in the presence of Umklapp scattering. 
In the spin gapped phase the Hamiltonian can be conveniently expressed in 
terms of new refermionized spin fields ^r(^)- 



= Frexp[-tJ7T/2{es - 2r0,)] 



:ii8) 



Luther and Emery observed that at the point Kg = 1/2 the Hamiltonian in 
terms of these new spin fields becomes that of noninteracting free fermions 
having a mass gap = gi/{2na). 

Hs^ [dxYl [-ivsr'ifld^'^r + A,*]:*-,] , (119) 

r=±l 



leading to the spin excitation spectrum 

Es = sj^l + viik-kpf . (120) 



^.7 Spin charge separation 



As in many other physical problems, the availability of an exact solution to 
the one-dimensional electron gas problem is intimately linked to the existence 
of additional conservation laws or symmetries. One may attack the Luttinger 
Liquid problem by looking for its symmetries. 

The U{1)l ® U{1)r symmetry in the absence of Umklapp scattering may be 
exhibited by considering the effect of the separate left and right rotations by 
angles on the fermion variables: 

i^La{x, t) ^ e'^^i^Laix, t) , 'iI:r,{x, t) ^ e^^^^'iia (a^, t) . (121) 
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All of the currents are trivially invariant under this transformation as the 
ip"^ fields transform with opposite phases. Physically this corresponds to the 
conservation of the number and net spin of left and right moving particles. 

As discussed by Metzner and Di Castro [167], these separate conservation laws 
for the left and right moving charge and spin currents lead to Ward identities 
which enable the computation of the single particle correlation functions. 

In the absence of Umklapp scattering, spin and charge are conserved about 
each individual Fermi point. The net total charge density p = p+ + p- and 
charge density asymmetry p = p+ — p_ in the Tomonaga-Luttinger Hamilto- 
nian satisfy the continuity equations 

drp = [n, p] = -qj , drP = [n, p] = -qj , (122) 

where 

= Uc[p+ - P-] , j = Uc[p+ + P-] , (123) 

and where the velocities are given by 

.91-92 ~ ,94 + 92 /10/l^ 
Uc — vf-\ , Uc — vf-\ ■ (124) 

TT TT 



These results follow straightforwardly from the form of Ti.^ ^ in combination 
with the fact that the only nonzero commutator is 

[Pr,a{Q)^Pr',a'{-q')] ^ Sqq'^rr'^aa/ (^^^ ■ (125) 



Let us illustrate simply how many of the results derived via bosonization 
may also be directly computed by employing these conservation laws. The 

existence of gapless charge modes is a direct consequence of the right-left 
charge conservation laws. The two first-order continuity equations given above 
lead to 

[dl + UcUcq'']p^O , (126) 

from which we can read off a linear charge dispersion mode 

u;^Vc\q\ (127) 



with velocity Vc = in agreement with the earlier result (101). Thus, 

collective charge excitations propagate with a velocity Vc- A similar relation 
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may be found for the spin velocity Vg which in general is different from Vc- 
This spin-charge separation also becomes clear from the explicit form of the 
expectation values of the charge and spin densities: 



{0\^r{xo)pr{x,t)^jJl{xo)\0) = S{x - Xq - TVct) , 

{0\il)r{xo)S^{x,t)il)l{xo)\0) ^ S{x -Xq- rVgi) , 



(128) 



where |0) denotes the ground state. 

The separate Right-Left conservation laws cease to hold if (back-scattering) 
impurities are present. Accordingly, as shown by Giamarchi and Schulz, [104] 
spin-charge separation then no longer holds. 



4-8 Spin-charge Separation in more than one- dimension? 

Spin-Charge separation in one dimension requires extra conservation laws. 
Can something analogous happen in more than one-dimension? No extra con- 
servation laws are discernible in the generic Hamiltonians in two-dimension, 
although such Hamiltonians can doubtless be constructed. Are there condi- 
tions in which generic Hamiltonians become dynamically equivalent to such 
special Hamiltonians (because the unwelcome operators arc "irrelevant"). No 
definite answers to these questions are known. In section (5.2)and later in this 
section we shall briefly review some interesting attempts towards spin-charge 
separation in higher dimensions. First we present qualitative arguments point- 
ing out the difficulty in this quest. 

There is a simple caricature given by Schulz [218] for qualitatively visualizing 
charge-spin separation for a special one-dimensional case: the C/ — > oo Hub- 
bard model. This model is characterized (at half-filling) by algebraic decay of 
Spin-Density correlations, which at short distances appear as almost antifer- 
romagnetic alignments of spins. Let us track the motion of a hole introduced 
into an antiferromagnetically ordered chain. The hole is subject to only the 
lattice kinetic term which enables it to move by swapping with a nearby spin. 

An initial configuration will be 



(129) 



After one move the configuration is 



(130) 
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After two additional moves to the left the configuration reads 

...uommu- 



(131) 



Thus the initial hole surrounded by two spins of the same polarization has 
broken into a charge excitation ("holon" or "chargon" — a hole surrounded 
by antiferromagnetically aligned spins) and a spin excitation ( "spinon" ) com- 
posed of two consecutive parallel spins in an antiferromagnetic environment. 
The statistics of the localized spinous and holons in this model must be such 
that their product is fermionic. 

The feasibility of well defined spin and charge excitations hinges on the com- 
muting nature of the right and left kinetic (hopping) operators TjnghtjTLeft 
which are the inverse each other. Any general term of the form 

{TLeft) ^ {Trnght) ^ {TLeft) ^ {TRight) ^ 

= {n^ftTmgntf-Tlf,''- = I X (132) 

where Np>x = J2i We have just shown that terms of the form T^^ lead to 
a cartoon of the sort depicted above which gives rise to spin-charge separation 
and therefore our result holds for the general perturbative term. The proof 
of spin-charge separation for the one dimensional electron gas rests on the 
existence of separate conservation laws for the left and right moving fronts, as 
a result of the fact that the operators Tmght and Ti^eft commute. 

Such a simple "proof" cannot be extended to higher dimensions. In higher 
dimensions this suggestive illustration for spin-charge separation is made im- 
possible by the noncommuting (frustrating) nature of the permutation op- 
erators TupjTcown, TRight, Tieft, ■■■ Moreover, even if the exchange operators 
commuted we would be left with terms of the form T^^fC'T^p'' which when 
acting on the single hole state will no longer give rise to states that may be 
seen as a direct product of localized holon and spinon like entities. 

Let us simply illustrate this by applying a sequence of various exchanges on 
the planar state \ip): 

- + - + - + - + 
+ - + - + - + - 

- + - + - + - + , (133) 
+ - + + - + - 

- + - + - + - + 
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where + and — denote up and down spins, respectively. Applying 
ToownT^ightTupTtft we arrive at jV'') 



- + - + - + - + 
+ - + - + - + - 

- -++- + - + , (134) 
+ + - + - + - 

- + - + - + - + 

a state which obviously differs from ToownTupT^ightTleftl'^) = 1"^)- Unlike the 
one dimensional case, damage is not kept under check. Note the extended 
domain wall neighboring the hole, enclosing a 2 x 2 region of spins of the 
incorrect registry. Note also that the hole is now surrounded by a pair of an- 
tiferromagnetically aligned spins along one axis and fcrromagnetically aligned 
spins along the other. A path closing on itself docs not lead to the fusion of the 
a "holon" and "spinon" like entities back into a simple hole. As the hole contin- 
ues to further explore both dimensions damage is continuously compounded. 
The state Tl^^^TRightTDownTRightTupTl^ft\i^) = 1-0) contains a string of eight 
spins of incorrect orientation surrounded by a domain wall whose perimeter 
is 16 lattice units long: 

- -++- + - + 
+ +-++- + - 

--- + - + -+ (135) 
+ + - + - + - 

- + - + - + - + 

As seen, the moving electron leaves a string of bad magnetic bonds in its 
wake. The energy penalty of such a string is linear in its extent. It is therefore 
expected that this (magnetic string) potential leads, in more than one dimen- 
sion, to a confining force amongst the spin and charge degrees of freedom. As 
this caricature for the single hole makes clear, the notion of localized "spinous" 
and "holons" is unlikely to hold water for the U ^ oo Hubbard model in more 
than one dimension. The well-defined SFL solutions for special models with 
nested Fermi-surface in two-dimensions should however be noted [94]. 

A certain form of spin-charge separation in two dimensions may be sought in 
the very special hole aggregates (or stripes) that have been detected in some of 
the cuprates [242] and the nickelates [58] . Here, holes arrange themselves along 
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lines which concurrently act as antiferromagnetic domain walls (i.e. behave like 
holons) in the background spin texture. Charge and spin literally separate and 
occupy different regions of space. In effect the two-dimensional material breaks 
up into one- dimensional lines with weak inter- connect ionj^. Related behavior 
is also found in numerical work [163] in the so called, t'-t-J model[^. An 
important question for such models is the extent to which the interconnections 
between stripes are "irrelevant" — i.e. the coupled-chains problem, which we 
briefly allude to in section 4.10. 



4-9 Recoil and the Orthogonality Catastrophe in Id and higher 

Here we show how the SFL behavior in one-dimension is intimately tied to 
the issue of orthogonality we discussed in sections 2.2 and 3.4. This line of 
thinking is emphasized by Anderson [20,24] who has also argued that this 
line of reasoning gives SFL behavior in two dimensions for arbitrary small 
interactions. 

We consider the effect of interactions through the explicit computation of our 
old friend from section 2.2, the quasiparticle weight 

Zk = {^^-'Vk\^''). (136) 

As we have seen, this indeed vanishes in all canonical one-dimensional models. 
Consider the model [49] of a Hamiltonian describing N fermions interacting 
with an injected particle via a delta function potential 

fl--^S^-^^VYj(^.--<i^- (137) 



This observation has led to a hne of thought which is of some interest in the 
context of the issues discussed here. If one focusses on the quantum mechanics of a 
single hne of holes by formulating it as a quantum mechanical lattice string model 
[85], the string traces out a two-dimensional world sheet in space-time. Quantum 
mechanical particles in one dimension, on the other hand, trace out world lines in 
space-time. It is claimed that one can recover most of the power law correlation 
functions of one-dimensional interacting fermions from the classical statistical me- 
chanics of fluctuating lines, and along these lines approach stripe formation as some 
form of spin-charge separation in two dimensions [285] . 

In this paper it was further observed that the kinetic motion of single holes may 
scramble the background spin environment in such a way that, on average, the 
holes may become surrounded by antiferromagnetically ordered spins on all sides 
(i.e. both along the horizontal and along vertical axis) — this is claimed to be a 
higher dimensional generalization of the holons encountered thus far. 
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The calculation for the small s-wave phase shifts for all single particle states 
{(f)i{xi)} is relatively straightforward. The quasiparticle weight Z reduces to 
an overlap integral between two (A^ + l)-particle Slater determinants, and one 
finds 

Z ~ Ar-2('5i^A)' (138) 

with 

5f^ -tan-^[UkF/2]. (139) 

In the thermodynamic limit Z = and no quasiparticles exist. As we see, 
the scattering phase shifts must conspire to give rise to anomalous behavior 
(exponents) for the electronic correlation function in such a way that they 
lead to a vanishing density of states at the Fermi level. We have already 
given explicit expressions for the anomalous exponent (s) under the presence 
of general scattering terms. 

As indicated in sections 2 and 3, an identically vanishing overlap integral 
between two (N+1) particle states could be a natural outcome of the emer- 
gence of additional Quantum numbers labeling orthogonal states. These states 
could correspond to different topological excitation sectors (e.g. solitons). Each 
quantum number corresponds to some conserved quantity in the system. In the 
one-dimensional electron gas this may be derived as we saw as a consequence 
of separate conservation law for left and right movers. 

An illustrative example of how singular Fermi liquid behavior due to orthog- 
onality of the wavefunction is robust in one dimension but easily destroyed in 
higher dimensions, is provided by the X-ray edge singularity problem, already 
discussed in section 3.5. As sketched in Fig. 19, we consider the transition of 
an electron from a deep core level to the conduction band through absorption 
of a photon. This problem is essentially the same as that of optical absorption 
in degenerate semiconductors, and from this point of view it is natural to an- 
alyze, following Nozieres [186], the effect of dispersion in the hole band, the 
analogue of the deep level state. For optical absorption in a semiconductor, the 
transition conserves momentum; hence in the absence of final state interac- 
tions, the threshold absorption is associated with the transition indicated with 
the arrowed line in Fig. 19, and absorption starts discontinuously above the 
threshold energy cud proveded the hole mass is infinite. For finite hole mass 
the threshold gets rounded on the scale of the dispersion of the hole band. 
But in one dimension, the edge singularity does survive because low-energy 
electron-hole excitations in one dimension have momenta only near and near 
2/ci?; see Fig. 10, electron-hole pairs can not carry away arbitrary momenta. 
This is seen in the following calculation [186]. 
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Assume a simple featureless final state potential V, and consider first the case 
without recoil. The relevant quantity to calculate is the transient propagator 
for the scatterer, 

G{t) = {0\de^'^'S\0) , (140) 

as the spectrum is the Fourier transform of G{t). In (140), the potential V is 
turned on at time and turned off at time t. In a linked cluster expansion, we 
may write G{t) = e*^*^*\ where C{t) is the contribution of a single closed loop. 
In lowest order perturbation theory, C{t) becomes 

t t 

C{t) = jdrj dr'V^gi^, r - r')g{r' - r) , (141) 



where g{0,r) is the free conduction electron propagator at the origin. For 
large times, one has (7(0, r) ^ —iN{0)/T, and when this is used in (141) one 
immediately find that for large times 

C{t) = V^N{Oyint , =^ G{t) ocl/r , (142) 

where 

n = V^N{Oy = S^n^ + 1 (143) 

is the phase shift exponent due to the orthogonahty effect, compare Eq. (56). 
A power law decay of G{t) ~ l/t" at long times corresponds to power law 
dependence ~ (a; — ujd)^~^ — {u — ujdY ^'^ just above the absorption edge. 

If we now take into account the recoil effect; the dispersion of the lower band 
implies that the hole in this band can hop from site to site. The propagator 
G{t) is then obtained as a sum over all trajectories -R(t) of the scattering hole 
which begin and end at i? = 0. For a given history, wc can extend the above 
analysis to lowest order by replacing the propagator 5'(0,t) by g{p,T), where 
p(r) = R{t) — R{t'). For positive time difference, we can then write 

g{p,T)^ ^ e-'^--*-". (144) 

k>kp 

For small hopping rates and large times, the integration over the modulus is 
dominated by the energy term, and this yields a term proportional to — l/r 
as in the recoilless case. The trajectory of the hole enters through the average 
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exp(— ik • p) over the Fermi sphere. A simple calculation yields 



exp(— ik ■ p) 



FS 



cos{kpp) d = I 
Joijipp) d = 2 



(145) 



In order to calculate the large time behavior of the Green's function, we finally 
have to average the square of this result over the distribution function of the 
trajectories p for large times. Using the large-p behavior of the expressions 
found above, one then easily finds that for large times 



|^exp(— ik ■ p) 



FS 




iiV2(0)/r2 d^l 
^\nptypiT)/{ptyj,{Ty) d = 2 ,(146) 



oc l/(p; 



typ 



d 



where ptyp is the typical distance the hole trajectory moves away from the 
origin. In one dimension, we see that still falls of as l/r^ and hence in 
analogy to (142) that G{t) has power law long time behavior: in the presence 
of recoil, an edge singularity persisits but the exponent n is now only half that 
in the absence of dispersion of the lower state, due to the averaging over p. 

1 

Since ptyp{T) grows diffusively as r2 for large r, the integrand in the expression 
(141) for C{t) converges faster than 1/t^, and hence C{t) converges to a finite 
limit for large times. The singular X-ray edge effect is washed out in two and 
three dimensions due to the recoil. 

If the above argument is extended to arbitrary dimension by analytically con- 
tinuing the angular average over the Fermi surface to continuous dimensions, 
one finds that the orthogonality and concomitant singular behavior is de- 
stroyed for any dimension d > 1. Nevertheless, the subdominant behavior of 
the integrals will contain noninteger powers of time, and this gives rise to sub- 
dominant nonanalytic terms in the spectrum for noninteger d. This behavior 
is completely in agreement with an analysis of the dimensional crossover from 
Luttinger liquid behavior to Fermi liquid behavior as a function of dimension 
[50]. 
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4-10 Coupled One- dimensional Chains 

The two-coupled chain problem has been thoroughly considered [88,93,217] 
[32,181] following earlier one-loop calculations [249] on a related model. The 
two-chain or ladder is especially interesting both theoretically and experi- 
mentally. In general interchain coupling is a relevant parameter, changing the 
behavior qualitatively. In the model of coupled Luttinger chains, the weight 
of the massless bosons characteristic of one dimension goes down with the 
number of chains. 



4.11 Experimental observations of one- dimensional Luttinger liquid behavior 

There has, of course, been a long-standing interest to observe the fascinating 
one-dimensional Luttinger liquid type SFL behavior experimentally, but the 
possibility of clear signatures has arisen only in the last few years. The clearest 
way to probe for Luttinger liquid behavior is to measure the tunneling into 
the one-dimensional system: Associated with the power law behavior (106) in 
one dimension, one has a power law behavior for the single electron tunneling 
amplitude into the wire. For fixed voltage, this leads to a differential conduc- 
tance dl/dV ~ l^", with the exponent a determined by the charge stiffness 
Kc, the geometry, and the band structure. Hence, from the measurement of 
the tunneling as a function of temperature or voltage, a and thus Kc can be 
extracted. Recent experiments on resonant tunneling [31] of small islands em- 
bedded in one-dimensional quantum wires in semi-conductors, grown with a 
so-called cleaved edge overgrowth method, do indeed yield a power law tem- 
perature behavior of the conductance which is consistent with Luttinger liquid 
behavior, but the value of the exponent is substantially different from the one 
expected theoretically. 

It has recently also been realized that nature has been kind enough to give 
us an almost ideal one-dimensional wire to study one- dimensional electron 
physics: the wave functions of carbon nanotubes turn out to be coherent over 
very large distances [73] . Though the circumference of the nanotubes is rather 
large, due to the band structure of the graphite-like structure the conduction 
in nanotubes can be described in terms of two gapless one-dimensional bands. 
Moreover, it was realized by Kane et al. [133] that due to the special geometry 
the backscattering in nanotubes is strongly suppressed, so that they are very 
good realizations of the Tomonaga- Luttinger model of section 4.2, with an 
interaction constant which is determined by the Coulomb energy on a 
cylinder. Their calculation based on this idea gives a value ~ 0.2. 

Fig. 32 shows recent nanotube data [197] which beautifully confirm the pre- 
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dictions by Kane et al. [133]: the differential conductance dl/dV is found to 
vary as V"" as a function of the bias voltage at low but fixed temperatures, or 
as T" as a function of temperature at fixed bias. For tunneling into the bulk 
of a carbon nanotube, the relevant density of states 

Nbuik{uj) ~ cu"''""^ with ai,uik = 7c (147) 

The fact that the exponent here is only half the value attained in the simple 
Luttinger Liquid — see Eq. (106) — is due to the fact that there are two Lut- 
tinger Liquid-like bands present in the carbon nanotube^^. Only one linear 
combination of the two associated charge modes attains a nontrivial stiffness 
^ 1 [133]. 

By contrast, at the tips of the nanotubes, surplus electronic charge can prop- 
agate in only one direction and, as a consequence, the tunneling is more re- 
stricted, 

Nup{u) u"'^'' with aup = {K,-l)/4 . (148) 



The electrostatic charging energy depends only on the symmetrized band mode, 
which in bosonic variables can be written as 6c,+ = {0 c,band=i + 0c,band=2 

)/V2. The 

essential reason that exponents can change depending on the number of bands is 
that the normalization factor l/\/2 in this bosonic variable enters in the exponent 
when the electron variables are written in terms of the bosonic modes, as discussed 
in section 4.2. A simple way to illustrate the halving of the exponents in the context 
of the various results we have discussed is by considering the difference between 
a the spinfull Luttinger liquid that we have discussed and the spinless Luttinger 
liquid (a system only having charge degrees of freedom): the spinfull model has 
a density of states exponent which is half of the spinless one. A calculation pro- 
ceeds along the following lines: for the spinfull case, the relevant Greens function 
is GR^=+i{x,t) = {e'^^cR{x,t)+^sR{x,t))/V2)^-i{^,Rm+^sRm)/V2)y Note the factor 

1 / \/2 in the exponent, coming from the projection onto the proper bosonic variables. 
This expection value can be written as GRa=+i{x,t) = (^e''^cR{x,t)/V2^-i'S'cR{o,o)/V2^ 
X ^e**-«(^'*)/^e-**-«(°'°)/^), which with the aid of the results of section 4.3 be- 
comes Gr„-_i_i (x, t) = 1 hrTToT^ — X 1 ^TTTjn — • From this result, 

one immediately obtains the density of states N{uj) = J dt G{x = 0, t) e*"^*. A sim- 
ple integration then yields N{uj) = [v'^^'^v-'^^'] J dt t'"" t~'^hc+is) ^iu>t ^ ^2(7^+7=). 
Since Kg renormalizes to 1 for repulsive interactions [see the remark just after (115)], 
one usually has 7<j = and so in this case the density of states exponent is simply 
. This is precisely (106). Now consider what one gets for the spinless Luttinger 
liquid. In this case neither the bosonic spin modes nor the projection factor \/y/2 
are present in the above expression for G. This results in a spatial decay with an 
exponent 27c instead of 7c, and hence an exponent 47c in the density of states. In 
other words, for 7^ = the density of states exponent in the spinful case is half of 
what it is in the spinless case. The same mechanism is at work in the nanotubes. 
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bla^ voltage (my) eV/I^T 

Fig. 32. Differential conduction dl/dV measured by Postma et al. [197] for carbon 
nanotubes. At low voltages or temperatures, Coulomb blockade effects dominate, 
but at higher temperatures or bias voltages, one probes the one-dimensional SFL 
behavior. Panel (a) shows the differential conductance as a function of bias voltage 
for various temperatures (note that these temperatures are relatively high, reflecting 
the fact that the electronic energy scales of the nanotubes is high). The effective 
exponent a for the large V behavior is 0.48; since these data are for tunneling 
between two nanotubes, a = 27c, so 7c ~ 0.24 and Kc ~ 0.27. The predicted value 
is Kc ~ 0.2 [133]. The data in panel (h) show the differential conductance as a 
function of temperature at fixed bias for two nanotubes which cross, as well as for 
a single nanotube with a bend. 

By Fermi's golden rule, the relevant exponents for (tip-tip) or (bulk- bulk) tun- 
neling arc = {Kc — l)/2 and ah-h ~ 27c respectively. Bulk-bulk tunneling 
is achieved by arranging the nanotubes according to the crossing geometry 
depicted in the inset of Fig. 32 above. By extracting the value of the charge 
stiffness K^. from each of the independent measurements of at-t-, O-b-h foi' the 
two different geometries, a single consistent stiffness Kc ~ 0.27 was found 
[197], in good agreement with theoretical prediction. This is a beautiful real- 
ization of Luttinger Liquid physics! 

Other, older, canonical realizations of Luttinger liquids include the Quan- 
tum Hall edge states. These represent a chiral spinless Luttinger hquid. Here, 
low lying energy states can only prevail at the edge of the sample, and, con- 
currently, disperse linearly about the Fermi energy. Edge states can attain 
macroscopic linear (perimeter) extent, and the tunneling experiments between 
such states [170,110] have observed several features predicted theoretically 
[271,132]. We refer the interested reader to the review articles by Schulz et al. 
[218] and Fisher and Glazman [96] for further details. 
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5 Singuleir Fermi-liquid behavior due to gauge fields 



5. 1 SFL behavior due to coupling to the electromagnetic field 



Almost thirty years ago Holstein, Pincus and Norton [120] (see also Reizer 
[207]) showed that the coupling of electrons to the electromagnetic fields gives 
rise to SFL behavior. Since the typical temperatures where the effects become 
important are of the order of 10""^^ K, the effect is not important in practice. 
But the theory is of considerable general interest. 

If we work in the Coulomb gauge in which V • A = for the electromagnetic 
A field, then the transverse propagator D^^ in free space is given by 

Dli,Ko.) = (A,A,)(^^.) = J^^-^^, (149) 

where e is an infinitesimal positive number. The interaction of the electrons 
with the electromagnetic field is described by the coupling term 

yA + pA\ (150) 



where j is the electron current operator and p the density operator. 

Quite generally in the Coulomb gauge, one finds from perturbation theory, 
or phenomenologically from the Maxwell equations, that the electromagnetic 
propagator in a metal can be written as 

D-^ = - M [5,^ - kkj)'' . (151) 



The perturbative diagrammatic expansion of M is indicated in Fig. 33. The 
first term leads to a term proportional to the density while the second term 
is the first correction due to particle-hole excitations; in formulas, these terms 
yield 

M(k,-) ^^"fnl ^ / [p'-(P-k)-][/.-A:/2- W] 
my mJ {2ny uj - {ep-k/2 - ep+fc/2) 



Here n is the electron density, and a — 1/137 is the fine structure constant. 
For Lo/k ^ the two terms combined yield 



M{k,uj) 



Snan 



m 



kvp 



+ 2i 



UJ 

kvp 



(153) 
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Fig. 33. (a) The diagrams for M; (b) The diagram for the self-energy from the in 
the discussion of the SFL effect arising from the couphng of the electrons to the 
electromagnetic gauge fields; (c) Anharmonic interaction of fluctuations such as 
(c) are non-singular in the SFL problem of coupling of the electrons in metals to 
electromagnetic fields. 

At small frequencies, this yields for the propagator D 



Dij{\<i,uj) 



nkvF 



(154) 



which corresponds to an overdamped mode with dispersion uj k^. 

Before discussing how such a dispersion gives SFL behavior in three dimen- 
sions, it is instructive to point out that although (154) was obtained perturba- 
tively. Maxwell's equations ensure that the field propagator must generally be 
of this form at low frequencies and momenta. Indeed, for a metal we can write 
the current j as j = a{k,uj)'E; if we combine this with the Maxwell equation 
V X H = j + d'E/dt we easily find that the general form of the propagator is 

Ai(k,a;) = . (155) 

Here /i is the diamagnetic permeability. For pure metals, the low-frequency 
limit is determined by the anomalous skin effect [5] and cT(k, 0) ~ k~^ . Accord- 
ing to the expression (155), this 1/k behavior then implies that the dispersion 
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at small frequencies generally goes as ~ fc^ for a pure metal. For dirty met- 
als, the fact that a{k, 0) approaches a finite limit o"o at small wavenumbers, so 
for dirty metals there is according to (155) a crossover to a behavior u; ~ fc^ 
at small wavenumbers. 

Gauge invariance of the theory requires that the photon cannot acquire a 
mass (a finite energy in the limit /c — > 0) in the interaction process with 
the electrons, and hence the form of Eq. (154) remains unchanged. Thus the 
anharmonic corrections to the photon propagator due to processes such as 
shown in Fig. 33(c) do not change the form of Eq. (154). The self-energy 
of the electrons due to photon exchange. Fig. 33(b), may now be calculated 
with confidence given the small coupling constant in lowest order. The leading 
contribution in c? = 3 is 

S(fcF, to") ~ 0(0; lnc<j + icusgncj) (156) 



The momentum dependence, on the other hand, is nonsingular as a function 
of k — kp. 

The nonanalytic behavior of the self-energy as a function of frequency implies 
that the resistivity of a pure metal in d = 3 is proportional to T^/^ at low 
temperature^^. 

The simplicity and strength of the above example lies in the fact that the 
theory has no uncontrolled approximations, and the gauge-invariance of the 
photon field dictates the low-energy low-momentum behavior of the photon 
propagator Djj(k, cj). Moreover, vertex corrections are not important because 
the Migdal theorem is valid [195] when the frequency of fluctuations is very 
small compared to their momenta, as in a; ~ g^. 

However, the SFL behavior as a result of the coupling to the electromagnetic 
field is not relevant in practice. This can most easily be argued as follows. For 
a Fermi gas, the entropy per particle is 

S_ _ vr^mfclT 



This follows from the fact that in the quasi-elastic approximation, the transport 
relaxation rate rj~^ is related to the single-particle relaxation rate t~^{6) due to 
scattering through an angle 9 near the Fermi-surface 



'tr 



J dn{l-cos{e))T-\9). (157) 



For small T, w ~ A;'^ ~ T, and hence the characterestic angle of scattering, 9 ~ 
{k/kp) ~ {T/EpY/^ is small. Upon expanding (1 — cos^) ~ ^^/2, one finds that 
the effective transport scattering rate goes as / d?k k'^f{k^/T) ~ T^/'^. 
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while for the entropy for the electrons interacting with the electromagnetic 
field one finds from the above results [120] 



S 



27r2 



3 ck} 



■f \T ) 



(159) 



N 



with 



cep 



(160) 



Upon comparing these two results, one concludes that the SFL effects start to 
become important for temperatures 



Since the numerical factor in the exponent is typically of order 10^//"^ for 
ordinary metals, the temperature range one finds from this is of order 10~^^ 
Kelvin for values of of order unity, according to this estimate. Note however 
that for pure ferromagnetic metals fi can be as large as 10^. Possibly, in some 
ferromagnets, the effects can become real. 

5.2 Generalized gauge theories 

The example of coupling to the electromagnetic field identifies one possible 
theoretical route to SFL behavior, but as we have seen the crossover tempera- 
tures that one estimates for this scenario are too small for observable physical 
properties. The smallness of the estimated crossover temperature is essentially 
due to the fact that the couphng to the electromagnetic field is determined by 
the product avp/c, where a — 1/137 and typically vp/c = C(10~^). Motivated 
by this observation, many researchers have been led to explore the possibility 
of getting SFL effects from coupling to different, more general gauge bosons 
which might be generated dynamically in strongly correlated fermions. For 
a recent review with references to the literature, see [150,179,180]. The hope 
is that if one could consistently find such a theory in which the small factor 
avp/c arising in the electromagnetic theory is replaced by a term of order 
unity, realistic crossover temperatures might arise. Much of the motivation 
in this direction comes from Anderson's proposal [24] of spin-charge separa- 
tion and resonating valence bonds in the high-temperature superconductivity 
problem which we discuss in section 7. 

The essence of approaches along these lines is most easily illustrated by consid- 
ering electrons on a lattice in the case in which strong on-site (Hubbard-type) 




(161) 
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repulsions forbid two electrons to occupy the same site. Then each site is ei- 
ther occupied by an electron with an up or down spin, or by a hole. If we 
introduce fictitious fermionic creation and annihilation operators and / for 
the electrons and fictitious bosonic hole creation and annihilation operators 
and b for the holes, we can express the constraint that there can only be one 
electron or one hole on each lattice site by 

jZflfio^ + blbi^l for each i. (162) 

a=l 

With this convention, the real electron field ipia can be written as a product 
of these fermion and boson operators 

Ipia = fiabi . (163) 

This expresses the fact that given the constraint (162), a real fermion annihi- 
lation at a site creates a hole. This way of writing the electron field may be 
motivated by the physics of the one-dimensional Hubbard model: there a local 
excitation may indeed be expressed in terms of a charged spinless holon and 
an uncharged spinon. In general, in a transformation to boson and fermion 
operators as in (163) there is some freedom with which operator we associate 
the charge and with which one the spin. 

Whenever we split a single electron operator up into two, as in (163), then 
there is a gauge invariance, as the product is unchanged by the transformation 

fiait) e*'(*V«a(^) , b,{t) e'^^%{t) . (164) 

We can promote this invariance to a dynamical gauge symmetry by intro- 
ducing a gauge field a and writing the Hamiltonian in the continuum limit 

as 

j dvb^ [-'^;^^\+{<t>- Mb . (165) 

Here is a Lagrange multiplier field which is introduced to implement the 
constraint (162). Note that the a field enters in much the same way in the 
Hamiltonian as the electromagnetic field usually does — indeed, a change 
in the individual /'s and 6's by a space dependent phase factor as in (164) 
can be reabsorbed into a change of a. Note also the presence of the chemical 
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potentials / and fif, to enforce that for a deviation of the (average) density n 
from one per site, the density of holes is 1 — n and of fermions is n. 

The next step in the theory is to find the fluctuation propagator for the a field, 
as a function of /// and //{,. For finite /// and negative (bosons uncondensed) , 
the fluctuation propagator is similar to that of the previous section but with 
{avp/c) replaced by a term of (9(1). Spinon and ho Ion self-energies can now be 
calculated and composition laws [126] are derived to relate physical correlation 
functions to correlation functions of spinous and holons. 

Unfortunately, this very attractive route has turned out to be less viable than 
had been hoped. It is not clear whether the difficulties are purely technical; 
they are certainly formidable. The essential reason is that while photons have 
no mass and are not conserved, and hence can not Bose condense, in a gauge 
theory obtained by introducing additional bosons, the bosons generally can 
and will Bose condense because they do have a chemical potential. Bose con- 
densation leads to a mass term in the propagator for the gauge fields. The 
singularities in the fermions due to the gauge fields then disappear. It is the 
analogue of the fact that superconductivity leads to the Meissner effect — 
there the emergence of the superconducting field breaks gauge invariance and 
leads to the expulsion of the magnetic field from the superconductor. The lat- 
ter effect can also be thought of as being due to the generation of a mass term 
for the gauge fields. 

Several variants of these ideas have been proposed with and without attempts 
to suppress the unphysical condensation through fluctuations [176]. The trou- 
ble is that such fluctuations tend to bind the spinous and holons and the happy 
situation in one dimension where they exist independently — being protected 
by (extra) conservation laws, see section 4.9 — is hard to realize. As usual, it 
appears that introduction of new quantum numbers requires new symmetries. 

Interesting variants using the idea of spin-charge separation have recently 
appeared [220,176]. 



6 Quantum Critical Points in fermionic systems 

As mentioned in the introduction, quantum critical behavior is associated 
with the existence of a T = phase transition; of course, in practice one can 
only study experimentally the behavior at nonzero temperatures, but in this 
sense the situation is no different from ordinary critical phenomena: one never 
accesses the critical point itself, but observes the critical scaling of various 
experimentally accessible quantities in its neighborhood. 
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In practice, the most common situation in which one observes a quantum crit- 
ical point is the one sketched in Fig. 1, in which there is a low temperature 
ordered state — a ferromagnetic state, antiferromagnetic state or charge den- 
sity wave ordered state, for instance — whose transition temperature to the 
disordered state or some other ordered state goes to zero upon varying some 
parameter. In this case, the quantum critical point is then also the end point 
of a T = ordered state. However, sometimes the "ordered" state really only 
exists at T = 0; for example in metal-insulator transitions and in quantum 
Hall effect transitions [232,198]. A well-known example of this case in spin 
models is in two-dimensional antiferromagnetic quantum Heisenberg models 
with "quantumness' as a parameter g [55], which do not order at any finite 
temperature, but which show genuine ordered phases at T = below some 
value g < Qc- 

Although the question concerning the origin of the behavior of high tempera- 
ture superconductors is not settled yet, there are strong indications, discussed 
in the next section, that much of their behavior is governed by the proximity 
to a quantum critical point. 

One of the first formulations of what we now refer to as quantum critical be- 
havior was due to Moriya [177,178] and Ramakrishnan [174], who did an RPA 
calculation for a model of itinerant fermions with a Stoner-type instability to 
a ferromagnetic state. In modern language, their approach amounts to a 1/iV 
expansion. Various other important contributions were made [37,77]. The now 
standard more modern formulation, which we will follow, is due to Hertz [117]. 
A nice introduction can be found in the article by Sondhi et al. [232], and for 
a detailed expose we refer to the book by Sachdev [210]. See also the review 
by Continento [66]. 



6.1 Quantum critical points in ferromagnets, antiferromagnets, and charge 
density waves 

A clear example of quantum critical behavior, and actually one for which 
one can compare with theoretical predictions, is summarized in Figs. 34-37. 
The figures show various data from [239] on the magnetic compound MnSi 
[191,192,152,161,239]. Fig. 34 shows that for low pressures and temperatures, 
this compound exhibits a magnetic phase whose transition temperature 
vanishes as the pressure is increased up to Pc — 14.8 khar. This value of the 
pressure then identifies the quantum critical point. Fig. 35 shows that when 
the same data are plotted as T^'^ versus pressure, the data fall nicely on a 
straight line except close to the critical pressure. This observed scaling of T^^"^ 
with p— Pc away from pc is in accord with the behavior predicted by the theory 
discussed below. Actually the transition is weakly first order near so very 
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Fig. 34. Magnetic phase diagram as a function of pressure of MnSi [191,192,239]. 
Prom [239]. 
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Fig. 35. Power law dependence of the Curie-temperature as a function of pressure 
for MnSi. Prom [239]. 



detailed verification of the theory is not possible. 



Fig. 36 shows data for the temperature dependence of p — po near pc, where po 
is the residual low temperature resistivity. In the presence of a field of 3 Tesla, 
one observes the usual p — po ~ Fermi-liquid scaling, but at zero field the 
results are consistent withp — po ~ T^/^ behavior predicted by the theory. But 
if we write the low temperature resistivity behavior as 

p = po + AT^ (166) 
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Fig. 37. Evolution of po and A under pressure, when the temperature dependence 
of the resistivity is fitted to {p = po + AT"^). Prom [239]. 



then both the residual resistivity po and the amphtude A are found to show 

a sharp peak at pc as a function of pressure — see Fig. 37. This behavior 
is not understood nor is the fact that the exponent 6 does not appear to 
regain the Fermi-hquid value of 2 for significant values of p > Pc (at H = 0) 
and in a temperature regime where the theory would put the material in the 
quantum-disordered Fermi-liquid regime. 
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6.2 Quantum critical scaling 



Before discussing other experimental examples of quantum critical behavior it 
is expedient to summarize some of the essential quantum critical scaling ideas. 

As is well known, at a finite temperature transition, the critical behavior is 
classical and we can use classical statistical mechanics to calculate the correla- 
tion functions. This is so because due to critical slowing down the characteristic 
time scale r diverges with the correlation length, 

T ~ r ■ (167) 

Near a critical point the correlation length ^ diverges as 

^~|r-7;r''- (168) 

The combination of these two results shows that critical slowing down implies 
that near any finite temperature critical point the characteristic frequency 
scale LUc goes to zero as 

a;e ~ |r - Tcr' ■ (169) 

Therefore near any phase transition cUc <^Tc, and as a result the phase transi- 
tion is governed by classical statistical physics: the Matsubara frequencies are 
closely spaced relative to the temperature, the thermal occupation of bosonic 
modes is large and hence classical, etc. 

In classical statistical mechanics, the dynamics is a slave to the statics; usually, 
the dynamical behavior is adequately described by time-dependent Landau- 
Ginzburg type of equations or Langevin equations which are obtained by 
building in the appropriate conservation laws and equilibrium scaling behavior 
[112]. At a quantum critical point, the dynamics must be determined a priori 
from the quantum-mechanical equations of motion. 

The general scaling behavior near a T = quantum critical point can however 
be discussed with the formalism of dynamical scaling [112,156], just as near 
classical critical points. Consider for example the susceptibihty for the case of 
MnSi that we considered above. The scaling ansatz for the singular part of 
the susceptibility 

x{k,u;,p)^{MM){k,u;,p) (170) 
implies that near the critical point where the correlation length and time scale 
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diverge, the zero-temperature susceptibility x is a universal scaling function 
of the scaled momentum and frequency, 



Xik,u;,p) = r'''r{k^,cu^,) , (171) 

where now 

e~b-Pcr% 6 = r = r. (172) 

This is just like the classical scaling with T — Tc replaced by p — pc- The reason 
for writing instead of r is that in quantum statistical calculations, the "time- 
wise" direction becomes like an additional dimension, so that plays the role 
of a correlation "length" in this direction. However, the time-direction has both 
a long-time cutoff given by l/ksT and a short-time cutoff given by the high 
energy cutoff in the problem — exchange energy or Fermi-energy, whichever 
is smaller in the ferromagnetic problem. The short-time cutoff has its analog 
in the spatial scale. The long-time cutoff, which determines the crossover from 
classical to quantum behavior, plays a crucial role in the properties discussed 
below. The crucial point is that when 2; 7^ 1, there is anisotropic scaling 
between the space and time-wise direction, and as we shall discuss below, this 
implies that as far as the critical behavior is concerned, the effective dimension 
of the problem is d + z, not d + 1. 

The exponent dM in (171) reflects that a correlation function like x has some 
physical dimension which often is inevitably related to the spatial dimension. 
The dependence of critical properties on spatial dimensions must be express- 
ible purely in terms of the divergent correlation length ^. Often d^i is fixed by 
dimensional considerations (in the language of field theory, it is then given by 
the "engineering dimension" of the field), but this may not be true in general. 
It must be so, however, if x is a correlation function of a conserved quant ityp^. 

Let us now address the finite temperature scaling, taking again the case of 
MnSi as an example. The various regimes in the T-p diagram discussed below 
are indicated in Fig. 38. To distinguish these regimes it is necessary to define 
an additional quantity, the thermal length 




(173) 



Lt corresponds dimensionally to a timescale. It marks the crossover between 
phenomena at long time scales which can be treated essentially classically from 

E.g., if we consider the free energy per unit volume at a classical transition, the 
energy scale is set by ksT, and cIm = d; likewise, if we consider the surface tension 
of an interface, whose physical dimension is energy per unit surface area, dM = d—1. 
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Fig. 38. Generic phase diagram and cross-overs for quantum critical points with the 
various regimes indicated. ^ is the correlation length at T = 0; Lt is the "thermal 
correlation length". 

those on a shorter scale which are inherently quantum mechanical. Whenever 
< Lt, the correlation length and time are finite and quantum mechanics 
begins to dominate. This is the regime on the right in the figure. Fermi-liquid 
behavior is expected in this regime. However, if one approaches the critical 
point (T = 0, p = pc) from above along the vertical line, then Lt diverges 
but so does ^j. Moreover, diverges faster than Lt since z is usually larger 
> 1. This means that the characteristic fiuctuation energy and temperature 
are similar. So the behavior is quasiclassical throughout each correlated region 
down to zero temperature (In a path integral formulation [232,210] one consid- 
ers the model on a infinite strip whose width is finite in the timewise direction 
and equal to Lt- Hence for > Lt the model is fully correlated across the 
strip in this direction). This regime is therefore characterized by anomalous 
T dependence in the physical quantities up to some ultra-violet cutoff. It is 
important to stress that this so-called "quantum critical scaling behavior" 
is expected in the observable properties in the whole region above pc- So SFL 
can be observed over a whole range of parameter p and temperatures between 
the left and right cross-over lines. 

If one approaches the line of phase transitions to the ordered phase, which is 
marked with a solid line in the figure, one has a region with SFL properties 
dominated by classical fiuctuations close to the transition. Millis [171] has 



The term is somewhat problematic; it refers to the quasiclassical fluctuation 
regime around a quantum critical point. 
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corrected Hertz's results [117] on this point, and has found that the critical 
temperature of the phase transition scales as ~ \p — PcY/^'^'^^~'^\ Estimates 
of the classical critical region are also given. Results along similar lines may 
also be found in [65,66]. 

If we include both the temperature and the parameter p, the scaling ansatz 
for the imaginary part of x becomes 

x"{k,uj,p,T) ~ r'^^Ti(A;e,^r,^) • (174) 

This can be rewritten in other forms depending on which experiment is being 
analyzed. For example, the above form is especially suitable for analysis as a 
function of (p — Pc)- For analyzing data as a function of temperature, we may 
instead rewrite 

x"{k,u,p,T) = L^'-'^'T, (kL'J^uLT, ^] , (175) 



and for analyzing data as a function of frequency 

X%k,u,p,T)=T~'^^'/^T,(^-±^,^,±-^ . (176) 

Moreover, the scaling of the free energy JF can be obtained from the argument 
that it is of the order of the thermal energy /c^T per correlated volume 
Moreover, since Lt acts as a finite cutoff for C,t in the timewise direction, we 
then get the scaling 

jr _ TC'^ ~ r(6)-^/" ~ Ti+°'/" . (177) 

By differentiating twice, this also gives immediately the specific heat behavior 
at low temperatures. In writing the above scaling forms, we have assumed that 
no "dangerously irrelevant variables" exist, as these could change u/T scaling 
to uj/T^ scalingf^. 

"Dangerously irrelevant variables" are irrelevant variables which come in as pref- 
actors of scaling behavior of quantities like the Free Energy [95] . Within the Renor- 
malization Group scenario, the hyperscaling relation di' = 2— a is violated above the 
upper critical dimension because of the presence of dangerously irrelevant variables. 
Presumably, dangerously irrelevant variables are more important than usually at 
QCP's, since the effective fluctuation dimension is above the upper critical dimen- 
sion for d = 3 and z > 1. Some examples are discussed in [210]. 
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In order to get the critical exponents and the cross-over scales, one has to turn 
to a microscopic theory. The theory for this particular case of the quantum 
critical point in MnSi is essentially a Random Phase Approximation and 
proceeds along the following lines: (i) One starts with a model of interacting 
fermions; (ii) An ordering field M{k,uj) is introduced; (in) One assumes the 
fermions can be eliminated near the critical point to get a free energy in terms 
of M of the form 



Note that this is essentially an extension of the usual Landau-Ginzburg- Wilson 
free energy to the frequency domain. Indeed, from here on one can follow the 
usual analysis of critical phenomena, treating the frequency ut on an equal 
footing with the momentum k. 

The important result of such an analysis is that the effective dimension as far 
as the critical behavior is d + z, not d + 1 as one might naively expect. Since 
z > 1 in all known examples, the fact that the effective dimension is bigger 
than d + 1 reflects the fact that the correlation "length" in the timewise di- 
rection grows as i.e. faster than the spatial correlation lengthp^. Moreover, 
the fact that z > 1 implies that the effective dimension of a c? = 3 dimen- 
sional problem is always larger or equal than four. Since the upper critical 
dimension above which mean field behavior is observed equals four for most 
critical phenomena, one thus arrives immediately at the important conclusion 
that most quantum critical points in three dimensions should exhibit classical 
fluctuations with mean field scaling exponents! It also implies that the critical 
behavior can typically be seen over a large parameter or temperature range — 
the question of the width of the critical region, which normally is determined 
by the Ginzburg criterion, does not arise. On the other hand, questions con- 
cerning the existence of dangerously irrelevant variables, due for example to 
the scaling of the parameters V in Eq. (178), do arise. 

In order to judge the validity and generality of these results, it is important 
to keep in mind that they are derived assuming that the coefficients of the 
M^, terms are analytic functions of k, oj and the pressure p, etc. This is 
completely in line with the usual assumption of analyticity of the bare coupling 
parameters in a renormalization group approach. This assumption may well 

This has important consequences for a scaling analysis of numerical data, aimed 
at determining the critical behavior. For it implies that the finite size scaling has 
to be done anisotropically, with the anisotropy depending on the exponent z which 
itself is one of the exponents to be determined from the analysis 




+ J d{cu} J d{k} VM{ki,u;i)M{k2,uj2)Mk3,uj3)M{ki,uj4) 

X 5{uJi +UJ2 + UJS + uj^)5{ki + k2 + ks + k^) + --- (178) 
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be violated — in fact none of the impurity models discussed earlier can be 
treated along these lines: the fermions can not be integrated there, and if one 
attempts to apply the above procedure, one finds singular contributions to 
the bare coupling parameters. Later on we shall discuss a three-dimensional 
experimental example where this assumption appears to be invalid. Secondly, 
it is inherently an expansion about the nonmagnetic state, which can not 
apply in the ordered phase: In the ordered phase with nonzero magnetization, 
M 7^ 0, there is a gap for some momenta in the fermionic spectra. This gap 
can not be removed perturbatively. 

In a ferromagnet, the ground state susceptibility on the disordered side is given 

by 

X-'(A;,a;)= [(p-pe) + /c' + £^] . (179) 



In the first two terms we recognize the usual mean field type behavior with a 
correlation length that diverges as ^ ~ {p ~ PcY^'^y hence the critical exponent 
u = 1/2. The last term, which describes Landau damping of the spin wave 
modes, is very special here as it arises from fiuctuations of magnetization, a 
quantity which is conserved, i.e., commutes with the Hamiltonian. Therefore 
the characteristic damping rate must approach zero as /c — > 0. Since at the 
critical point (179) leads to a damping uj ~ k^, the critical exponent z = 3. 
According to the theory described above, the critical behavior at the quantum 
critical point (T = 0, p=Pc) is therefore of the mean-field type for any physical 
dimension d > 1, since the effective dimension d-|-2;>4for(i>l. 

The scattering of electrons off the long wavelength spin waves is dominated 
by small angle scattering, and it is easy to calculate the resulting dominant 
behavior of the self-energy of the electrons. Near the critical point the behavior 
of X is very similar to the electromagnetic problem we discussed in section 5. 
Analogously, one also finds SFL behavior here: in c? = 3, 'L{kp,uj) ~ a;lncj -|- 
i\uj\ while in = 2 one obtains Tjikp^uj) ~ cj^/^ + i|cjp/^. Furthermore, for the 
resistivity one finds in three dimensions p ~ T^/^ — this is consistent with the 
behavior found in MnSi in the absence of a field, see Fig. 36. Moreover, as we 
mentioned earlier, according to the theory, near the critical point should 
vanish as \p — p^|^/('^+2-2)- -v^rith d = z = ?> this yields Tc ~ \p — Pc\^^^- As we 
saw in Fig. 35, this is the scaling observed over a large range of pressures, 
except very near pc- ZrZn2 [HI] is an example in which the Ferromagnetic 
transition is shifted to T = under pressure continuously. The properties are 
again consistent with the simple theory outlined. There is however trouble on 
the horizon [153]. The asymptotic temperature dependence for p > is not 
proportional to T^, as expected. We will return to this point in sections 6.4. 

For antiferromagnets or charge density waves the critical exponents are dif- 
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ferent. In these cases, the order parameter is not conserved, and the inverse 
susceptibihty in these cases is of the form 



lOJ 



{p-Pc) + {k-koY + — 



(180) 



where is the wavenumber of the antiferromagnetic or charge density wave 
order. From this expression we immediately read off the mean field exponents 
z = 2 and z/ = 1/2. Since the effective dimension d + z is above the upper 
critical dimension for o? = 3, the mean field behavior is robust in three dimen- 
sions. In d = 2, on the other hand, the effective dimension d + z = 4 is equal 
to the upper critical dimension, and hence one expects logarithmic corrections 
to the mean field behavior. Indeed, in two dimensions one finds for the self- 
energy [119] S(k^,a;) ~ cij InciJ -|- i|t(;| (only)for those from which spanning 
vectors to other regions of the Fermi-surface separated by ko can be found; the 
resistivity goes as p{T) ~ In T in this case. However, if several bands cross 
the Fermi-surface, as often happens in heavy fcrmions, Umklapp-type scat- 
tering may enforce the same temperature dependence in the resistivity as in 
the single particle self-energy, except at some very low crossover temperature. 
The physical reason for this dependence despite the fact that the soft modes 
are at large momentums (and therefore vertex corrections do not change the 
temperature dependence of transport relaxation rates) is that the set of k^? 
usually covers a small portion of the Fermi-surface. 

Fig. 39 shows the phase diagram of two compounds that order antiferromag- 
netically at low temperatures. The first one, Celn^ is a three-dimensional 
antiferromagnct. A superconducting phase intervenes at very low tempera- 
tures (note the different scale on which the transition to the superconducting 
phase is drawn), covering the region around the quantum critical point at 
a pressure of about 26 kbar. At this pressure, the normal state resistivity 
is found to vary as p{T) ~ j'i-6±o.2 ^jiicj^ jg consistent with the theoretical 
prediction that at a quantum critical point dominated by antiferromagnetic 
fluctuations the resistivity should go as p ~ T^'^. The right panel in Fig. 39 
shows the phase diagram and resistivity data of the three dimensional anti- 
ferromagnct CePd2Si2, the data in this case are best fitted by p ~ T^-'^; this 
is consistent with the theoretical prediction p ~ T^-^^ which results if one has 
a (/c — dispersion around the AFM vector ko in one direction and the 
usual {k — koY dispersion in the other two. No independent evidence for such 
dispersion is however available yet. In both of these cases, part of the region 
of superconductivity, in a region bounded by a line emanating from the QCP 
and going on to the antiferromagnetic to the normal state line, is expected to 
be antiferromagnetic as well. 

These two compounds are also of interest because the phase diagram bears a 
resemblance to the phase diagram of the high- T^. copper-oxide based super- 
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Fig. 39. Left panel: Temperature-pressure phase diagram of high-purity sin- 
gle-crystal Celn^. A sharp drop in the resistivity consistent with the onset of 
superconductivity below Tc is observed in a narrow window near pc the pressure 
at which the Ncl temperature Tat tends to absolute zero. Upper inset: this tran- 
sition is complete even below Pc itself. Lower inset: just above pc, where there is 
no Nel transition, a plot of the temperature dependence of d(lnA/9 )/d(ln T) is 
best able to demonstrate that the normal state resistivity varies as T^-^^^-^ below 
about 3 K. Ap is the difference between the normal state resistivity and its residual 
value (which is calculated by extrapolating the normal-state resistivity to absolute 
zero). For clarity, the values of Tc have been scaled by a factor of ten. Right panel: 
Temperature-pressure phase diagram of high-purity single-crystal CePd2Si2- Su- 
perconductivity appears below Tc in a narrow window where the Nel temperature 
Tjv tends to absolute zero. Inset shows the normal state a-axis resistivity above the 
superconducting transition varies as T^-'^^'^-^ over nearly two decades in tempera- 
ture. The upper critical field Bc2 at the maximum value of Tc varies near Tc at a 
rate of approximately —6T/K. For clarity, the values of Tc have been scaled by a 
factor of three, and the origin of the inset has been set at 5 K below absolute zero. 
Both plots are from Mathur et al. [161]. 



conductors in which the conduction electron density is the parameter varied 
— see. Fig. 49. Unlike the heavy fermion compounds where the ordered phase 
is antiferromagnctic, the order in copper-oxides near the QCP is not AFM. Its 
nature is in fact unknown. In the heavy fermion compounds superconductivity 
promoted by antiferromagnctic fluctuations is expected to be of the (i-wave 
variety [175] as it is in the high- Tc copper-oxide compounds. 
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6.3 Experimental Examples of SFL due to Quantum- Criticality: Open Theo- 
retical Problems 

We have discussed the observed quantum critical behavior in some system 
which is largely consistent with the simple RPA-like theoretical predictions. 
There are however quite a few experimentally observed signatures of singu- 
larities near QCPs, especially in heavy fermion compounds, which are not 
understood theoretically by the simple RPA theory of the previous subsec- 
tion. In this section we present some prominent examples of these. 

The experimental observations fall in two general categories, in both of which 
the low-temperature resistivity does not obey the power laws expected of 
Fermi-liquids: Compounds in which resistivity decreases from its limit at T = 
and those in which it increases. In both cases the C^/T is singular for 
T ^ 0. It is reasonable to associate the former with the behavior due to 
impurities and the latter with the QCP properties of the pure system. But 
as we discussed in section 3.9 the QCP due to impurities requires tuning to 
special symmetries unlikely to be realized experimentally. As we will discuss, 
the effect of impurities without any special symmetries but coupling to the 
order parameter is expected to be quite different near the QCP of the pure 
system compared to far from it. Under some circumstances, it is expected to 
be singular and may dominate the observations. 

We start with experiments in the second category. Figs. 40-42 shows several 
data sets for the heavy fermion compound CeCuQ^^Aux for various amounts of 
gold. This compound shows a low-temperature paramagnetic phase for a; < 0.1 
and a low temperature antiferromagnetic phase for x > 0.1. The specific heat 
data of Fig. 40 and the resistivity data of Fig. 42 show that while without Au, 
i.e. in the paramagnetic regime, the behavior is that of a Heavy Fermi-liquid 
metal, the alloy near the quantum critical composition CeCu^,^AuQ,i exhibits 
a specific heat with an anomalous 

C^^TXnT. (181) 

At the same composition, the resistivity shows a linear temperature depen- 
dence, and the susceptibility data in Fig. 41 show a v^T cusp. The anomalous 
behavior is replaced by Fermi-liquid properties both by a magnetic field and 
by increasing the substitution of copper by gold. The compound YbRh2S2 
seems to have similar properties [243]. Related properties have also been 
found in U2Pt2ln [86,87] and in UPt^.^Pd^ [74], UBei3,CeCu2Si2, CeA^isGcs 
[235]. The SFL properties observed at the Mott-Insulator-Metal Transition in 
BaVSs [97] are also of related interest. 

None of the quantum critical properties of the CeCuAu compounds is con- 
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Fig. 40. The specific heat C /T of CeCu^-^cAu X versus log T. From Lohneysen et al. 
[266-268]. 
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Fig. 41. Susceptibility data for CeCuAu. From Lohneysen et al. [269]. 

sistent with any of the models that we have discussed. Information on the 
magnetic fluctuation spectra for CeCu5,gAuo,i is available through neutron 
scattering experiments [215,236]. The data shown in Fig. 43 show rod-hke 
peaks, indicating that the spin fluctuations are almost two-dimensional at 
this composition. The neutron scattering data can be fitted by an expression 
for the spin susceptibility of the form 



X-\k,u;)^C[f{5k) + {-iu + aTr] 



(182) 



with a function / which is consistent with an effectively two-dimensional scat- 
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Fig. 42. Resistivity data for CeCuAu. From Lohneysen et al. [269]. 
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Fig. 43. Neutron scattering data for CeCurj,%AuQ i, a compound which is close to 
a QCP. The figure shows g-scans along three different crystallographic directions, 
from top to bottom in the a, h and c directions for Tiu = O.lmeV. The figures shows 
that there is only a weak g-dependence along the rods (g±), while transverse scans 
shows well defined peaks with nearly the same line width. Prom Stockert et al. 
[236]. 
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tering, 



f{dk) = b{dk±y + c{6kii)^ . (183) 

6k\\ and 6k± are the deviations from the AFM Bragg vector parallel and per- 
pendicular to the c-axis in these tetragonal crystals. Further 

a^l, q; = 0.74 ±0.1. (184) 

At present, there is no room for this anomalous exponent a within known the- 
oretical frameworks. However, if one accepts this particular form of x ^ giving 
an adequate fit, the observed specific heat follows: at the critical composition, 
we expect the temperature scaling relation 

-p ^ rjn^-id-l)^-l ^ rprp{d-l)a/2rpa/4: ^ rpl+{d-l/2)a/2 ^ ^^gg^ 

which immediately gives Cy ~ (PJ- /dT"^ ~ T for a = 4/5 with a plausible 
logarithmic multiplying factor. 

The observed resistivity does not follow directly from the measured x; a further 
assumption is needed. The assumption that works is that fermions couple to 
the fluctuations locally, as in an effective Hamiltonian ~ c\^a'^i,<T'^i,a,a' ■ Then 
if the measured fluctuation spectra are that of some localized spins Si, the 
single-particle self-energy is that due to exchange of bosons with propagator 
proportional to 

^X(k,a;) - ln(a;) +i sgn (a;). (186) 

k 

This ensures that the single-particle relaxation rate as well as the transport 
relaxation rate is proportional to T. A major theoretical problem is why the 
non-local or "recoil" terms in the interaction of itinerant fermions are irrele- 
vant; i.e why is the effective Hamiltonian not 

^kA^,,a(^kAS, + Sl^).',.- (187) 

Why has momentum conservation been legislated away? 

The singularity of x also raises the question whether the anharmonic processes. 
Fig. 33 (^c^, which are benign and allow the elimination of fermions in the RPA 
theory, give singular contributions to %. Also can fermions really be eliminated 
in calculating the critical behavior. 
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Fig. 44. The phase diagram of UCu^.,-xPdx- At low dopings and temperatures the 
system is in an antiferromgnetic phase. In the undoped sample (x = 0): T/v= 15 K 
with a magnetization ^ ~ l/U^. For doping x = 1 and 1.5 the specific heat C^jT ~ T 
(for X =V) while displaying weak logarithmic characteristics for x = 1.5. Similarly, 
the susceptibility x(^) ^ InT and % ^ T~°-^^ (for a; = 1, 1.5 respectively). Courtesy 
of M. C. Aronson. 

As an example in the second category, we show in Fig. 44 the phase diagram of 
and some of the resistivity data of the heavy fermion compound UCu^-xPdx- 
There are several other compounds in this category also; for a review we refer 
to [160,165]. For a theoretical discussion of the scaling properties of some of 
this class of problems see [27]. 

For X < 1, there is an (antiferromagnetic) ordered state at low temperatures, 
while for X > 2, a spin-glass phase appears. At first sight, one would therefore 
expect possible SFL behavior only near the critical composition x = 1 and 
near x = 2. The remarkable observation however is that over a whole range 
of intermediate compositions, one observes anomalous behavior of the type 



The data for x that show this power law behavior for U Cu^ ^Pdi ^ and U CuiPd 
are shown in Fig. 45. The fact that this is genuine scaling behavior is inde- 
pendently confirmed [29,30] from the fact that the frequency dependent sus- 
ceptibility, measured by neutron scattering, shows a very good collapse of the 
data with the scaling assumption 



In passing, we note that this as well as the result for x in CeCuAu are examples 
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Fig. 45. The temperature dependence of the static susceptibihty of Xo{T) for both 
UCu^Pd and UCu^,^Pdi,^, showing that for both compounds ^o(^) has a low tem- 
perature divergence as T~^/^. The measuring field is 1 Tesla. From Aronson et al. 
[29,30]. 

of an anomalous dimension, as the engineering dimension of the susceptibility 
X is 1/ energy — see the remark made just after Eq. (172). 



Again, none of this behavior finds a clear explanation in any of the well-studied 
models. One is tempted to use the critical points of impurity models (see for 
example [69] and references therein), but runs into the difhculty of having 
to tune to special symmetries. The ideas of critical points of metallic spin- 
glasses (see for example [209] and also [219]), although theoretically appealing, 
are also not obviously applicable over such a wide range of composition. It 
must be mentioned however that NMR evidence does show clear evidence 
for the inhomogeneity in the singular part of the magnetic fluctuations in 
several heavy fermion compounds [41]. This has inspired models of varying 
sophistication (see [173] and references therein, and also [54]) in which the 
Kondo-temperature itself has a inhomogeneous distribution. It is possible to 
fit the properties with reasonable distributions but there is room for a deeper 
examination of the theoretical issues related to competition of disorder, Kondo 
effects and magnetic-interaction between the magnetic moments. 



6.4 Special complications in heavy fermion physics 



In heavy fermion compounds, there is often an additional complication that 
besets treating a QCP as a simple antiferromagnetic transition coupled to itin- 
erant electrons. Often such materials exhibit magnetic order of the /-electrons 
(with magnetic moments of the 0{1^b) per /-electron). Thus, such materi- 
als have local moments in the ordered phase; so the disappearance of the 
(anti)ferromagnetic order at a quantum critical point is accompanied by a 
metal-insulator transition of the /-electrons. This means that the volume of 
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Fig. 46. Schematic occupation number function for i/Pta and for UPd^. 
the Fermi-surface changes in the transition. 

We may illustrate the above scenario by comparing UPts, a Heavy fermion 
compound with effective mass of the order of 100, with UPd^, an "ordinary" 
metal with effective mass of 0{1) in which the f-electrons are localized. A 
schematic summary of the momentum occupation rik for the two cases is 
shown in Fig. 46: in the former Uk is shown with two discontinuities, one 
small C(10~^) representing the large renormalisation in the effective mass of 
"/-electrons" while the other is close to 1 representing the modest renormal- 
isation of s and (i-electrons. The other case, representative of UPd^ has just 
one Fermi surface with a jump in of close to 1. The Fermi surface in the 
former encloses the number of electrons equal to the sum of the / and s-d elec- 
trons while the latter includes only the s-d electrons. This is consistent with 
de Haas-van Alphen measurements as well as the band structures calcula- 
tions of the two compounds; but the band-structure calculation must be done 
with the f-electrons assumed as itinerant in the former and part of the local- 
ized core in the latter. The magnetic transitions in heavy fermion compounds 
(with ordered moment of 0{1hb)) occurs through the conversion of itinerant 
/-electrons to locahzed electrons. So the Fermi surface on the two sides of 
the transition must switch between the two schematic representations in Fig. 
46. The problem couples the "metal-insulator transition" of the /-electrons to 
the magnetic fluctuations — those of itinerant electrons on one side and of in- 
teracting local moments on the other. The fluctuations of the metal-insulator 
transition and the Fermi surfaces is an important part of the problem. Some 
theoretical work with these ideas in mind is available [227]. Another possible 
approach is to generate an effective Hamiltonian for the heavy fermion lattice 
from a pairwise sum of the effective Hamiltonian deduced from the two-Kondo 
impurity problems discussed in section 3.9 and study its instabilities. The two 
impurity problem contains the rudiments of some of the essential physics. 

In connection with the data in CeCuQ-x^Ux, we have discussed two important 
puzzles: the non-trivial exponent dM/z measured by x(A;,a;), and the coupling 
of fermions to the local fluctuations alone for transport properties. In the other 
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category (impurity-dominated) the first puzzle reoccurs; the second puzzle 
may be explained more easily since the measured fluctuation spectrum x{k, ui) 
is in fact fc— independent. Both puzzles reoccur in the SFL phenomena in the 
cuprate compounds to be discussed in section 7. 

6.5 Effects of impurities on Quantum Critical Points 

As is well known, randomness can have an important effect on classical phase 
transitions. Two classes of quenched disorder are distinguished: First, im- 
purities coupling quadratically to the order parameter [116] or, equivalently, 
impurities which may be used to define a local transition temperature Tc(r); 
the second class concerns impurities coupling linearly to the order parameter 
[123]. The so-called Harris criterion, for the first class, tells us that the disorder 
is relevant, i.e. changes the exponents or turns the transition to a crossover, if 
the specific heat exponent a of the pure system is positive or, equivalently, if 

du-2<0 . (190) 

For application to QCP phenomena the value of u to be used is different 
in the quantum fluctuation regime and the quasi-classical regimeP^ For the 
latter u should be defined by the correlation length ^ ~ (T — Tc)~'^^ for a 
fixed {p — pc) while deep in the former, near T = it should be defined by 
^ ~ (p ~ Pc)~'^^- Accordingly the effect of disorder depends on the direction 
one approaches the QCP. Similarly the celebrated Imry-Ma argument [123] 
for linearly-coupled disorder can be generalized to QCP's. 

In problems of fermions, additional effects of disorder arise because the ver- 
tices coupling the impurity to the fermions can be renormalized due to the 
singularity in the fluctuation of the pure system [256]. Not too much work 
has been done along these lines. A simple example is the effect of magnetic 
impurities near a ferromagnetic transition [147] . The growth of the magnetic 
correlation length converts a single-channel Kondo effect to a multi-channel 
Kondo effect with a regime in which the singularities discussed in section 3.8 
for the degenerate multi-channel Kondo effect may be realized without tun- 
ing any parameters [157]. This may be relevant to the deviations from the 
predictions of the pure case discussed here in the properties near the QCP in 
MnSi. Extensions of these ideas to antiferromagnetic and other QCP's would 
be quite worthwhile. 

Actually, the Harris criterion is derived in the form (190) for v and d, not in terms 
of the exponent a. This is particularly important at QCP's, since as we discussed 
in section 6.2 for QCP's one is often above the upper critical dimension where the 
hyperscaling relation dv — 2 = a breaks down. 
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Fig. 47. (a) Schematic structure of the copper-oxide a6-planes in La2Cu04. Ba 
or Sr substitution for La in the parent compound La2Cu04 introduces holes in 
the CuOi planes. The structure of other high Tg materials differs only in ways 
which do not affect the central issues, e.g., it is oxygen doping in YBa2Cu306-^.x 
that provides planar holes. The magnetic moments of the planar copper atoms 
are ordered antiferromagnetically in the ground state of the undoped compounds. 
Prom [127]. (b) The "orbital unit cell" of the Cu-O compounds in the ab plane. 
The minimal orbital set contains a dj.2_y2 orbital of Cu and Px and a Py orbital of 
oxygen per unit cell. 

7 The High-Tc Problem in the Copper-Oxide Based Compounds 

About 10^ scientific papers have appeared in the field of high- Tc superconduc- 
tors since their discovery in 1987. For reviews, sec the Proceedings [164] of the 
latest in a series of Tri- annual Conferences or [106]. Although no consensus 
on the theory of the phenomena has been arrived at, the intensive investi- 
gation has resulted in a body of consistent experimental information. Here 
we emphasize only those properties which are common to all members of the 
copper-oxide family and in which Singular Fermi-liquid properties appear to 
play the governing role. 

The high-Tc materials are complicated, and many fundamental condensed mat- 
ter physics phenomena play a role in some or other part of their phase diagram. 
As we shall see, the normal state near the composition of the highest Tc shows 
convincing evidence of being a weak form of a SFL, a Marginal Fermi Liquid. 
Since the vertices coupling fermions to the fluctuations for transport in the 
normal state and those for Cooper pairing through exchange of fluctuations 
can be derived from each other, the physics of SFL and the mechanism for 
superconductivity in the cuprates are intimately related. 
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Fig. 48. a) The crystal structure of La2Cu04. From [187]. Electronic couplings 
along the c direction are weak, b) Schematic of the Cu02 plane. The arrows indicate 
the alignment of spins in the antiferromagnetic ground state of La2Cu04. Speckeled 
shading indicates oxygen orbitals; coupling through these leads to superexchange 
in the insulating state. 

7.1 Some basic features of the high-T^ materials 

A wide variety of Cu-0 containing compounds with different chemical for- 
mulae and different structures belong in the high-Tg family. The common 
chemical and structural features are that they all contain two-dimensional 
stacks of Cu02 planes which are negatively charged with neutralizing ions 
and other structures in between the planes. The minimal information about 
the structure in the Cu — O planes and the important electronic orbitals on 
Cu and on Oxygen is shown in Fig. 47. The structure of one of the simpler 
compound La2-xSrxCuOi is shown in Fig. A&(a) with the Cu02 plane shown 
again in Fig. A%(h). For x = the Cu02 plane has a negative charge of 2 per 
unit cell which is nominally ascribed to the Cm^^(0^~)2 ionic configuration. 
Since has a filled shell while Cv?^ has a hole in the 3d shell, the Cu-0 
planes have a half-filled band according to the non-interacting or one-electron 
model. But at x = the compound is an antiferromagnetic insulator with 
S — 1/2 at the Copper sites. This is well known to be characteristic of a 
Mott-type insulator in which the electron-electron interactions determine the 
ground state. Actually [282,251,252] copper-oxide compounds at x = belong 
to the Charge- Transfer sub-category of Mott- insulators. But at a; = the 
ground state and low-energy properties of all Mott-insulators are qualitatively 
the same. By substituting divalent Sr for the trivalent La in the above exam- 
ple, "holes" are introduced in the copper-oxide planes with density x per unit 
cell. 

Fig. 49 is the generic phase diagram of the Cu-O compounds in the T-x 



105 



T 



N 



I 



N SlofuUx Feraol litviild y 



U 




\ / m 




Fig. 49. Generic phase diagram of the cuprates for hole doping. The portion labeled 
by AFM is the antiferromagnetic phase, and the dome marked by SC is the su- 
perconducting phase. Crossovers to other characteristic properties are marked and 
discussed in the text. Not shown is a low-temperature "insulating phase" in region 
II due to disorder. 

plane. In the few compounds with electron doping which have been synthesized 
properties vary with doping density in a similar way. 

Antiferromagnctism disappears for x typically less than 0.05 to be replaced 
by a superconducting ground state starting at somewhat larger x. The super- 
conducting transition temperature is peaked for x typically between 0.15 and 
0.20 and disappears for x typically less than 0.25. We will define Xm to be 
the density for maximum T^. Conventionally, copper-oxides with x < Xm are 
referred to as underdoped, with x = x„i as optimally doped and with x > Xm 
as overdoped. Superconductivity is of the "c?-wave" singlet symmetry. 

The superconducting region in the T-x plane is surrounded by three distinct 

regions: a region marked (111) with properties characteristic of a Landau Fermi- 
liquid, a region marked (1) in which (marginally) SFL properties are observed 
and a region marked (11) which is often called the pseudo gap region whose cor- 
relations in the ground state still remain a matter of conjecture. The topology 
of Fig. 49 around the superconducting region is that expected around a QCP. 
Indeed it resembles the phase diagram of some heavy fermion superconductors 
(sec e.g. Fig. 39) except that region 11 has no long-range antiferromagnetic or- 
der — the best experimental information is that, generically, spin rotational 
invariance as well as (lattice) translational invariance remains unbroken in the 
passage from (I) to (II) in the Cu-O compounds. 

The quantity 7(T) = Cv{T) /T and the magnetic susceptibility x{T), which are 
temperature independent for a Landau-Fermi-liquid begin to decline rapidly 
[154] in the passage from region I to region II, which we will call Tp{x), but 
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Fig. 50. The electronic contribution to the specific heat as a function of temperature 
for underdoped, optimally doped and overdoped samples of Yo,8Cao.2Ba2Cu307_a;. 
For optimally doped and overdoped samples the heat capacity remains constant 
as the temperature is lowered, then shows the characteristic features at the super- 
conducting transition temperature Tc and rapidly approaches zero in the supercon- 
ducting state. For underdoped samples, however, the heat capacity starts to fall 
well above Tc as the temperature is reduced, and there is only a small peak at 
Tc indicating much smaller condensation energy than the optimal and overdoped 
compounds [154]. From [34]. 

without any singular feature. However, the transport properties — resistivity. 
Nuclear relaxation rate (NMR), etc. — show sharper change in their tem- 
perature dependence at Tp{x). The generic deduced electronic contribution to 
the specific heat for overdoped, optimally-doped and underdoped compounds 
is shown in Fig. 50. The generic behavior for an underdoped compound for 
the resisitivity. Nuclear relaxation rate and Knight shift — proportional to 
the uniform susceptibility — is shown in Fig. 51. ARPES (Angle Resolved 
Photoemission) measurements show a diminution of the electronic density of 
states starting at about Tp{x) with a four- fold symmetry: no change along 
the (tt, tt) directions and maximum change along the {tt, 0) directions. The 
magnitude of the anisotropic "pseudogap" is several times Tp{x). 

It is important to note that given the observed change in the single-particle 
spectra, the measured specific heat and the magnetic susceptibility in the 
pseudogap region arc consistent with the simple calculation using the single- 
particle density of states alone. Nothing fancier is demanded by the data, at 
least at in its present state. Moreover, the transport properties as well as the 
thermodynamic properties at different x can be collapsed to scaling functions 
with the same Tp{x) as a parameter [275]. 
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Fig. 51. Signatures of the pseudogap in various transport properties for the un- 
derdoped compound YBa2Cu40s- At high temperatures the resistivity (sohd hne) 
decreases hnearly with temperature. In the pseudogap region it drops faster with 
temperature before faUing to zero at the superconducting transition temperature 
(about 85K) . Similarly the NMR relaxation time displays characteristics of the op- 
timum doped compounds above about 200 K (squares on dashed line) but deviates 
strongly from it in the pseudogap region. The NMR shift (top squares) also de- 
viates from the temperature-independent behavior (not shown) below the inset of 
the pseudogap. Note that the pseudogap expresses itself as a sharper change with 
temperature in the transport properties compared to the equilibrium properties - 
specific heat and magnetic susceptibility. [45,278,10]. Prom [34]. 

7.2 Marginal Fermi Liquid behavior of the normal state 

Every measured transport property in Region I is unlike those of Landau 
Fermi-liquids. The most commonly measured of these is the dc resistivity 
shown for many different compounds at the "optimum" composition in Fig. 
52 including one with Tc ~ lOK. The resistivity is linear from near Tc to 
the decomposition temperature of the compound. As shown in Fig. 51, in the 
"under-doped" region the departure from linearity commences at a temper- 
ature ~ Tp{x) marked in Fig. 49. Similarly, the cross-over into region (III) 
shown in Fig. 49 is accompanied by Fermi-liquid like properties. Wherever 
measurements are available, every other measured transport property shows 
similar changes. 

The different measured transport properties study the response of the com- 
pounds over quite different momentum and energy regions. For example, the 
Raman scattering studies the long wavelength density and current response at 
long wavelength but over a range of frequencies from low 0{lcm^^) to high, 
(9(10^cm~^). On the other hand nuclear relaxation rate ^ depends on the 
magnetic fluctuations at very low frequencies but integrates over all momenta, 
so that the short wave-length fluctuations dominate. In 1989 it was proposed 
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Fig. 52. Resistivity as a function of temperature for various high temperature su- 
perconductors. Prom [38]. 

[253,254] that a single hypothesis about the particle-holc excitation spectra 
captures most of the diverse transport anomalies. The hypothesis is that the 
density as well as magnetic fluctuation spectrum has an absorptive part with 
the following property: 



= -X^' oj/T , for u; « T , 
—XoSgn{uj) , for Uc > uj T 



(191) 



Here x" order of the bare single-particle density of states A^(0) and 

cUc is an upper cut-off. The fluctuation spectrum is assumed to have only 
a weak momentum dependence, except at very long wavelength, where a 
dependence is required for fluctuations of conserved quantities like density or 
spin (in the absence of spin-orbit interactions). A form which implements these 
requirements for the conserved quantities with a rather arbitrary crossover 
function to get the different regimes of lo/T is 

2 

X"(q,^) ~ . o , — for VFq<^yA^, (192) 



where a; = a; for uj/T < T and a; = T for cu/T > 1. 

Using the Kramers-Kronig relations, one deduces that the real part of the 
corresponding correlation functions have a log(a;) divergence at all momenta 
except the conserved quantities where the divergence does not extend to vpq > 
X. Thus compressibility and magnetic susceptibility are flnite. (Besides (191), 
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(192) analytic fluctuation spectra of the Fermi-liquid form is of course also 
present.) 

The spectral function (191) is unlike those for a Landau Fermi-liquid dis- 
cussed in section 2, which always have a scale — the Fermi-energy, Debye- 
frequency, or spin-wave energy, etc. — which is obtained from parameters 
of the Hamiltonian. Such parameters have been replaced by T. As we have 
discussed in section 6, this scale- invariance of (191) is characteristic of fluctu- 
ations in the quasi-classical regime of a QCP. Eq. (191) characterizes the fluc- 
tuations around the QCP: comparing with Eq. (174), the exponent (Im/z — 
and 1/z = 0. These are equivalent to the statement that the momentum de- 
pendence is negligibly important compared to the frequency dependence. This 
is a very unusual requirements for a QCP in an itinerant problem: the spatial 
correlation length plays no role in determining the frequency dependence of 
the critical properties. 

The experimental results for the various transport properties for compositions 
near the optimum are consistent in detail with Eq. (191), supplemented with 
the elastic scattering rate due to impurities (see later). We refer the reader 
to the original literature for the details. The temperature independence and 
the frequency independence in Raman scattering intensity up to ou of 0(leV) 
directly follows from (192). Eq. (192) also gives a temperature independent 
contribution to the nuclear relaxation rate Tf^ as is observed for Cu nu- 
clei. The transport scattering rates have the same temperature dependence as 
the single-particle scattering rate. The observed anomalous optical conductiv- 
ity can be directly obtained by using the continuity equation together with 
Eq. (192), or by first calculating the single-particle scattering rate and the 
transport scattering rate. The single-particle scattering rate is independently 
measured in ARPES experiments and provides the most detailed test of the 
assumed hypothesis. 

To calculate the single-particle scattering rate, assume to begin with a con- 
stant coupling matrix element g for the scattering of particles by the singular 
fluctuations. We shall return to this point in the section on microscopic the- 
ory. Then provided there is no singular contribution to the self-energy from 
particle-particle fluctuations, the graph in Fig. 53 represents the singular self- 
energy exactly. It is important to note that for this to be true, Eq. (192) 
is to be regarded as the exact (not irreducible) propagator for particle-hole 
fluctuations; it should not be iterated. 

E(q, cu) can now be evaluated straight forwardly to flnd a singular contribu- 
tion. 
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Fig. 53. (a) Diagram for the singular contribution to the one-particle self-energy 
with the fluctuating x{q^^)- ff's are the vertices which in microscopic theory [255] 
is shown to have important momentum dependence, but which gives negligible mo- 
mentum dependence to the self-energy, (b) For XoiQ^'^) which is momentum inde- 
pendent, a total vertex F may sometimes be usefully defined, which has the same 
frequency dependence as Eq. (193) and which is also ^'-independent. The lines are 
the exact single-particle Green's functions. 

for X <^ cOc and vp\{q — kp)\ ^ Uc- The noteworthy points about (193) are: 

(1) The single particle scattering rate is proportional to x rather than to x"^ 
as in Landau Fermi-liquids. 

(2) The momentum independence of the single-particle scattering rate. 

(3) The quasi-particle renormalization amplitude 

fl-Aln —V' (194) 



scales to zero logarithmically as x 0. Hence the name Marginal Fermi- 
liquid. 

(4) The single-particle Green's function 

G{u;,q) = J ^—-r, r (195) 



has a branch cut rather than a pole. It may be written as 

Z{x) 



UJ- {Iq- p) - i/f ' 



(196) 



where Iq is the renormalized single-particle energy: 

Iq-fi ^Z{eq-fi) ^ ZYF-icL-kp) (197) 

for small \q — kp]- Also f~^{x) — Z ImE(a;), and the effective Fermi-velocity 

vp = Zvp has a frequency and temperature dependent correction. 
(5) The single-particle occupation number has no discontinuity at the Fermi- 
surface, but its derivative does, see Fig. 54. So the Fermi-surface remains a 
well-defined concept both in energy and in momentum space. 
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Fig. 54. The T = distribution of bare particles in a marginal Fermi liquid. No 
discontinuity exists at kp but the derivative of the distribution is discontinuous. 

The predictions of (195) have been tested in detail in ARPES measurements 
only recently. ARPES measures the spectral function 



In ARPES experiments, the energy distribution curve at fixed momentum 
(EDC) and the momentum distribution curve at fixed energy (MDC) can 
both be measured. It follows from Eq. (198) that if S is momentum inde- 
pendent perpendicular to the Fermi-surface, then an MDC scanned along 
for a; ~ // should have a Lorentzian shape plotted against (k — kp)± with a 
width proportional to E"(a;)/vF(k). For a Marginal Fermi Liquid (MEL), this 
width should be proportional to x. The agreement of the measured lineshape 
to a Lorentzian and the variation of the width with temperature are shown in 
Figs. 55. Fig. 56 shows the momentum width measured in different directions 
k. The Fermi- velocity v(k) is measured through the EDC with the conclusion 
that it is independent of k to within the experimental errors. The data shows 
that the temperature dependence is consistent with linearity with a coefficient 
independent of k although the error bars are huge near the (vr, 0) direction. 
Besides the MFL contribution, there is also a temperature independent contri- 
bution to the width which is strongly angle-dependent, to which we will soon 
turn. The ambiguity of the temperature (and frequency) dependence near the 
(tt, 0) direction is removed by the EDC measurements. In Fig. 56 EDC at 
the Fermi-surface in the (vr, tt) direction and the (vr, 0) directions are shown 
together with a fit to the MFL spectral function with a constant contribution 
added to E". EDC's have the additional problem of an energy-independent 
experimental background. This has also been added in the fit. In both direc- 
tions S" has a contribution proportional to uj with the same coefficient within 
the experimental uncertainty. In summary, the ARPES experiments give that 



A{q,u;) 



1 ^"{q,io) 

TT [a;-(e,-/x)-E'(g,u;)]2+[E"(g,u;)]2 " 



(198) 



i:"{k,u]T) 



r(M +T."^p^{uj,T). 



(199) 
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Fig. 55. (a) Momentum distribution curves for different temperatures. The solid 
lines are Lorentzian fits, (b) Momentum widths of MDCs for three samples (circles, 
squares, and diamonds). The thin lines are T-linear fits which show consistency with 
Eqn.(198) and the MFL hypothesis. The resistivity (solid black line) is also shown. 
The double-headed arrow shows the momentum resolution of the experiment. Prom 
Johnson et al. [246]. 



The (a;, T)-independent contribution r(ki;') can only be understood as due 
to impurity scattering [2]. Its dependence on can be understood by the 
assumption that in well-prepared samples, the impurities lie between the Cu- 
O planes and therefore only lead to small angle scattering of electrons in the 
plane. The contribution r(k^) at k^ then depends on the forward scattering 
matrix element and the local density of states at k^;' which increase from the 
(tt, tt) direction to the (tt, 0) direction. This small angle contribution has several 
very important consequences: (i) relative insensitivity of residual resistivity 
to disorder, (ii) relative insensitivity of (i-wave superconductivity transition 
temperature to the elastic part of the single-particle scattering rate [134], 
and (in) most significantly relative insensitivity to the anomalous Hall effect 
and magneto-resistance. Such anomalous magneto-transport properties follow 
from a proper solution of the Boltzmann equation including both small-angle 
elastic scattering and angle- independent MFL inelastic scattering [260] . 

The momentum independence of the inelastic part of S" is crucial to the SFL 
properties of the cuprates. Because the inelastic scattering to all angles on 
the Fermi-surface is the same, i.e. s-wave scattering, the vertex corrections 
to transport of vector quantities like particle current and energy current are 
zero. It follows that the momentum transport scattering rate measured in 
resistivity or optical conductivity and the energy transport rate measured in 
thermal conductivity have the same (cc',T) dependence as the single-particle 
scattering rate 1/t{ijj,T). 
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Fig. 56. Fits of the MFL self-energy r + A?ia; to the experimental data, according to 
(199). Estimated uncertainties of ±15% in F and ±25% in A. (a) A scan along the 
(1,0) direction, F = 0.12, A = 0.27; (b) A scan along the (1,1) direction, F = 0.035, 
A = 0.35. Prom Kaminsky and co-workers, [131]. 




Fig. 57. Left panel: Fermi surface of the optimally doped Bi2212, measure in the 
superconducting state. Indicated are the lines (1) to (5) on which the temperature 
dependence is measured in the experiments by Valla et al. [247] . Typical spectra are 
shown for line (2) in the normal (b) and superconducting state (c), as well as for 
line (4), in the normal (d) and superconducting state (c). Right panel: momentum 
widths as a function of temperature for different positions on the Fermi surface, 
obtained by fitting the Momentum Distribution Curves with Lorentzian lineshapes 
in the same experiments. Widths are measured at the Fermi level and at the leading 
edge, in the normal and in superconducting (grey region) state, respectively. 
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Recently far-infrared condTictivity measurements [67] have been analyzed and 
shown to be consistent with l/r(a;,T) deduced from MFL including the log- 
arithmic corrections. 

As already discussed no singular correction to the magnetic susceptibility is 
to be expected on the basis of (191). But the specific heat should have a 
logarithmic correction so that 

7(T) =7o[l + Aln(u;e/T)] . (200) 

Such a logarithmic correction has not yet been deciphered in the data pre- 
sumably because the electronic specific heat in the normal state is less than 
C(10~^) of the total measured specific heat and must be extracted by a sub- 
traction procedure which is not sufficiently accurate. 



1.3 General Requirements on a Microscopic Theory 

The MFL self-energy, Eq. (193), has been verified in such detail in its (a;, T, q) 
dependence that it is hard to see how any theory of CuO compounds can be 
relevant to the experiments without reproducing it (or a very close approxima- 
tion to it) in Region I of the phase diagram of Fig. 49. Such a scale-invariant 
self-energy is characteristic of the quasi-classical regime of a QCP and indeed 
the topological features of the phase diagram are consistent with their being 
a QCP at Xc near the composition for the highest Tc (Alternatively, a QCP in 
the overdoped region where Tc vanishes is predicted in some approaches, like in 
[220]). To date, no way has been found to obtain Eq. (192) except through the 
scale-invariant form of fluctuations which is momentum independent {z ~ oo) 
as in Eq. (191). 

A consistent and applicable microscopic theory of the copper-oxides must show 
a QCP with fluctuations of the form (192). This is a very unusual requirement 
for a QCP in a homogeneous extended system for at least two reasons. First, 
the fluctuations must have a negligible q-dependence compared to the fre- 
quency dependence, i.e. z ^ oo and second, the singularity in the spectrum 
should just have logarithmic form; i.e. there should be no exponentiation of 
the logarithm giving rise to power laws. Such singularities do arise, as we dis- 
cussed in section 3, in models of isolated impurities under certain conditions 
but they disappear when the impurities are coupled; recoil kills the singulari- 
ties. The requirement of negligible q/uj dependence runs contrary to the idea 
of critical slowing down in the fluctuation regime of the usual transitions, in 
which the frequency dependence of the fluctuations becomes strongly peaked 
at zero-frequency because the spatial correlation length diverges. 
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Another crucial thing to note is that any known QCP (in more than one 
dimension) is the end of a hne of continuous transitions at T = 0. Region II 
(at least at T = 0) must then have a broken symmetry (This includes part of 
Region III, which is also superconducting). Experiments appear to exclude any 
broken translational symmetry or spin-rotational symmetry for this regionj^ 
although as discussed in section 7.1, a sharp change in transport properties is 
observed along with a four-fold symmetric diminution of the ARPES intensity 
for low energies aX T ^ Tp{x). If there is indeed a broken symmetry, it is of a 
very elusive kind; experiments have not yet found it. 

A related question is how a momentum-independent x(q, co',T) can be the 
fluctuation spectrum of a transition which leads to an anisotropic state as 
in Region II. Furthermore, how can such a spectrum lead to an anisotropic 
superconducting state because it is unavoidable that xi^^^^T) of Eq. (192) 
which determines the inelastic properties in Region I also be responsible for the 
superconductive instability. After all, the process leading to the normal self- 
energy. Fig. 53, the superconductive self-energy, and the Cooper pair vertex. 
Fig. 58, are intimately related. Given Imx(q, Ci;,T) the effective interaction in 
the particle-particle channel is 



Rex(g; uj) is negative for all q and for all —Uc < uj < ujc- So we do have a mech- 
anism for superconductive pairing in the Cu-0 problem given by the normal 
state properties just as the normal state self-energy and transport properties 
of say Ph tell us about the mechanism for superconductivity in Ph. In fact 
given that the normal state properties give that the upper cut-of frequency is 
of (9(10^)K and the coupling constant A ~ g'^N{0) is of C(l), the correct scale 
of Tc is obtained. The important puzzle is, how can this mechanism produce 
(i-wave pairing given that x(q, a;) is momentum independent. How can one 
obtain momentum-independent inelastic self-energy in the normal state and a 
(i-wave superconducting order parameter from the same fluctuations? 

We summarize in the next section a microscopic theory which attempts to 
meet these requirements and answer some of the questions raised. 



A new lattice symmetry due to lattice distortions or antiferromagnetism, if signif- 
icant, would change the fermi-surface because the size of the Brillouin zone would 
decrease. This would be visible both in ARPES measurements as well as in hall 
effect measurements. 




(201) 
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Fig. 58. The Cooper pair-vertex and the normal state self-energy, Fig. 53, are inti- 
mately related. 

7.4 Microscopic Theory 



There is no consensus on even the minimum necessary model Hamiltonian 
to describe the essential properties in the phase-diagram, Fig. 49, of the CuO 
compounds. It is generally agreed that, since other transition metal compounds 
do not share the properties of CuO compounds, a model Hamiltonian with 
some rather special features is called for. Two such features are 

(A) They are two-dimensional with an insulating antiferromagnetic ground 
state and spin S = 1/2 per unit cell at half filling. Although not unique this 
feature is rare. If this is the determining feature, a two-dimensional Hubbard 
model is adequate. [24]. Even this model is not soluble in ci = 2. 

(B) CuO are the extreme limit of charge transfer compounds [282] in which 
charge fluctuations in the metallic state occur almost equally on Oxygen and 
Copper. Then the longer-range ionic interactions, which in magnitude are 
comparable to the on-site interactions have a crucial role to play in the low- 
energy dynamics in the metallic state through excitonic effects. A model with 
both Cu and O orbitals, hopping between them, and the excitonic interactions 
besides the on-site repulsions is then required [251,252]. This is of course even 
harder to deal with than the Hubbard model. 

Numerous attempts have been made using one or the other such models to 
obtain SFL behavior. We briefly discuss the motivations for the pursuit of 
model (A) before summarizing in a little more detail the only attempt to 
obtain the phenomenlogical fluctuation spectrum of Eq. (191) , and which relies 
on a model of type (B). 



7.4.1 The Doped Hubbard Model 

The investigations of the copper-oxide problem from this point of view asks 
some valid and deep questions [24]. How does a low concentration of holes move 
through the spin configurations in a two-dimensional model with restriction 
of zero double occupancy? We have sketched in section 4.8 the difficulties of 
connecting to the same problem in one dimension when the ground state at 
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zero doping is an antiferromagnet. In fact analytic [213] and numerical [162] 
answers to the question for a single hole show the spectral weight of a heavy 
particle with an incoherent part composed of multiple spin-wave polaronic 
cloud. Simply extrapolated (a dangerous thing to do), a Fermi-liquid is ex- 
pected. The larger zero-point fluctuations of the S — 1/2 model, compared 
to a large spin model only change the relative weight of the coherent and 
the incoherent parts. But more subtle possibilities must be considered. The 
antiferromagnetic ground state of a Heisenberg S = 1/2 model (or the un- 
doped non-degenerate Hubbard model) in two-dimensions is close in energy 
to a singlet ground state. A possible description of such a state is as linear 
combination on a basis of singlet-bonds of pairs of spins. As noted earlier, such 
itinerant bond states have been termed Resonating Valence Bonds (by analogy 
to the ground state of Benzene like molecules). The massive degeneracy of the 
singlet bond-basis raises interesting possibilities. If the quantum fluctuations 
of spins were (significantly) larger than allowed by = 1/2, such states would 
indeed be the ground state, as they arc in the one-dimensional model or two- 
dimensional models with additional frustrating interactions [6]. It is possible 
that by doping with holes in the S = 1/2 Heisenberg model, the additional 
quantum fiuctuations induce a ground state and low-lying excitations which 
utilize the massive degeneracy of RVB states. Especially intriguing is the fact 
that resonating valence ground state may be looked on as the projection of 
the BCS ground state to a fixed number of particles [24]. Furthermore in the 
normal state this line of reasoning is likely to lead to a SFL. 

A specific proposal incorporating the RVB idea [136] rehes on the ground state 
of the half-filled model to be localised dimers. Then defects in this state due to 
deviation from half-filling can plausibly support excitations which are charged 
spinless bosons. Further work on this idea may be found in [210,176]. Related 
ideas were put forth in [80]. 

These are very attractive set of ideas and no proof exists that they are disal- 
lowed. We have already considered an implementation of these ideas in section 
5.2 on generalized gauge theories. As discussed, controlled calculations using 
these ideas are hard to come by. Moreover, what theoretical results do ex- 
ist do not correspond in a persuasive way to the experimental results on the 
copper-oxide materials. 

One should take special note here of the idea of Anyon superconductivity 
which besides being a lovely theoretical idea, is founded on the solution to 
a well defined model, and has clear experimental predictions. Laughlin and 
collaborators [148,130,89] found a specific model with long-range four spin 
interactions for which his quantum hall wavefunction is the ground state. 
Therefore Time-reversal and Parity are spontaneously broken in this state. 
This state is superconducting. The predicted Time-reversal broken properties 
have not been observed experimentally [233]. 



118 



Fig. 59. Three bands result from the orbitals shown in Fig. 4:7 (b) in a one-electron 
calculation; two of these are shown. The chemical potential lies in the "anti-bonding" 
band and is varied by the doping concentration. The other band shown is crucial 
for the theory using excitonic effects as in [255]. 

An alternative idea from the microscopic characterization of these materi- 
als as doped Hubbard models is that a dilute concentration of holes in the 
Hubbard model is likely to phase separate or form ordered one-dimensional 
charge- density wave/spin- density wave structures [283,284,84]. There exists 
both empirical [242] and computational [272] support for this idea at least for 
a very dilute concentration of holes. For concentrations close to optimum com- 
positions these structures appear in experiment to exist only at high energies 
with short correlation lengths and times and small amplitudes in the majority 
of copper-oxygen compounds. Their relation to SFL properties is again not 
theoretically or empirically persuasive. 

7.4-2 The excitonic Interactions Model 

This relics on a model of type {B). A brief sketch of the calculations leading to 
a QCP and a MFL spectrum is given here. We refer the reader to [251,253,260] 
for details. 

At half-filling the ground state and the low-lying excitations of such models 
are identical to the Hubbard model. But important differences can arise in the 
metaUic state. Consider the one-electron structure of such models. The O-O 
hopping in the lattice structure with d^^^yi orbitals in Cu and Px, Py orbitals 
on O, as shown in Fig. 47, produces a weakly dispersing "non-bonding" band 
while the Cu-0 hopping produces "bonding" and "anti-bonding" bands — see 
Fig. 59. We need consider only the filled non-bonding band and the partially 
filled anti-bonding bands shown in Fig. 59. 

In the mean-field approximation such an electronic structure together with 
the on-site interaction and the ionic interactions is unstable to an unusual 
phase provided the latter summed over the nearest neighbors is of the order or 
larger than the bandwidth. In this phase translational symmetry is preserved 
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Fig. 60. The current pattern in the time-reversal breaking phase predicted for Region 
II of the phase diagram. 

but time-reversal symmetry and the four-fold rotation symmetry about the 
Cu sites is broken. The ground state has a current pattern, sketched in Fig. 
60, in which each unit cell breaks up into four plaquettes with currents in 
the direction shown. The variation of the transition temperature with doping 
X is similar to the line Tp{x) in Fig. 38, so that there is a QCP at x = Xc- 
Experiments have been proposed to look for the elusive broken symmetry 
sketched in Fig. 60 [259]. 

The long-range interactions in the model also favor other time-reversal break- 
ing phases which also break translational invariance. This is known from cal- 
culations on ladder models [189]. Such states have also been proposed for 
copper-oxide compounds [140,121,57]. 

Our primary interest here is how the mechanism of transition to such a phase 
produces the particular SFL fluctuation spectrum (192) in Region 1 of the 
phase diagram. The driving mechanism for the transition is the excitonic sin- 
gularity, due to the scattering between the states of the partially filled con- 
duction band c and the valence band v of Fig. 59. This scattering is of course 
what we considered in sections 3.5 and 4.9 for the problem of X-ray edge 
singularities for the case that the interband interaction V in Eq. (55) is small 
and the valence band is dispersionless (i.e., the no recoil case). Actually, the 
problem is exactly soluble for the no- recoil case even for large V [64] . For large 
enough V the energy to create the exciton, cug^, is less than the v-c splitting 
A. The excitonic lineshape is essentially the one sketched in Fig. 19(b) and 
given by Eq. (56) for ur > ujex, but 6 is now the phase shift modulo n which 
is the value required to pull a bound state below A. The excitonic instability 
arises when cupr — > 0. 



120 



(a) No Recofl 



WfthHecoIl 




Fig. 61. Sketch of the development of the particle-hole spectra in the microscopic 
model for the cuprates. 

The effect of a finite mass or recoil on the excitonic spectra is to smoothen the 
edge singularity on the scale of the valence band dispersion between A; = and 
k — kp- The phase shift 5 or the interaction energy V no longer determines 
the low-energy shape of the resonance. V docs determine its location. The 
low-energy fluctuation spectra is determined by the following argument: Let 
us concentrate on q = which is obviously where Imxex{<l, i-^) is largest. The 
absorptive part of a particle-hole spectra must be odd in frequency: 



Im x(q,u;) = -Im x(q, -uj) 



(202) 



As V increases, lmXex{^,<^) must shift its weight towards zero-frequency as 
shown in Fig. 61. Let us continue to denote by \ujf.x\ the characteristic energy 
of the maximum in Im^;. For \uj\ small compared to {uj^xl, Ihix(0,c<j) ~ u while 
for |ct;| large compared to \uJex\, it is very slowly varying up to a cut-off Uc on 
the scale of the Fermi-energy. Then by Kramers-Kronig transform the leading 
term in 

Rex(0, a;) ~ In (a;c/ max(a;a;, to")) . (203) 



For any finite \uJex\, Rex is finite and there is no instability. Only for \ujgx\ 0, 
i.e. Imx(0, lu) — >• sgna;, Rex(0, uj) is singular ~ In |a;| and there is an instability. 
Thus in an excitonic instability of a Fermi-sea with a dispersive valence band, 
the zero-temperature spectrum has the form x(<^, 0) ~ In |a;| -|- i sgn cu at the 
instability. Given a parameter p such that the instability occurs only at pc, i.e. 
^^exip Pc) 0, the generalization for finite temperature T and momentum 
q and p ^ Pc is 



x(q,'^) 



+ ln- 



^max(a;, T, ouexip)) max(a;, T, Uexip)) J 

+aY + {Pc{T)-p)r' . (204) 
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Fig. 62. Singularity of interaction Taaaa between states 'a' near the chemical poten- 
tial generated by the excitonic singularity between the states 'a' of the conduction 
band and states '6' of the valence band. The excitonic singularity is indicated by 
the shaded block. 

Here ujc is the frequency cut-off of 0{A). Since the binding energy is O(leV), 
the size of the exciton, a, is of the order of the lattice constant. The q' de- 
pendence of (204) is negligible compared to the frequency dependence. The 
exponent z is effectively infinite. At p ^ pc to logarithmic accuracy, the above 
expression (204) is identical to the phenomenological hypothesis (192). 

The effective low-energy interaction for states near the chemical potential, 
which is sketched in Fig. 62, is singular when the excitonic resonance is at 
low frequency. Here is an example of the mechanism mentioned under (v) in 
section 2.6 where the irreducible interaction obtained by integrating over non- 
perturbed high-energy states is singular. This is of course only possible when 
the interactions represented by the shaded block in Fig. 62 are large enough. 

In relation to some of the questions raised about the phenemenology at the end 
of the last subsection, the momentum dependence of the vertex coupling the 
low energy fermions to the fluctuations in. Figs. 53 and 58 has been evaluated 
[258]. It is non-local with a form depending on the wave functions of the 
conduction and valence bands and the leading result is 

g{k, k') - (sin(/c^a/2) sin(A;>/2) - sin(A;ya/2) sin(/c^a/2)) . (205) 

Note that at {k — k') = this is proportional to [cos(A;a;a) — cos{kya)]. This is 
intimately related to the d-wave current distribution in the broken-symmetry 
phase predicted for region (II), shown in Fig. 60. Eq. (205) is such that that 
when diagram 53 is evaluated, the intermediate state momentum integration 
makes the self-energy depend very weakly on the incoming momentum. But 
when the pairing kernel of Fig. 58 is evaluated, its momentum dependent and 
exhibits attraction in the d-wecve channel. 

Similarly, as has been shown [258], the vertex of Eq. (205) leads to an anisotropic 
state with properties of the pseudo-gap state of region II below a temperature 
Tp{x). The principal theoretical problem remaining with this point of view is 
that a transition of the Ising class occurs at Tp(a;)at least in mean- field the- 
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ory. This would be accompanied by a feature in the specific heat unhke the 
observations-!^ 

The microscopic theory reviewed above reproduces the principal features of 
the phase diagram Fig. 49 of the copper-oxide superconductors, and of the SFL 
properties. It also gives a mechanism for high temperature superconductivity 
of the right symmetry. Further confidence in the applicability of the theory to 
the cuprates will rest on the observation of the predicted current pattern of 
Fig. 60 in region II of the phase diagram [^. 



8 The Metallic State in Two-Dimensions 

The distinction between metals and insulators and the metal-insulator transi- 
tion has been a central problem in condensed matter physics for seven decades. 
Despite the accumulated theoretical and empirical wisdom acquired over all 
these years, the experimental observation made in 1995 of a metal-insulator 
transition in two-dimensions [141] was a major surprise and is a subject of 
great current controversy. The theoretical work in the 1980's [11,90-92,149,14] 
on disordered interacting electrons pointed to a major unsolved theoretical 
problem in two dimensions. Infrared singularities were discovered in the scat- 
tering amplitudes which scaled to strong coupling where the theory is uncon- 
trolled (The situation is similar to that after the singularities in the one- or 
two-loop approximations in the Kondo problem were discovered revealing a 
fascinating problem without providing the properties of the asymptotic low 
temperature state). However the problem was not pursued and the field lapsed 
till the new experiments came along. 

The 1980's theoretical work shows that this problem naturally belongs as a 
subject in our study of Singular Fermi-Liquids. We will first summarize the 
principal theoretical ideas relevant to the problem before the 1995 experi- 
ments. We then briefiy summarize the principal results of these and subse- 
quent experiments. Reviews of the experiments have appeared in [3,16,17]. 
There are two types of theoretical problems raised: the nature of the metallic 
state and the mechanism of the metal-insulator transition. We will address 
the former and only briefiy touch on the latter. 



^ One might appeal to disorder to round off the transition, but this appears im- 
plausible quantitatively. More likely, the nature of the transition is strongly affected 
by the Fluctuation spectra of the form of Eq. (191) and is unlikely to be of the Ising 
class. 

As already mentioned in section 2.6, ferromagnetism in some compounds has an 
excitonic origin. The dynamics near such a transition should also exhibit features 
of the edge-singularity as modified by recoil. 
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Fig. 63. Sketch of a MOSFET. Holes (or electrons) are trapped at the interface of 
the semi-conductor and the insulator due to the band gap difference between them, 
the dipole layer and the applied electric field. Two-dimensional electrons (holes) 
may also be found by layered structures (heterostructures) of semi-conductors with 
different band gaps such as GaAs and AlGaAs. 



?. i The two-dimensional Electron Gas 



We consider an electron gas with a uniform positive background with no com- 
plications arising from the lattice structure — this is how the many-electron 
problem was originally formulated: the Jellium model. This situation is indeed 
realized experimentally in MOSFETS (and heterostructures) in which an in- 
sulator is typically sandwiched between a metallic plate and a semiconductor 
— see Fig. 63. By applying an electric field a two-dimensional charge layer 
accumulates on the surface of the semiconductor adjacent to the insulator, 
whose density can simply be changed by varying the field strength (For de- 
tails see [26]).. Similar geometries have been used to observe the Quantum 
Hall Effects and the metal-insulator transition by varying the densityp^. 

Typically we will be interested in phenomena when the average inter-electron 
distance is (9(10) — (9(10^) nm. The thickness of the insulating layer is typically 
more than 100 nm, so that the positive (capactive) charge on the insulator 
provides a uniform background to a first approximation. In Si samles, surface 
roughness is the principle source of disorder at high densities, while at low 
densities (in the regime where the transition takes place) ionized impurity 
scattering dominates due to the fact that there is much less screening. In 
GaAs, remote impurity scattering dominates, and this scattering is mostly 
small-angle. This is the main reason the mobility in these samples is large. 



Recently, novel "high temperature superconductivity" was observed in such a 
FET geometry. Gate induced doping in such a Field Effect Transistor was employed 
to introduce significant hole densities into Gqq which became superconducting at 
Tc = 52KI [214]. 
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Neglecting disorder, the problem is characterized by r^, defined as the ratio of 
the potential energy to the kinetic energy: 



Here n is the electron density, m the band mass, and e the background static 
dielectric constant. We can also write 



which expresses that is the radius of the circle whose area is equal to the 
area per conduction electron, measured in units of the Bohr radius qq. For a 
two- valley band structure, as on the 110 surface of Si, the kinetic energy is 
reduced and is twice that defined by Eqs. (206) and (207). 

For ^ 1, (the dense electron limit) the kinetic energy dominates and metal- 
lic behavior is expected. For ^ 1, the potential energy dominates and a 
crystalline state (Wigner crystal) is expected. The best current numerical esti- 
mates put the transition to the crystalhne state at Vgc ~ 37 [238]. The entropy 
at the transition is tiny, indicating that the radial distribution function for 
the liquid state at low densities is similar to that of the crystal for distances 
up to a few times TsOo. 

It is important to note for our subsequent study that magnetism is always 
lurking close by. Reliable numerical calculations show that the magnetic state 
in the Wigner crystal near the critical density is determined by multiple- 
particle exchanges [46]. On the metallic side the energy of the ferromagnetic 
state is only a few percent above the unpolarized metallic or crystalline states 
for Ts ~ Tsc [238]. Disorder is expected to make the metal- insulator state 
continuous. On the insulating side at T — > 0, the disordered Wigner crystal 
is expected to be glassy and have low-energy properties of a Coulomb glass 
[224]. On the metallic side fluctuations in the local density of electrons might 
be expected to lead to locally polarized magnetic states or possibly to some 
unusual frustrated magnetic states [56]. The perturbative calculations with 
disorder and interactions, already alluded to [90-92,53] also hint at formation 
of magnetic moments in the metallic state. It is the interplay of such magnetic 
fluctuations with itinerancy which is one of the principal theoretical problem 
in understanding the metallic state. 




(206) 




(207) 
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8.2 Noninteracting Disordered Electrons: Scaling Theory of Localization 



Detailed reviews on the material in this section may be found in [241,149] 
[14,125]. 

The concept of localization of non-interacting electrons for strong enough dis- 
order was invented in 1958 by Anderson [18]. In one dimension all electronic 
states are localized for arbitrarily small disorder while in three dimension a 
critical value of disorder is required. That two dimensions is the marginal 
dimension in the problem was discovered through the scaling theory of local- 
ization. 

The conceptual foundations for the scaling theory of localization were laid by 
Thouless and co-workers [240,241] and by Abrahams et al. [1], and were de- 
veloped formally by Wegner [270]. Abrahams et al. [1] also made predictions 
which could be tested experimentally. Thouless noted first of all that the con- 
ductance G of a hypercube of volume in any dimension d is dimensionless 
when expressed in units of (e^//i), thus defining a scale-independent quantity 

g = G/ {e^h) . (208) 



Next, he argued that g for a box of linear size 2L may be obtained from the 
properties of a box of size L and the connection between two of them. The 
conductance of a box of size L itself increases with the transition amplitude t 
between energy levels in the two boxes and decreases with the characteristic 
width of the distribution of the energy levels in the boxes AW{L) due to the 
disorder. 



For weak-Gaussian disorder, the bandwidth may be expected to be propor- 
tional to the square-root of the number of impurities in the box, so AW{L) ~ 
The transition amplitude t is obtained by the hopping between near- 
neighbors near the surface of the boxes of size L. It is therefore proportional 
to the surface area L'^"^. Thusp^ 

g{L) = f . (210) 



Now, in three dimensions the conductivity should approach a constant for large 
L (Ohm's law!), and hence the conductance should scale as L. This implies 

This line of reasoning of course breaks down when we include electron-electron 
interactions 
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that the scahng function f{x) should go for large L as f{x) ~ x^. Note that 
while for d > 2 the g therefore increases with increasing L while for d < 2 the 
large L behavior is determined by the small argument behavior of the scaling 
function, while = 2 is the marginal dimension. 

In a very influential paper, Abrahams et al. [1] analyzed the /^-function of the 
RG flow 

(3(g) = d(\ng)/d(\nL) , (211) 

and showed by a perturbative calculation in 1/g that 

l3(g)^(d-2)-^ - , (212) 
TT^ g 

where the first part comes from Eq. (210) with f(x) ~ x'^. 

For small enough g (i.e. for large disorder) we expect exponential localization 
g{L) ~ e^^/'', where ^ is the localization length, so that P{g) ~ {—L/Q. The 
smooth connection between the perturbative result (212) for large g and the 
exponentially localized solution at small g is shown in Fig. 64. While for d = 3 
(or any d > 2), a critical disorder gc is required for localization, va.d — 2 states 
are asymptotically localized for any disorder for non-interacting fermions. The 
characteristic value of the localization length in o? = 2 is estimated from the 
perturbative solution: 

9{L) -9o-^,\n (^) , (213) 

where g^, is the dimensionless conductance dX L k, ^. In conventional Boltz- 
mann transport theory go = {e^ /2'n-h)kFt C is of the order of the value of L 
at which the correction term is of order g^, so that 

Cf«£exp . (214) 

At T — > 0, the sample size of a sample with kpi ^ 1 has to be very large 
indeed for weak localization to be observable. 

The theory described above must be modified at finite temperatures due to 
inelastic scattering. If the inelastic scattering rate is much less than the elastic 
scattering rate, r^"^ -C r~^, localization effects are cut off at a length scale 
LTh{T), the Thouless length scale: 

Lth = {Dnnf^ (215) 
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Fig. 64. The scaling function for non-interacting electrons with disorder deduced by 
Abrahams et al. [1]. 

where D = (vpr/d) is the (Boltzmann) diffusion constant. However as noted 
by Altshuler, Aronov and Khmelnitslci [12,13], the correct scale for the cut-off 
is T^^, the phase breaking rate. In an individual collision the energy change AE 
may be such that the phase changes only by a very small amount, TinAE <C 
27r. The phase breaking time is then longer and is shown to be given by 
~ {AE Tin)-'^'^ Tin. The T = theory with 



replacing L then gives the finite temperature scaling behavior to which exper- 
iments may be compared. 

The characteristic temperature at which weak-localization effects become 
prominent may be estimated in a manner similar to (213), 



This expression puts useful bounds on the temperatures required to observe 
weak localization. 

Eq. (212) is derived microscopically by considering repeated backward scat- 
tering between impurities. It can also be derived by considering quantum 
interference between different paths to go from one point A to another B [40] . 




(216) 




(217) 
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Fig. 65. Interfering (time-reversed) parts in elastic scattering off a fixed set of im- 
purities. The probability for the particle to arrive at B is reduced because of the 
enhance probability for the particle to arrive back at A, as a result of interference. 

The total probability O, for this process is 



where aj is the amplitude of the i-th path. The second term in Eq. (218) is 
non-zero only for classical trajectories which cross, for example at the point 
O in Fig. 65. The probability of finding a particle at the point O is increased 
from 2|aip to 



because the two paths are mutually time-reversed. Increasing this probability 
of course leads to a decrease in the probability of the particle to arrive at B, 
and hence to a decrease in the conductivity. 

This argument makes clear why the interfering paths must be shorter than the 

phase relaxation rate due to inelastic processes and why magnetic impurities 
or a magnetic-field which introduces phase-shift between two otherwise time- 
reversed paths suppress weak-localization. In two- dimensions [14] 



2 



(218) 



oil^ + |a2|^ + 2 Re 0^02 — 4|ai 



(219) 



a{H,T)-a{0,T) 




(220) 



where tjj is the digamma function and 



X = 4LJ eH/hc = {L^/Lnf . 
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The quantity in brackets in (220)is equal to /2A for x — > and to ln(x/4) — 7 
for X — > 00. 

Spin-orbit scattering preserves time-reversal symmetry but spin is no longer 
a good quantum number. The spins are rotated in opposite directions in the 

two self-intersecting paths of Fig. (65) if the impurities are spin-orbit scatters 
[14,149]. This has been shown to lead to an average overlap of the spin-function 
of — I (because a rotation by 2ti of wave function of a spin 1/2 particle leads 
to a wave function of opposite sign) . The correction to the /3-function of Eq. 
(212) due to this effect is 



This effect tends to enhance the conductivity. 



8.3 Interactions in Disordered Electrons 



Fermi-liquid theory for interacting electrons survives in three dimensions in 
the presence of a dilute concentration of impurities [42]. Some noteworthy 
differences from the pure case are: 

(1) Because of the lack of momentum conservation, the concept of a Fermi- 
surface in momentum space is lost but it is preserved in energy space, i.e., 
a discontinuity in particle occupation as a function of energy occurs at 
the chemical potential. Momentum of particles may be defined after im- 
purity averaging. General techniques for calculating impurity-averaged 
quantities are well developed; sec for example [4,42]. Here and subse- 
quently in this chapter the self-energies, vertices etc. refer to their form 
after impurity averaging. 

(2) In the presence of impurities, the density-density correlation (and spin- 
density correlation, if spin is conserved) at low frequencies and small mo- 
mentum must have a diffusive form (This is required by particle-number 
conservation and the continuity equation) 

7r(g, u)^K , ? < and a; < r"^ . (223) 

loj -\- Dq'^ 

Here k — dn/d/i is the compressibility and D is the diffusion constant. 

For non-interacting electrons D = ^VpT. Interactions renormalizc D and 
K [42] . In the diagrammatic representation used below, the diffusive prop- 
agator is shown by a cross-hatched line connecting a particle and a hole 
fine as in Fig. 66. 
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(3) Because of statement 1, the impurity-averaged single-particle spectral 
function at a fixed k is spread out over an energy F, so that for frequencies 
within a range F of the chemical potential, it has both a hole-part (for 
cu < /i) and a particle part (for cu > fj). This is an important technical 
point in microscopic calculations. 

(4) The Ward-identities relating the coupling of vertices to external pertur- 
bations change for the coupling to unconserved quantities (For the pure 
case they are given in section 2.6). For example no Ward-identity can 
be derived for the vertex needed for the conductivity calculation, i.e. 
Lim^^o LirUq^o A^"***"^^, because current is not conserved. 

Lim,_o Lim^-.o A^"''^ = - - (224) 

holds because after impurity averaging momentum is conserved. But mi- 
croscopic calculations show that at least when Fermi-Liquid theory is 
valid (Cf. section 2.4) 

Lim^_o Lim,_o AjT^"''^ = Lim,_o Lim^^o AL'"''"'' ■ (225) 

Indeed, if this were not so, one would not get a finite d.c. conductivity 
at T = for a disordered metal m. d — ?>. An argument for this is 
as follows: Normally, we calculate the conductivity by first taking the 
limit g — * and then the limit c<j ^ 0, as on the left side of (225). In 
practice, however, even when we apply a homogeneous field to a system, 
the electrons in a disordered medium experience a field which varies on 
the scale of the distance between the impurities, and so the physically 
relevant limit is the one on the right hand side of (225), where the limit 
a; — is taken first. But the validity of (225) appears not to extend to 
the case of Singular Fermi Liquids, at least for the present case where 
the singularities are g-dependent. This is one of the important difficulties 
in developing a consistent theory for disordered interacting electrons in 
d = 2. 

The diffusive form of the density correlation function and spin-density corre- 
lation is the genesis of the singularities due to interactions in two-dimensions. 
For example the elementary effective vertex in Fig. 66 due to a bare frequency- 
independent short-range interactions in two-dimensions is 

£-1 

v'^ I dqq I ) ^ N(0) ln(ujT) . (226) 



The singularity arises because 7r(u;,g) = f{uj/Dq^). Recall that for pure elec- 
trons 7r(a;,g) = f{ui/vpq) leading to a logarithmic singularity for the second- 
order vertex in one-dimension and regular behavior in higher dimensions. Sim- 
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Fig. 66. Elementary processes important in the problem of two-dimensional disor- 
dered interacting electrons and referred to in the text, (a) Representation of the dif- 
fusion propagator due to impurity scattering vertices and corresponding self-energy. 
The particle lines and hole lines should be on opposite sides of the chemical potential. 
(b) Singular second-order interactions, (c) Singular vertex in the density channel 
(and in the spin-density channel for the spin-conserving problem), (d) Singular (ir- 
reducible) first-order interactions, (e) Elementary singular polarization propagator. 

ilarly 7t{uj, q) = f{uj/q^) leads to a logarithmic singularity in the second-order 
vertex in three-dimensions, as we saw in section 5.1 on SFL's due to Gauge 
interactions. 

Note that in Eq. (226) and other singular integrals in the problem have ultra- 
violet cutoffs at g ~ and a; ~ since the diffusive form is not applicable 
at shorter length scales or time scales. It also follows that Boltzmann transport 
theory is valid at temperatures larger than t~^. 

Actually even the first order interaction dressed by diffusion fluctuations is 
singular. Consider first the diffusion correction to the vertex shown in Fig. 15 

^ ^- {iu + Dq'Y' , (227) 

provided e < 0,e — uj > Oor vice- versa. The restriction is a manifestation 
of point (3) and arises because in the diffusion process, only intermediate 
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(a) 



(b) 

{q,uj), S = or 1 

I 
I 

Fig. 67. Effective interactions can be split up into singlet and triplet channels. 
In the singlet-only channel (a), the density-density interaction is screened by the 
Coulomb interaction and is universal at long wavelengths. In the triplet channel and 
in the singlet channel for large momentum, the screened density-density interaction 
appears only in the cross channel and is therefore non-universal. 

states with one line above (particle) and the other below (hole) the chemical 
potential contribute as they alone define the physical density. This leads to 
the first order irreducible interaction and the polarization graph shown in Fig. 
67 to be logarithmically singular. 

For the small q of interest for singular properties, one need consider interac- 
tions only in the s-wave channel. One then has two interaction parameters, 
one in the singlet channel and the other in the triplet channel. 

Consider the problem with Coulomb interactions. Then the effective inter- 
action in the singlet channel sums the polarization bubbles connected by 
Coulomb interactions. Using (223) for the polarization bubble, it is shown [14] 
that for small momentum transfer the interaction in the spin singlet {S = 0) 
channel. Fig. 67, becomes 

^singlet = 2« . (228) 

In the non-interacting limit k — N{0), independent of density. Consider now 

the ladder-type interactions illustrated in Fig. 68. These involve both the 
singlet and the triplet interactions. The momentum carried by the interaction 
lines is however to be integrated over. Therefore the triplet-interactions do not 
have a universal behavior, unlike the singlet interactions. 

Altshuler, Aronov and collaborators [14] calculated the logarithmic correc- 
tions to first-order in the interactions for various physical quantities. To these 
one can add the contribution already discussed due to weak-localization. The 
corrections to the single-particle density of states, the specific heat and the 
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(a) 





Fig. 68. Simplest processes contributing to the singular self-energy, (a) Exchange 
process, (h) Hartree process. 



conductivity over the non-interacting values are respectively: 
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The compressibility has no logarithmic correction. In these equations, s is the 
screening length and F is a parameter which is of the order of the dimensionless 
interaction r^. 



The first terms in (230) and (231) are due to exchange processes and the second 
due to the Hartree processes. The exchange process, of which the contribution 
to the self-energy is shown in Fig. 66, use the interaction in the singlet channel; 
hence the universal coefficient. The second contribution uses both triplet and 
the (large q part of the) singlet interactions. In first-order of interaction the 
difference in sign of the two processes is natural. In pure systems, the Hartree 
process does not appear as it involves the g = interaction alone which is 
exactly canceled by the positive background. For disordered systems, due to 
the fluctuation in the (ground state) density, a first order Hartree process. Fig. 
68 contributes. 



In the presence of a magnetic field the = ±1 parts of the triplet interac- 
tions acquire a low-energy cut-off. Therefore the logarithmic correction to the 
resistivity is suppressed leading to negative magneto-resistance proportion to 
F{H/kTY for small H/kT but qhbH » t~^^,t^^ where and t~'^ are 
spin-orbit and spin-scattering rates respectively, for appropriate impurities. 
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8-4 Finkelstein Theory 



Finkelstein [90] has used field-theoretical methods to generalize Eqs. (229)- 
(231) beyond the Hartree-Fock approximation. His results have been rederived 
in customary diagrammatic theory [51,39]. The interference processes leading 
to weak localization are again neglected. The theory may be regarded as first 
order in l/kpi- In effect the method consists in replacing the parameter F by 
a scattering amplitude 7t for which scaling equations are derived. The equiva- 
lent of the parameter is fixed by imposing that the compressibility remains 
unrenormalized, i.e., does not acquire logarithmic corrections. A second im- 
portant quantity is a scaling variable z, which is analogous to the dynamical 
scaling exponent z which we discussed in the section 6, which gives the rela- 
tive scaling of temperature (or frequency) with respect to the length scale. A 
very unusual feature of the theory is that z itself scales! Scaling equations are 
derived for ■jt and z to leading order l/kpi. As T — 0, both 7t and z diverge. 
The divergence in z (see the discussion in section 7) usually means that the 
the momentum dependence of the fluctuations is unimportant compared to 
their frequency dependence. The divergence in 7^ as T — > in such a case has 
been interpreted to imply divergent spatially localized magnetic-fluctuations; 
in other words, it implies the formation of local moments [91,53]. At the same 
time, the scaling equations show conductivity flowing to a finite value. 

The scaling trajectories of Finkelstein's theory are shown schematically in 
Fig. 69. While in the non-interacting theory with disorder, one always has an 
insulator, this theory flows always towards a metal. However the theory cannot 
be trusted beyond 7^ ~ 1 beyond which it is uncontrolled. The theory also 
cannot be trusted for large disorder, kpi ~ C'(l) , even for small interactions. 

It is worth emphasizing that Finkelstein's theory gives an effect of the inter- 
actions in a direction opposite to the leading perturbative results. The pertur- 
bative results themselves of course are valid only for small while Finkelstein 
theory is strictly valid only for < 0{1). One possibility is that the Finkel- 
stein result itself is a transient and the correct theory scales back towards an 
insulator. Another possibility is that it correctly indicates (at least for some 
range of and disorder) a strong-coupling singular Fermi-liquid metallic fixed 
line. The new experiments discussed below can be argued to point to the latter 
direction. 

It is hard except in very simple situations (the Kondo-problem, for instance) to 
obtain the approach to a strong-coupling fixed point analytically. In that case 
one may usually guess the nature of the fixed point and make an expansion 
about it to ascertain its stabilityP^. 

Although the theory breaks down in the strong coupling regime, this situation is 
somewhat comparable to the hints that the weak coupling expansion gave in the 
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Fig. 69. Schematic renormalization group flow for the disordered interacting electron 
problem according to the Finkelstein theory, dashed lines represent the effect on the 
solid lines on applying a magnetic- field which couples to spins alone. 

In making such a guess, the SFL properties towards which the Finkelstein 
solution flows should be kept in mind: 



(i) The conductivity flows towards a flnite value in the theory as T 

(ii) The density of single-particle states flows towards zero 



0. 



N{uj) 



(232) 



(Hi) The magnetic susceptibility diverges at a flnite length scale (the effect of 
a diverging z) indicating the formation of local moments. 

The last point appears to be crucial. As may be seen from Eq. (239) below, 
the growth of the triplet scattering overrides the exchange processes which 
favor the insulating state. Indeed if the triplet divergence is suppressed by an 
applied magnetic fleld, the theory reverts to the pcrturbativc form of Eq. (220). 
The scale of the magnetic-fleld for this effect is given by the temperature. The 

early phase of the work on the Kondo problem: these weak coupling expansions 
broke down at temperatures comparable to the Kondo temperature, but did hint at 
the fact that the low temperature regime was a strong coupling regime. 
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formation of localized regions of moments may be linked to the fact already 
discussed that the ferromagnetic state is close in energy to the paramagnetic 
fluid (and the crystalline states) as density is decreased. The experiments 
discussed below have a significant correspondence with this picture, although 
there are some crucial differences. 

A possible strong coupling fixed point is a state in which the local moments 
form a singlet state with a finite spin stiffness of energy of 0{Hc). This elimi- 
nates any perturbative instability of the triplet channel about the fixed point. 
The state is assumed to have zero density of single-particle states at the chem- 
ical potential. This eliminates the localization singularity as well as well as the 
singularity due to the singlet channel. This state is then perturbatively stable. 
The issue of whether such a state has finite conductivity is unresolved. The 
occurrence of a new characteristic scale He observed in the magneto-resistance 
experiments discussed below with He ^ as the metal-insulator transition in 
zero field as n —> ric is also in correspondence with these ideas. 



8.5 Compressibility, Screening length and a Mechanism for Metal- Insulator 
Transition 

Suppose the metallic state in two dimensions is described by a fixed point 
hinted by the Finkelstein theory and expansion about it. Such a description 
must break down near the critical where a first order transition to the 
Wigner transition must occur in the limit of zero disorder. General arguments 
suggest that the transition for finite disorder must be continuous [124]. 

A suggestion for the breakdown of the Finkelstein regime follows from the 
calculation of the correction to the compressibility due to disorder [226]. As 
already mentioned, no perturbative singularity is found in the compressibil- 
ity due to interactions. However, the correlation energy contribution of the 
zero-point fluctuations of plasmons is altered due to disorder with a mag- 
nitude which depends also on r^. The leading order contribution in powers 
of {kpt}'^ can be calculated for arbitrary Vg- Including this contribution the 
compressibility k may be written in the form 

+ O.Url/M + O {{rt)/{u;orr) . (233) 

K K pure 

Here Hpure is the compressibility for zero disorder, kq = iV(0), and ljq is the 
Rydberg. In the Hartree-Fock approximation 

^ = 1 - {V2/n)rs (234) 

K pure 
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The best available numerical calculations also give kq/k varying slowly enough 
with Tg that the correction term (233) dominates for r.^ of interest near the 
metal-insulator transition even for modest disorder; for example for uuqT ~ 10, 
the disorder contribution in Eq. (233) is larger than the pure contribution 
for Vs > 10. This has an important bearing on the Metal-Insulator transition 
because the screening length s is given by 

s/so = Ko/k, (235) 

where Sq = ao/2. Strictly speaking s is the screening length for an exter- 
nal immobile charge and the screening of the electron-electron interactions is 
modified from (233) due to vertex renormalizations. But in this case they do 
not change the essential results. From Eq. (233) it follows that the screening 
length s{£) > i, the mean-free path, for 

Ts > 3(a;or) . (236) 

Suppose the condition s{£) > L > i is satisfied. Here L again is the size of the 
box for which the calculation is done, defined through DL~'^ T. The assump- 
tion of screened short-range interactions, with which perturbative corrections 
leading to results of Eqs. (229), (230), and (231) are obtained, is no longer 
valid. In this regime the calculations must be done with unscreened Coulomb 
interactions. The correction proportional to F in Eqs. (229), (230), and (231) 
is not modified but the singlet contributions are more singular (due to the 
extra in the momentum integrals). For instance, Eq. (231) is modified to 

5a/a ^-iV2/7rysj ■ (237) 

This imphes a crossover to strong localization. It is therefore suggested that 
the metaUic state ceases to exist when condition (236) is satisfied. 

The above line of reasoning is of particular interest because as discussed below 
a sharp variation in the compressibility is indeed observed to accompany the 
transition from the metallic-like to insulating-like state as density is decreased 
(as shown in Fig. 79 below). 

8. 6 Experiments 

Soon after the publication of the theory of weak localization its predictions 
were seemingly verified in experiments on 5'i-MOSFETS [47,245]. The ex- 
periments measured resistivity on not very clean samples of high density with 
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Fig. 70. Resistivity data on a fine scale for the two highest densities in Fig. 75 below, 
showing correspondence with the theory of weak localization at such high densities. 
From Pudalov et al. [202]. 

resistivity of O(10~^/i/e^). In a limited range of temperature the predicted log- 
arithmic rise in resistivity with decreasing temperature with about the right 
prefactor was found [40]. In view of the perturbative results of Altshuler and 
Aronov [14] and the knowledge that electron-electron interactions alone lead 
to a Wigner insulator at low- densities, one w^as led to the conviction that the 
metallic state does not exist in two dimensions. It was expected that sam- 
ples will larger will simply show logarithmic corrections to the resistivity 
at a higher temperature and pure samples at a lower temperature. Not too 
much attention was paid to Finkelstein's results which pointed to the more 
interesting possibility of corrections in the opposite direction. 

The more recent experiments on a variety of samples on a wider range of den- 
sity and of higher purity than earlier have re- focused attention on the problem 
of disorder and interactions in two dimensions and, by implication, in three di- 
mensions as well. Several reviews of the experiments are available [3,16,17]. We 
will present only a few experimental data to highlight the theoretical problems 
posed, and will focus on the behavior of the data as function of temperature. 
The scahng of the data as a function of the electron density n — ric or field 
E will not be discussed; there is a considerable body of data on nonlinear 
^'-dependence (see e.g. [223] and references therein) but the significance of 
the data is not clear at present. 

The first thing to note is that results consistent with the earlier data [47,245] 
are indeed obtained for high enough densities. Fig. 70 shows the resistance 

versus temperature in Si for ~ The magnitude of the temperature 

dependence is consistent with the predictions of weak localization corrections. 
As we will show below in the same region of densities, the negative magneto- 
resistance predicted as correction to weak-localization, discussed above, is also 
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Fig. 71. Resisitivity as a function of temperature for a wide range of densities (and 
Fermi-energy) in a disordered Si MOSFET. The inset shows accurate measurements 
of p{T) close to the speparatrix for another sample. From Sarachik and Kravchenko 
[211,3]. 

observed. 

Fig. 71 shows the resistivity as a function of T over a wide range of densities. 
Similar data from [142] over a large small of densities is shown in Fig. 72, 
and data over a large range of densities are plotted as a function of T/Ep 
in Fig. 73. The resistivity clearly shows a change of sign in the curvature as 
a function of density at low temperatures. The resistivity at the cross-over 
density as a function of temperature is shown down to 20mK in the inset to 
Fig. 72 and is consistent with temperature independence. The true electron 
temperature in these samples is a question of some controversy [16,3], but more 
recent experiments, whose data are shown in Fig. 74, have corroborated these 
results by studying this issue very carefully down to 5mK. In the high density 
region the resistivity does rise with decreasing temperature logarithmically, 
consistent with earlier measurements. The consistency of these data sets for 
two very different types of samples therefore gives strong evidence that these 
are genuine effects in both types of systems. 

The data shown in Figs. 71-73 is for 5i-M0SFET samples. The data for 
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Fig. 72. Resistivity versus temperature at five different electron densities in the 
experiments of Kravchenko and Klapwijk [142]. The inset shows that the middle 
curve {ug = 7.25 x 10-*^° cm"^) changes by less than ±5% in the entire temperature 
range. 

GaAs heterostructures, and Si in other geometries is qualitatively similar 
[115,61,62,196,114,172]. Fig. 74 shows data on high quality gated GaAs quan- 
tum wells with densities on the metallic side of the metal-insulator "transition" 

taken to temperatures as low as 5mK. The resistivity is essentially temper- 
ature independent at low temperatures; the logarithmic corrections expected 
from weak-localization (calculated using the measured resistivity and the the- 
oretically expected r,/,) is also shown. 

Experiments in a Parallel Magnetic Field: A magnetic field applied parallel 
to the plane couples primarily to the spin of the electrons. For small-fields 
and for n 3> ric, a positive magneto-resistance proportional to is observed 
as expected from perturbative calculations in the interactions. For fields such 
that /ibH ~ Ep, the electrons are fully polarized and the resistivity satu- 
rates as expected. The temperature dependence of the resistivity begins to be 
insulating-like at low temperatures with the cross-over temperature increas- 
ing as n decreases [228]. This is an indication that the metallic state becomes 
unstable as the spins are polarized. A complete set of data is shown in Fig. 75 
where resistivity versus temperature in a S'z-MOSFET with density varying 
across is shown together with the resistivity versus magnetic-field at the 
lowest temperature for some densities on the n > Uc side. It is noteworthy 
that the temperature dependence of the high-field data (not shown) appears 
to fall on the curve of resistivity versus temperature (at H=0) which the high 
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Fig. 73. Plot of the resistivity as a function of the scaled temperature T/Ep. The 
encircled region indicates the range of parameters explored in Fig. 72 and in [142]. 
The dash-dotted vertical line depicts the empirical temperature Tq = 0.007-Ei? 
below which the logarithmic temperature dependence like that of Fig. 70 sets in. 
Prom Prinz et al. [200]. 

field (low temperature) data saturates asymptotically. 

The parallel magneto-resistance has been examined carefully for n close to but 
larger than Uc, and is shown in Figs. 76 for p-type GaAs [279]. Qualitatively 
similar results are found in n-type Si [203]. It is discovered that a critical field 
as a function of density Hdji) exists such that for H < Hdji) the resistivity 
continues to be metallic-like dp/dT > and for H > Hc{n) it is insulating like 
dp/dT < 0. The field Hc{ri) tends to zero as n —> ric. The low temperature 
data on the high fields side is puzzling and should be re-examined to ensure 
that the electron-temperature is indeed the indicated temperature. 

Experiments in a Perpendicular Field: The behavior of a resistance in all 
but very small perpendicular fields, is dominated by the Quantum Hall Effect 
(QHE). The connection of the Quantum Hall transitions to the metal-insulator 
transition aX n = and iJ = is an interesting question which we will not 
touch on. At low fields and for n > Uc, outside the QHE regime, negative 
magneto-resistance predicted by weak-localization theory are observed. Data 
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Fig. 74. (a) Temperature dependence of the longitudinal resistivity of a 2D hole 
gas for various gate biases and associated densities in experiments on GaAs. The 
solid curves are estimated weak localization predictions, (h)-(f): Magnified view of 
the data \n (a) averaged over a 5% temperature interval. The estimated weak local- 
ization prediction has been shifted to coincide with the data curve at T = 50mi^. 
Prom Mills et al. [172] 

for n ^ ric is shown in Fig. 77 and agrees quite well with the theoretical 
curves as shown; similar results for n close to ric are also reported [114]. More 
recent low temperature data in GaAs hetcrostructures[172] is reproduced in 
Fig. 78 for densities n > ric but close to ric- A magnetoresistance two orders 
of magnitude smaller than the weak localization theory is estimated although 
the width of the negative magneto-resistance region is not in disaccord with 
the weak-localization correction. 



Compressibility measurements Compressibility (/t) measurements [79,122] in 
the region around n = Uc show a rapid change in from the negative value 
characteristic of high Vg metallic state to positive values — see Fig. 79. These 
are very important measurements which show that a thermodynamic quantity 
has a very rapid variation near n ~ Uc- We have already discussed that such 
changes were predicted [226] to occur through perturbative corrections due 
to disorder to the energy of interacting electrons. Some more recent measure- 
ments [276] of the local electrostatic potential show that in the region n ^ ric 
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Fig. 75. Results for resistivity versus temperature and versus magnetic field applied 
in the plane for a few densities on either side of ric ■ The magnetic field is shown 
on the upper axis and the data is taken at the lowest temperature for some of the 
densities shown in the resistivity versus temperature plots. From Pudalov et al. 
[201]. 

large scale density fluctuations (puddles) occur with very weak connections 
between them. Completely isolated puddles evidently have k.^^ = 0. 

8. 7 Discussion of the Experiments in light of the theory of Interacting Dis- 
ordered Electrons 

In comparing the experimental results with the theory, it is necessary to sepa- 
rate out the effects due to "customary-physics" — for instance electron-phonon 
interactions, creation of ionized impurities with temperature [15], change of 
screening from its quantum to its classical form as a function of temperature 
[72] change of single particle wavefunctions with a magnetic fields [137], inter- 
valley scattering [277] , etc. — from the singular effects due to impurities and 
interactions. The separation is at present a matter of some debate. However 
it seems that the following features of the experimental data in relation to the 
theoretical ideas summarized in sections 8.3, 8.4 are especially noteworthy. 
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Fig. 76. Plot of the magnetoresistancc. In (a) the T dependence of p in the zero 
field metallic phase is shown on a semilog plot for a hole density 3.7 x 10 



10 



cm for 

varying B\\ values. As S|| increases from zero, the strength of the metallic behavior 
measured by the total change in p from about IK to 50mK weakens progressively, 
and for _B[| > dp/dT becomes negative (i.e. the system becomes insulating). 
An alternate way of demonstrating the existence of a well defined -Bjj is to plot p 
against B\\ at several different temperatures. In. (h) p is plotted vs B\\ at a hole 
density 1.5 x 10^^cm~'^. 5|j is read off the crossing point marked by the arrow. In 
(c), the differential resistivity dV/dl is measured at 50mK across the B\\ induced 
Metal-Insulator transition is shown at magnetic field strengths similar to those in 
(a). From Yoon et al. [279]. 
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Magnetic Field (Gauds') 

Fig. 77. Plot of the magnetoresistance. The ehange in the resistivity 
^p{B) = pxx{B) — Pxx{^) vs magnetic field B at an electron density of 1.05 x 
10^^ cm~^ at various temperatures. The open circles denote the measurements and 
the full line is the best least square fit according to the single electron weak local- 
ization correction to the conductivity. From Brunthaler et al. [48]. 
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Fig. 78. Variation of the longitudinal resistance with perpendicular magnetic field 
for 2D sample at T = 9mK and at various indicated densities. The weak-localization 
correction is estimated to be O(IO^) larger than the observations at these densities. 
From Mills et al. [172]. 
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Fig. 79. Compressibility data. In this experiment, at low frequencies, Ix is directly 
proportional to Rx, the dissipation of the two dimensional hole system, while ly 
is proportional to the inverse compressibility. Ix and ly are shown as a function 
of density for five different temperatures ranging from 0.33 K to 1.28 K at an 
excitation frequency of lOOHz. The crossing point of the five dissipation channel 
curves corresponds to the metal-insulator transition at B = 0. The minimum of 
the inverse compressibility occurs at the same hole density of 5.5 x 10^^ cm~^. From 
Dultz et al. [79]. 



These must be read bearing in mind our earlier discussion that most of the 
interesting experiments are in a range of and disorder where the theoretical 
problems are unresolved and only hints about the correct form of a theory are 
available. 



• At < 0(1) and /ci^f ^ 1 a logarithmic increase in resistance with decreas- 
ing T consistent with weak-localization as well as with the perturbative inter- 
action correction is observed. A positive magneto-resistance consistent with 
the latter is also observed. Also observed is the correction to weak-localization 
due to phase-breaking of backscattering in a perpendicular magnetic field. The 
latter yields sensible values and temperature dependence for the phase relax- 
ation rate given by the theory. It appears that at high enough density, the weak 
localization theory supplemented by the perturbative theory of interactions is 
in excellent agreement with the experiments in the range of temperatures 
examined. 



• As is increased (and kpl decreased) the logarithmic resistance is lost in the 
observed temperature range, whereas weak-localization theory predicts that 
the coefficient of such terms (as well as the onset temperature for their oc- 
currence) should increase. For not too large the decreased logarithmic term 
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may be associated with the perturbative corrections (220) due to interactions. 

• Upon further increasing r^, the derivative dp/dT becomes positive in the low 
temperature region as in a metal. The magneto-resistance in a parallel field is 
positive ~ as is predicted by Finkclstein [although the variation is closer 
to H^/T rather than as {H/TY] [203]. The phase breaking correction in a 
perpendicular field continues to be observed. But quite curiously the deduced 

is larger than r deduced from resistivity — by definition a phase-breaking 
rate serves as a cutoff only if < r. 

• In the "metallic" regime for intermediate r^, a strongly temperature de- 
pendent contribution for T < Ep is found which may be fitted to the form 
p'(n) exp(— ii^a(n)/T). The magnitude of this term rapidly decreases as the 
density n decreases. No accepted explanation for this contribution has been 
given. In Si the change of resistivity dX n lOric due to this contribution 
is an order of magnitude larger than in GaAs. It has been proposed [14,201] 
that this contribution together with the weak-localization contribution may 
well account for all the data in the "metallic" regime since it pushes the min- 
imum of the resistivity below which the logarithmic temperature dependence 
is visible to lower temperature than the available data at lower densities. 

This issue can be resolved by experiments at lower temperatures. At this point, 
especially in view of the consistency of the recent results of Mills et al. — see 
Figs. 74 and 78 — with the earlier experiments, one can say that it requires an 
unlikely conspirancy of contributions to remove the temperature dependence 
over a wide range for different materials and with different degrees of disorder. 

A quite different scenario also consistent with the existing data is that the 
logarithmic upturn in the resistivity observed in high density samples is a 
transient that on further decreasing the temperature disappears to be replaced 
by ^ — > as T — > 0, at least above some characterestic density which is a 
function of diorder. We will come back to this issue when we discuss the 
possible phase-diagram. 

• As Vg approaches Tsc, dp/dT tends to zero (through positive values). A 
separatrix is observed with dp/dT ^ over about 2 orders of magnitude in 
temperature for Si and over an order of magnitude in GaAs. For > r^c, 
dp/dT is negative befitting an insulator. Vgc appears to be smaller for dirtier 
samples but not enough systematic data is available for drawing a functional 
relation. 

• The electronic compressibility rapidly changes near the "transition" and 
rapidly becomes small on the insulating side. Its value on the insulating side 
is consistent with approaching zero in the limit of zero temperature. Although 
this is in qualitative accord with the theoretical suggestion [226], further ex- 
periments measuring simultaneously the compressibility and the conductivity 
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at low temperatures are necessary to correlate the metal-insulator transition 
with the rapid variation of compressibility or the screening length. Note that 
it follows from the Einstein relation a = Dk, that if k is finite in the metallic 
state (cr finite) and zero in the insulating state (a = 0), k must go to zero at 
the transition. Otherwise we would have the absurd conclusion that D ^ oo 
at the transition. 

Interesting phenomenological connections between the transport properties 
and formation of "puddles of electron density" of decreasing size as the metal- 
insulator "transition' is approached have been drawn [166]. 

• Near = r^c the resistivity as a function of temperature on the insulating 
side appears to be a reflection of that on the metallic side about the dp/dT = 
line if the data is not considered to low temperatures [141]. Now with more 
complete data, we know that the resistivity flattens to zero slope at low tem- 
peratures on the "metallic" side of Uc ■ The most likely behavior appears to 
be that the resistivity approaches a finite value at low temperature on one 
side of Uc and an infinite value on the other. A one-parameter scaling ansatz 
[75] for the problem with interaction and disorder gave p —>■ oo for n < Uc 
and p ^ for n > -n-c as T ^ and reflection symmetry just as at any 
second-order transition with one scaling parameteip^. Does this necessarily 
imply that mult i- dimensional scaling is required near this transition? But an- 
other important point to bear in mind is that since resistance does not depend 
on a length scale in two- dimensions, it need not be a function, in particular, 
of the correlation length near the transition. The resistivity is allowed to be 
finite on one side of a metal-insulator transition and infinite on the other even 
though the transition may be continuous and the correlation length diverges 
on either side with the same exponent. The glassy nature (Coulomb glass) of 
the insulating state is also expected to change the critical properties. 

• The resisitivity at low temperatures for n -C ric has been fitted to an acti- 
vated form oc exp(A/T)" with a ~ 1/2 and with A —> as n —> ric. This is 
characteristic of a Coulomb glass [224] . Whether this is indeed the asymptotic 
low temperature form is not completely settled. 

• The resistance a,t n k, n^. appears to vary from sample to sample but is 
within a factor of 3 of the quantum of resistance. It is worth emphasizing that 
ric is close to the density expected for Wigner crystallization. With Coulomb 
interactions and disorder, the insulating state is indeed expected to be Wigner 
glass. In that case one might expect singular frequency dependent properties 
and hysteretic behavior near the transition. 

• For Ts ~ Tsc the resistance in a parallel-field is especially noteworthy. In 

The data also led to suggestions for a superconducting ground state on the metal- 
lic side [193], and to an anyonic state [286]! 
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a parallel field dp/dT decreases untill at a field H = Hc{n) it changes sign. 
Hf. vanishes at ric, the density where dp/dT — for — 0. In this regime, 
p{H, T, n) can be scaled as [279] 

J jn- Uc) {H - Hcjn - nc)) \ 

P[ rpa ^ ] ^ ' 



This means that the transition from the metallic state to the insulating state 
can be driven by a magnetic field. It appears that the "metallic" state owes 
its existence to low-energy magnetic fluctuations which are quenched by a 
magnetic fleld. This is in line with Finkelstein theory and the flow diagram of 
Fig. 69 yet the existence of a scale He is not anticipated by the calculations 
of Finkelstein (nor, of course, is the mere existence of ric). As is further 
increased dp/dT approaches the insulating behavior characteristic of n > tt-c 
at H\\ = 0. At a flxed temperature the resistivity saturates for gpBH\\ > Ep, 
i.e. for a fully polarized band. 

For small perpendicular flelds, negative magneto-resistance of the form of (220) 
continues to be observed at least for Si for n ^ Uc- In GaAs this contribution 
at least in the range nc> n> 2nc is negligible. 

• The Hall coefficient Rh is continuous across the transition, obeying the 
kinetic theory result Rh ~ On the "metallic" side this is not surprising. 
On the "insulating" side this is reminiscent of the properties of Wigner glasses 
[60,105] 



8.8 Phase Diagram and Concluding Remarks 



It is worthwhile to try to guess the T — phase-diagram of interacting disor- 
dered electrons on the basis of the data and the available theory, inadequate 
though it is. A convenient set of axes is (or n~^/^), as it parameterizes the 
dimensionless interaction, and the resistance in units oih/e^ as it parametrizes 
the dimensionless disorder, see Fig. 80. 

Reliable theoretical results are found only along the two axes of Fig. 80. States 
are localized all along the horizontal axis. Localized states at T = must be 
organized into one or another kind of magnetically ordered state. On the 
vertical axis a Fermi-liquid gives way via a first order transition to a Wigner 
crystal which may have various magnetic phases. 

In light of the experiments, the assumption made too long that the entire re- 
gion in Fig. 80 is an insulator has to be abandoned in all likelihood. There does 
appear to be a " metaUic state" . With disorder a crossover to a Wigner glass 
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Insulator 



Impurity density 

Fig. 80. A tentative phase diagram at T = for two-dimensional disordered electrons 
with interactions. 

must occur at large r^. It is also clear that at high densities weak- localization 
theory supplemented by perturbative corrections due to interactions works 
quite well in the range of temperatures examined. At moderate for small 
disorder the Finkelstein correction appears to take over and a "metallic" state 
takes over. The best evidence for this, paradoxically, is the magnetic filed (par- 
allel to the plane) dependence of the resistivity which appears to eliminate the 
"metallic" state. 

Based on these considerations, the phase diagram Fig. 80 is put forth. It is 
surmised that the weak localization correction flows to strong localization for 
sufficiently strong disorder and small enough r^, but that it gives way to a 
metallic state at weak-disorder and larger r^. What determines the boundary? 
A possible criterion is that on one side the Finkelstein renormalization is more 
important and on the other side localization due to disorder is more important. 
The cross-over to strong localization occurs at a length scale ( given by Eq. 
(214) where the resistivity doubles. 

The scaling equation for the triplet interaction parameter is [90,50,149] 

djt/dluL = n kpi (1 + 2jtf , (239) 

so that for small initial value 7° at L = £ one gets 'jt{L) = '~f" + {'n-kFi) ln(L/7). 
As L = 7( ^ 7° + 1. The boundary between the "metallic" and the "insu- 
lating" regions on this basis is linear at small and small g~^, as shown in 
Fig. 80. This is highly conjectural but the existence of the phase boundary at 
the point ^ and 1/(7 —> is more robust. 

This scenario can be tested in high density, low disorder samples by mea- 
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surements of resistivity at very low temperatures. If correct in some regime 
of parameters near the boundary, the logarithmic weak-localization correction 
should appear at high temperatures and disappear at lower temperatures. 

We have stressed that the "metaUic state" in 2-dimensions is likely to be a 
Singular Fermi-liquid with an interesting magnetic-ground state. Direct or 
indirect measurements of the magnetic susceptibility through, for instance 
the magnetic field dependence of the compressibility should yield very inter- 
esting results. Also interesting would be measurements of the single-particle 
density of states through tunneling measurements. Further systematic and 
careful measurements of the compressibility are also required to correlate the 
transition from the "metallic" state to the increase in susceptibility. 

The basic theoretical and experimental problem remains the characterization 
of the "metallic state" its low temperature entropy, magnetic susceptibility, 
single-particle density of states etc. The experimental and theoretical problems 
are many but one hopes not insurmountable. 
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